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RESEARCH  SUMMARY 

Regression  hypotheses  admitted  to  statistical  evalua- 
tion should,  of  course,  be  germane  to  relations  between 
the  variables  considered.  This  constraint  is  easily  violated 
with  the  screening  of  arbitrarily  adopted  sets  of  hypothe- 
ses. Sequential  polynomials  in  additive  models,  for  exam- 
ple, may  include  components  that  are  not  commensurate 
with  expectation  (prior  knowledge)  and  may  not  include 
those  that  are  (e.g.,  unique  interaction  terms)  and,  since 
the  analytical  results  from  the  screening  of  meaningless 
hypotheses  may  themselves  be  meaningless,  there  is 
substantial  incentive  for  serious  preanalysis  considera- 
tion of  the  hypotheses  to  be  evaluated. 

Where  prior  knowledge  is  reasonably  strong,  the  rela- 
tion under  consideration  is  highly  interactive  and/or  has 
unique  curve  form  or  intercept  constraints,  the  develop- 
ment of  a  fairly  specific  (structured)  hypothesis  is  prob- 
ably advisable.  This  can  be  evaluated  on  pertinent  data 
sets  not  involved  in  the  hypothesis  development. 

Prior  knowledge  information  sources  can  include 
theory  resulting  from  past  investigation  as  well  as  infor- 
mation derived  from  data  sets  or  parts  of  such  data  sets 
devoted  to  hypothesis  development.  Methods  for  develop- 
ing multidimensional  graphic  models  from  two- 
dimensional  effects  found  in  partitioned  data  and  for  in- 
corporating even  complex  constraints  are  shown  in  both 
contrived-  and  live-application  examples. 

Mathematical  descriptors  for  the  resulting  graphic 
models  are  assembled  using  two  families  of  curves:  X", 
n>  1;  and  e_K,  a  highly  manipulate  function  for  bell- 
shaped  curves  or  portions  thereof.  In  either  case,  the 
X-axis  for  each  two-dimensional  effect  is  altered  as 
needed  to  orient  the  curve  in  the  upper-right  quadrant. 
The  analyst  then  simply  compares  a  normalized-scale  ver- 
sion of  the  resulting  curve  with  graphed  standard  curves, 
Standards,  and  locates  the  two  adjacent  ones  most  nearly 


like  his  own  in  shape.  Interpolation  between  the  curves 
and  their  associated  parameters,  results  in  identification 
of  the  function  for  the  best  matching  curve  offered  by  the 
system. 

The  alternative  of  using  a  completely  automated 
X-transform  selection  system  from  the  Standards  was 
considered  at  length.  But,  inability  to  define  criteria  (in- 
cluding the  fitting  process)  that  always  result  in  the  best 
choice  of  descriptor,  regardless  of  data  distribution  or 
analytical  objectives,  resulted  in  retention  of  the  less 
sophisticated,  but  still  completely  flexible  graphic 
comparison  technique  shown  herein.  Methods  of  use  are 
fairly  well  documented  for  each  of  the  graphed-relation 
examples,  which  vary  from  two  to  five  dimensions. 

The  mathematical  descriptors  can  be  regarded  as  perti- 
nent hypotheses  to  be  evaluated  on  data  sets  not  used  in 
their  development.  The  evaluation  on  a  single  data  set 
can  involve  a  simple  scale  adjustment  in  Y  =  b  (hypoth- 
eses), as  shown  herein,  or  a  variety  of  more  complex 
testing  formats  (not  shown)  depending  on  complexity  of 
the  hypothesis  itself,  capacity  of  the  analyst  to  linearize 
it,  and/or  on  the  nature  of  the  fitting  system  adopted 
(linear,  nonlinear).  Also  evaluation  can  be  based  on  per- 
formance statistics  compiled  from  real-world  applications 
of  the  hypothesized  descriptors. 

Experience  gained  from  the  presentation  of  a  dozen 
short  courses  covering  methodology  through  the  fourth 
dimension,  suggests  that  this  material  can  be  assimilated 
by  researchers  with  only  college  algebra-level  mathe- 
matical background  and  brief  training  in  the  operation  of 
programable  electronic  calculators.  Also  note  that  the 
graphic  presentation  scheme  for  the  four-  and  five- 
dimensional  relations  used  in  this  paper  can,  with 
nominal  explanation,  be  readily  interpreted  by  land 
managers  having  little  or  no  interest  in  the  complex 
mathematical  descriptors  that  may  be  involved. 
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INTRODUCTION 

Applications  of  statistics  reported  in  biological  literature  were  found  by  Mead  and  Pike 
(1975)  to  contain  some  prevalent  and  potentially  damaging  oversights.  Among  other 
things,  they  noted  that,  when  transformations  were  used,  simple  polynomials  predomi- 
nated and  seem  to  be  used  as  the  simplest  readily  available  smoothing  curve,  without  any 
appeal  to  their  theoretical  properties  as  approximations  to  the  true  response  function. 
From  this,  it  appears  that  analysts  are  inclined  to  adopt  polynomials  as  a  convenient 
curve-form  panacea  for  all  regression  hypotheses  rather  than  undertake  purposeful 
preanalysis  development  of  hypotheses  specific  to  the  relation  under  study. 

In  the  presumed  presence  of  incentive  for  analysts  to  follow  technical  direction,  these 
findings  suggest  that  this  is  the  current  (1975,  at  least)  state  of  regression  art  insofar  as 
applied  statistics  is  concerned.  One  need  only  consult  related  text  material  to  verify  the 
notion.  Hypothesis  development  is  simply  not  a  topic  covered.  Rather,  integral 
polynomials  are  usually  presented  as  a  package  of  ready-made  hypotheses  that  can  be  ex- 
plored sequentially  and  mechanically  using  least  squares  fitting  procedure.  Discussion 
covering  the  inclusion  of  known  constraints  in  these  hypotheses  is  generally  limited  to  that 
of  the  simple  intercept.  It  is  seldom,  if  ever,  mentioned  that  the  polynomials,  or  log- 
transforms,  etc.  cannot  be  presumed  to  provide  adequate  curve-form  representation  for  all 
hypotheses.  And,  it  is  pertinent  to  note  that  failure  to  discuss  hypothesis  development  and 
its  potentially  complex  outcomes  in  text  material  constitutes  de  facto  disapproval  of  any 
but  simpler  hypotheses,  such  as  those  by  polynomials. 

All  told,  there  is  ample  cause  for  the  hypothesis-related  problems  identified  by  Mead 
and  Pike  (1975)  and  attributed  to  authors  of  biological  research  papers.  It  would  seem, 
however,  that  the  source  of  the  problem  might  more  nearly  be  associated  with  limitations 
in  text  material  and  statistical  curricula  than  with  oversights  by  users  of  currently  prof- 
fered statistical  methodology. 

An  effort  is  made  in  this  paper  to  initiate  solutions  to  the  hypothesis  development  prob- 
lem through  discussion  and  by  showing  analysts  how  to  develop  preanalysis  hypotheses 
graphically,  regardless  of  complexity,  from  accumulated  knowledge  and  existing  data.  For 
biological  subject  matter,  these  hypotheses  will  be  curvilinear  and  highly  interactive  as 
often  as  not. 

Mathematical  formulation  of  graphed  hypotheses  is  accomplished  through  use  of  two 
families  of  curves:  X",  n>  1  (revamped  herein  from  Jensen  and  Homeyer  1970);  and  e  K 
(Jensen  1979),  a  highly  manipulable  function  for  bell-shaped  curves  or  parts  thereof.  De- 
scriptor development  can  be  accomplished  by  analysts  with  a  knowledge  of  basic  algebra 
and  access  to  a  programable  calculator. 
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GRAPHIC  HYPOTHESIS  DEVELOPMENT 


HYPOTHESIS  DEVELOPMENT,  GENERAL 

A  condition  for  valid  statistical  evaluation  of  hypotheses  is  that  they  be  developed  inde- 
pendently of  data  used  for  evaluation.  So,  the  analyst's  first  undertaking  in  a  regression 
analysis  should  be  to  establish  the  strongest  possible  concept  of  the  underlying  form  of 
the  relation  from  accumulated  knowledge  on  the  relations  being  investigated.  It  can  be  an- 
ticipated for  a  broad  array  of  problems  that  these  concepts  will  vary  from  extremely 
weak,  with  little  or  no  basis  for  form  definition,  to  very  firm  expectations  for  even  spe- 
cific algebraic  forms. 

In  the  weak-concept  case,  there  is  little  reason  to  allow  for  any  but  simple  curve-form 
hypotheses.  Low-order  polynomials  in  a  linear  model  will  probably  suffice.  Where  well- 
established  algebraic  forms  already  exist  for  the  form  of  the  relation,  the  exact  hypothesis 
is  already  available  and  evaluation  can  be  initiated.  The  more  common  and  perhaps  more 
challenging  circumstance,  however,  involves  some  intermediate  level  of  prior  knowledge 
where  the  analyst  is  sure  of  some  elements  in  the  relation  being  considered  and  less  sure 
of  others.  His  or  her  deliberations  on  the  nature  of  the  relation  may  result  directly  in  a 
mathematical  hypothesis  or,  more  commonly,  they  may  result  in  a  graphed  form  wherein 
the  analyst  imposes  reasonable  intercept-,  form-,  and  scale-constraints.  These  elements  of 
the  resulting  graph,  then,  must  be  represented  with  suitable  mathematical  accuracy  to 
establish  the  final  intermediate-strength  hypothesis. 

If  the  graphed  hypothesis  involves  unique  curve  forms,  substantial  effort  might  be  re- 
quired to  search  through  existing  functions  like  those  in  the  list  provided  by  Mead  and 
Pike  (1975),  for  those  with  the  capacity  to  mimic  unique  forms  of  the  hypothesis.  (Note 
that  mathematical  descriptors  should  meet  the  acceptance  criteria  of  the  analysts  involved 
[Bartlett  1947;  Draper  and  Hunter  1969]). 

A  generally  more  efficient  alternative  may  be  that  of  analyst-developed  functions  to 
describe  graphed  hypotheses.  Relatively  simple  methods  for  describing  even  complex 
multidimensional  curvilinear  interactions  are  included  in  this  paper.  Curves  of  only  two 
classes  are  utilized:  X",  n>\  \  and  the  very  versatile  three-parameter  bell-shaped  function 
e~K.  Methods  and  related  examples  shown  include  a  variety  of  research  applications  to 
lend  realism  to  the  instruction.  Additional  research  applications  not  cited  in  the  text  are 
given  under  Additional  References  to  provide  the  reader  with  the  opportunity  for  ex- 
panded self-study  of  descriptor  development. 

Given  a  mathematical  hypothesis,  it  can  be  evaluated  using  data  sets  not  used  in  its 
development.  The  evaluation  can  vary  from  a  test  of  general  performance  \ny=b 
(hypotheses),  as  shown  herein,  to  more  complex  formats  designated  to  test  performance 
of  individual  components  of  the  hypothesis.  The  testing  format  adopted  is  likely  to  de- 
pend on  the  complexity  of  the  hypothesis  itself,  the  capacity  of  the  analysts  to  linearize  it, 
and/or  on  the  nature  of  the  fitting  system  to  be  used  (linear,  nonlinear).  Also  evaluation 
can  be  based  on  performance  statistics  compiled  over  repeated  real-world  applications  of 
the  hypothesized  descriptor.  See  Mosteller  and  Tukey  (1977)  and  Box,  Hunter,  and 
Hunter  (1978)  for  discussion  of  empirical  and  structural  (mechanical)  models,  and  fitting 
systems,  etc. 
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GRAPHIC  HYPOTHESIS  DEVELOPMENT  -  TWO  DIMENSIONS 


GRAPHIC  HYPOTHESIS  DEVELOPMENT  -  TWO 
DIMENSIONS 

Accumulated  knowledge  on  the  general  nature  of  the  regression  relation  under  consider- 
ation may  be  used  to  enhance  the  derivation  of  finite  graphed  hypotheses  from  data.  A 
common  analytical  ploy  is  to  dedicate  part  (perhaps  half)  of  a  current  data  set  to  the 
development  of  a  hypothesis,  reserving  the  rest  for  hypothesis  evaluation. 

To  illustrate,  accumulated  knowledge  regarding  the  relation  between  Yand  X  is  such 
that  the  curve  is  expected  to  pass  through  the  origin,  to  have  a  positive  slope  over  X,  and 
to  be  linear-  to  concave-upward  in  form.  Within  these  constraints,  a  specific  estimate  of 
the  true  shape  of  the  relation  is  developed  from  half  of  a  current  data  set  (fig  1). 

Here,  the  analyst  visually  estimates  the  curve  form  expressed  in  the  data  and  simultane- 
ously draws  and  hand-fits  the  form  to  the  data  by  approximate  least  deviations  (Karst 
1958).  All  of  this  is  accomplished  within  the  constraints  of  expectation. 

In  another  case,  prior  knowledge  leads  to  a  strong  expectation  for  a  positive  intercept 
and,  above  this,  a  sigmoid  upward  to  the  right  with  Y asymptoting  at  a  value  of  100  per- 
cent over  large  values  of  X.  Within  these  constraints,  a  specific  estimate  is  again  devel- 
oped from  half  of  a  current  data  set  (fig.  2). 


I  1  1  1  I  '  ' 
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Figure  2.— A  graphed  estimate  of  the  true 
form  of  the  relation  between  X  and  Y. 


The  foregoing  procedures  apply  to  any  curve  shape  in  any  quadrant  or  combinations  of 
quadrants.  They  also  apply  to  two-dimensional  segments  within  three-dimensional  relations. 
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GRAPHIC  HYPOTHESIS  DEVELOPMENT  -  THREE  DIMENSIONS 

-  FOUR  OR  MORE  DIMENSIONS 

GRAPHIC  HYPOTHESIS  DEVELOPMENT  -  THREE 
DIMENSIONS 

Extension  of  two-dimensional  procedures  to  a  three-dimensional  case  is  shown  in  figure  3. 

Note  that  XI  and  XI  are  invariably  correlated  to  a  greater  or  lesser  extent  so  that 
estimates  of  the  X\  effect  will  generally  improve  with  a  reduction  in  breadth  of  the  XI 
segments.  While  only  three  XI  segments  are  shown  in  figure  3  for  the  sake  of  simplicity, 
five  or  more  segments  would  have  been  more  desirable  if  more  data  points  had  been 
available.  The  same  procedure  applies  to  additional  independent  variables,  such  as  X3  in 
figure  4  (following). 


Figure  3.— A  graphed  estimate  of  the  true 
form  of  the  relation  developed  over  X1  within 
each  of  three  segments  of  the  range  in  X2. 


GRAPHIC  HYPOTHESIS  DEVELOPMENT  -  FOUR 
OR  MORE  DPMENSIONS 

Extension  of  three-dimensional  procedures  to  a  four-dimensional  case  is  shown  in  figure 
4.  The  same  procedures  can  be  extended  to  five  or  more  dimensions. 


X3 

I  1  1 
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Figure  4.— A  graphed  estimate  of  the  true  form  of  the  relation  developed 
over  X1  for  three  segments  of  the  range  in  X2  within  each  of  three  seg- 
ments in  the  range  of  X3. 
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MATHEMATICAL  DESCRIPTIONS  OF  GRAPHED  HYPOTHESES  -  TWO  DIMENSIONS 


MATHEMATICAL  DESCRIPTIONS  OF  GRAPHED 
HYPOTHESES  -  TWO  DIMENSIONS 

Regardless  of  the  number  of  dimensions  and  complexity  of  the  relation  to  be  described, 
it  is  ultimately  sorted  into  two-dimensional  effects  that  are  described  individually.  These 
are  finally  reassembled,  additively  or  multiplicatively  as  appropriate,  to  arrive  at  the  whole 
descriptor. 

Each  two-dimensional  element  is  described  with  either  one  or  both  of  two  families  of 
curves,  X",  n>\,  or  e~K,  a  three-parameter  sigmoidal  function  designed  for  easy  use  by 
the  modeler. 


Outline  for  "TWO  DIMENSIONS" 

PAGE 

Class X",  n>\  (Exponential  Curves)   10 

Curves  flat-  to  concave-upward,  positive  slope: 

In  the  right  half  of  two-space  (example  1)   10 

without  constraint  in  two-space  (example  2)   14 
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Curves  flat-  to  concave-upward,  negative  slope: 

In  the  right  half  of  two-space  (example  3)  

without  constraint  in  two-space  (example  4) 
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MATHEMATICAL  DESCRIPTIONS  OF  GRAPHED  HYPOTHESES  -  TWO  DIMENSIONS 
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Two-Component  X",  with  Inflection 
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MATHEMATICAL  DESCRIPTIONS  OF  GRAPHED  HYPOTHESES  -  TWO  DIMENSIONS 


Xn,  n  >  1,  Plus  Sigmoids 
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MATHEMATICAL  DESCRIPTIONS  OF  GRAPHED  HYPOTHESES  -  TWO  DIMENSIONS 


Class- Xn,  n  >  1,  Reference  is  made  to  Matchacurve-2  (Jensen  and  Homeyer  1971)  in  which  use  of  X" 

(Exponential  Curves)     covers  not  only  n>  1,  but  0<n<  1  and  n<0  as  well.  Our  experience  since  that  work  was 

published  has  indicated  that  almost  all  descriptor  needs  involving  X"  are  reasonably  well 
met  with  the  (n  >  1)  option;  so  methods  and  Standards  in  this  summary  publication  in- 
clude the  (n  >  1)  option  only,  a  great  simplification  over  the  Matchacurve-2  presentation 
(Jensen  and  Homeyer,  1971). 

Descriptor  development  for  graphed  hypotheses  accommodates  all  possible  orientations 
in  the  upper-right  quadrant  of  two-dimensional  space  (fig.  5),  and  those  involving  any 
combination  of  these  up  to  four  quadrants. 


-Y- 


Figure  5.— Alternative  orientations  of  X"-curves  in  the  upper-right 
quadrant. 


Regardless  of  the  form  alternative  or  quadrant(s)  involved,  the  X-axis  can  be  altered  as 
needed  to  orient  the  curve  in  the  upper-right  quadrant  and  through  the  origin.  The  analyst 
simply  compares  a  normalized-scale  version  of  this  reoriented  hypothesis  with  graphed 
standard  curves,  or  Standards  (appendix,  marked  for  cut-out),  and  selects  the  two  adja- 
cent ones  most  nearly  like  his  own  in  shape.  Interpolation  between  the  curves  and  their 
associated  powers  (n)  results  in  identification  of  the  best  alternative  offered  by  the  system. 

Note  that  if  duplicates  are  made  of  the  standards,  the  copies  should  be  exactly  the  same 
scale  as  the  Standards. 
Example  1 

Here  (fig.  6),  the  hypothesis  is  a  flat-  ^ 
to  concave-upward,  positive  slope, 
and  occupying  any  vertical  position 
in  the  right  half  of  two-dimensional 
space.  Mathematical  formulation  for 
the  middle  alternative  (zero  • 
y-intercept)  is  detailed  in  figure  7  Y 
and  in  table  1,  but  the  methods  will 
accommodate  any  ^-intercept.  0 

-a 

~0  + 
X 

Figure  6.— Y  =  a  +  bX", 

0<X<XP,n>1. 
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MATHEMATICAL  DESCRIPTIONS  OF  GRAPHED  HYPOTHESES  -  TWO  DIMENSIONS 


In  figure  7,  the  analyst  visually  esti- 
mates the  curve  form  expressed  in 
the  data  and  simultaneously  draws 
and  hand-fits  the  form  to  the  data 
by  approximate  least  deviations.  All 
of  this  is  accomplished  within  the 
constraints  of  expectation.  Mathe- 
matical formulation  procedures 
follow  in  table  1. 


0  200        400        600  800 

X 

Figure  7.—  Y  =  0  +  bX". 

0<X<XP,n>1. 

Table  1.— Mathematical  hypothesis  development 


Overlay  Perform-  Perform- 

Control  X       (Y-a)         control        b(X1)270    a  nee  of    b(X1)264    a  nee  of     a  +  Y2 

points  =  =  points  =  Y1  =  Y2  = 


X  Y  X1  Y1         XP1      Y1/YP1       Y1       (Y1-Y1)       Y2       (Y2-Y1)  Y 

(1)         (2)         (3)         (4)         (5)         (6)         (7)  (8)  (9)  (10)  (11) 


0 

0.0 

0 

0.0 

0.00 

0.00 

0.0 

0.0 

0.0 

0.0 

0.0 

200 

1.3 

200 

1.3 

.25 

.02 

1.6 

+  .3 

1.7 

+  .4 

1.7 

400 

10.6 

400 

10.6 

.50 

.16 

10.2 

-  .4 

10.6 

.0 

10.6 

600 

31.4 

600 

31.4 

.75 

.48 

30.4 

-  1.0 

30.9 

-  .5 

30.9 

800 

66.0 

800 

66.0 

1.00 

1.00 

66.0 

.0 

66.0 

.0 

66.0 

•  Columns  1  and  2:  Select  representative  control  points  (point  coordinates)  from  the 
graphed  hypothesis,  including  XP,  the  largest  value  of  Xto  be  represented  by  the  hypoth- 
esis (800,  in  this  case),  and  YP,  the  corresponding  smoothed  value  of  Y  (66  here).  Five 
control  points  have  been  selected  and  are  listed  in  these  columns,  along  with  computations 
necessary  to  the  search  for  a  suitably  similar  mathematical  curve-form.  Use  fewer  control 
points  where  the  desired  sensitivity  is  less;  use  more  for  greater  sensitivity. 

•  Column  3:  The  A'-scale  remains  unchanged,  but  is  relabeled  as  X\  for  consistency 
with  procedures  for  subsequent,  more  complex  examples.  XP  is  called  XPl  in  the 

A"  1 -scale. 

•  Column  4:  Subtract  the  intercept,  a,  from  all  control  point  7-values  to  produce  Y\. 
Whether  a  is  zero  (as  here)  or  any  positive  or  negative  value,  Y\ -values  will  have  zero  in- 
tercept and  will  be  in  the  upper-right  quadrant  of  two-dimensional  space.  Such  orientation 
of  the  curve  is  necessary  to  use  of  the  Standards.  YP  is  called  YPl  in  the  Yl -scale. 
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MATHEMATICAL  DESCRIPTIONS  OF  GRAPHED  HYPOTHESES  -  TWO  DIMENSIONS 


•  Columns  5  and  6:  Compute  ratios  of  XI  and  Yl  to  XP\  and  YP1,  respectively. 
These  arrays  of  ratios  are  the  normalized  (proportional)  scales  for  the  XI  and  Yl  control 
points.  The  resulting  overlay  points  can  now  be  plotted  on  a  sheet  of  conventional  graph 
paper  (10  x  10  to  Vi  inch  is  best),  at  the  exact  scale  specified  in  figure  8.  Although  not 
necessary,  it  may  help  to  draw  a  smooth  curve  through  the  points.  Call  this  the  overlay 
curve.  It  is  the  normalized  version  of  the  original  hypothesis  plotted  at  the  same  scale 
used  for  the  Standards  (fig.  8). 


7.5"  —  


-i  1  = — i  ■  1  1  1  1  1  1  1  i 

0  1  .2  .3  .4  .5  .6  .7  8  .9  1.0 


X/XP 

Figure  8.— Here  the  graphed  curve  has  been  normalized  and  plotted  at 
the  same  scale  as  the  Standards. 

•  Place  the  overlay  graph  over  the  Standards,  being  certain  that  the  X-  and  Y-axes  are 
matched  exactly.  Find  the  Standards  that  most  nearly  bracket  the  overlay  curve  (fig.  9).  In 
this  case,  n  =  2.5  and  3.0  bracket  the  overlay  curve  nicely.  Note  that  a  light  source  behind 
the  graphs  that  are  being  compared  will  assist  in  the  matching  process.  Where  the  hypoth- 
esized curve  is  not  suitably  matched  by  the  Standards,  a  segmented  or  multicomponent 
^-function,  or  some  portion  of  a  sigmoid  may  suffice  (see  Contents). 

Use  proportional  departure  of  the  overlay  curve  from  the  left  bracketing  Standard  to 
approximate  an  interpolated  value  between  n  =  2.5  and  3.0.  This  is  best  done  at  a  point 
where  curve-form  accuracy  is  most  crucial — the  point  of  sharpest  bend  was  adopted  in 
this  case.  From  figure  9,  the  interpolated  n  is  2.7. 

•  Column  7:  After  the  analyst  has  made  the  first  approximation  of  the  mathematical 
form  of  the  hypothesis,  Yl  =  b(Xl)210,  where6=  yPl/(A-pi)270  =  66/(800)2-70  =  9.576x  10"7, 

A 

he  can  make  associated  estimates,  Yl,  at  each  control  point  X\  value  as  shown. 

•  Column  8:  Differences  between  estimated  and  original  Yl  values  indicate  that 
b(X\)2J0  has  a  preponderance  of  negative  departures  (that  is,  too  much  curvature).  Try 
flatter  curves  (less  power)  iteratively. 
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MATHEMATICAL  DESCRIPTIONS  OF  GRAPHED  HYPOTHESES  -  TWO  DIMENSIONS 


XP 

Figure  9.— The  Standards  best  matched  to  the 
original  curve  bracket  the  overlay  curve. 

•  Column  9:  A  second  estimator,  Y2  =  b(X\)2(A,  where  b  =  YP\(XP\)2M  =  1.430x  10"6, 

A 

produces  Y2  values  as  shown. 

•  Column  10:  Differences  between  estimated  Y2  and  original  Y\  values  are  small  and 
well  balanced  ( +  and  - );  so  b(X\)2M  is  adopted. 

•  Column  11:  The  final  hypothesis  is  then,  Y=a+bX"  =  0+  1.430  x  10-6(AT)264,  or, 
since  X= X\,  Y=  1.430x  10"6(JO2  64- 


70  r 


Figure  10.— The  fitted  model  has  been  plotted 
over  the  original  graphed  form. 

Ntote  that  at  the  scale  shown  in  figure  10,  differences  between  the  graphed  hypothesis 
and  Y=  1.430x  10" 6  {X)2  64  are  barely  discernible.  This  example  produced  a  closer  mathe- 
matical match  for  the  objective  curve  than  one  might  expect  on  the  average,  but  serves  to 
demonstrate  the  procedures  very  well. 
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MATHEMATICAL  DESCRIPTIONS  OF  GRAPHED  HYPOTHESES  -  TWO  DIMENSIONS 


•  As  a  matter  of  fitting-methods  interest,  the  hypothesis  was  then  fitted  by  least 
squares  (forced  through  the  origin)  back  to  the  data  from  which  it  was  partially  derived 

ft 

and:     Y=  c{XT),  where  XT=  b(X\ )" 

c  =  L(XT*  Y)/L(XT)2  =  0.9748 

ft 

or:       Y  =  1 .394  x  10~  6(X)2  M  and  is  plotted  in  figure  10  as  a  dashed  line. 

The  least  squares  model  is  quite  close  to  the  graphed  hypothesis  as  might  be  expected, 
but  the  mathematical  hypothesis  should  not  be  evaluated  statistically  on  this  used  half  of 
the  data  set.  Rather,  the  same  fitting  procedures  can  be  applied  to  the  unused  half  of  the 
data  set  and,  for  this,  conventional  statistical  evaluations  are  appropriate. 

•  Given  the  same  general  hypothesis,  Y=a+  1.430  x  \0~6  (X)2  M,  but  with  a  being  a 
fixed  constraint  other  than  zero,  the  fitting  process  is  identical  to  that  just  specified. 
Assuming  for  example  that  a  is  fixed  at  -  20,  the  final  form  would  be  Y—  -  20  + 
1.394 x  10"6(A')2-64,  0<^T<800. 

•  Where  a  is  variable  (not  fixed),  linear  curve-fitting  techniques  can  be  used  to  estimate 
both  a  and  b  in  Y=a+  b(Xl)n.  Then,  letting  X=  (AT)",  b  =  Zxy/Lx2  and  a=Y-bX, asm 
conventional  notation. 

Example  2 

Assume  a  graphed  hypothesis  like 
that  in  figure  1 1  (flat-  to  concave- 
upward,  positive  slope,  with  negative 
Y- intercept,  a,  and  crossing  into  the 
left  half  of  two-dimensional  space). 
Mathematical  formulation  procedures 
follow  closely  those  of  example  1  and 
are  summarized  in  table  2. 


FigureH.— Y  =  a  +  b(X-kf, 
/c<X<XP,n>1. 

Table  2.— Mathematical  hypothesis  development 


Overlay  Perform-  Perform- 

Control  (X-/c)  (V-a)  control  b(X1)270  a  nee  of  b(X1)264  a  nee  of  a  +  Y2 
points  -  =  points  =  Y1  =  Y2  = 


X  V  X1  Y1       X1/XP1    Y1/YP1       Y1       (Y1-Y1)       Y2       (Y2-Y1)  Y 

(1)         (2)         (3)         (4)         (5)         (6)         (7)  (8)  (9)  (10)  (11) 


-400 

-10.0 

0 

0.0 

0.00 

0.00 

0.0 

0.0 

0.0 

0.0 

-10.0 

-200 

-  8.7 

200 

1.3 

.25 

.02 

1.6 

+  .3 

1.7 

+  .4 

-  8.3 

0 

.6 

400 

10.6 

.50 

.16 

10.2 

-  .4 

10.6 

.0 

.6 

200 

21.4 

600 

31.4 

.75 

.48 

30.4 

-  1.0 

30.9 

-  .5 

20.9 

400 

56.0 

800 

66.0 

1.00 

1.00 

66.0 

.0 

66.0 

.0 

56.0 
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MATHEMATICAL  DESCRIPTIONS  OF  GRAPHED  HYPOTHESES  -  TWO  DIMENSIONS 


Columns  1  and  2:  List  control  points  representative  of  the  curve. 

Shift  the  curve  to  the  upper-right  quadrant  to  match  the  position  of  the  Standards: 

Column  3:  Let  X\  =  (X-k),  where  k  =  the  smallest  value  of  A' to  be  represented  by 

the  equation  ( -  400  here).  Then,  X\  ={X-  [- 400])  =  (X+  400),  and  XPl  =  800. 

Column  4:  Y\  =  ( Y- a)  =  ( Y-  [  -  10])  =  ( Y+  10),  and  YPl  =  66.0. 

Columns  5-10:  Identical  to  example  1 ,  table  2,  columns  5-10. 

Column  11:  The  final  hypothesis  is  then,  Y=  a+  bX'  =  -  10+1 .430*10-  6  (X  +  400)2  M, 
-400<X<400. 
Least  squares  fit. 

With  a  fixed;  Y=c(XT),  where  XT=b(X\)n,  c  =  L(XT*Yl)/(LXT)2,  and 
£  ft 
Y=a  +  c(b(X\)n).  In  this  case,  Y=  -  10  + c{b[X\}"). 

ft 

With  a  variable,  7=0  +  6(^1)".  Letting  ^=(^1)"  and,  under  conventional  notation, 
b  =  Lxy/Zx2  and  a  =  Y-  bX. 


Example  3 

Assume  a  hypothesis  like  that  in 
figure  12  (flat-  to  concave-upward, 
negative  slope,  and  occupying  any 
vertical  position  in  the  right  half  of 
two-dimensional  space).  Mathemati- 
cal formulation  for  the  middle  alter- 
native (g  =  zero)  is  detailed  in  table 
3,  but  the  methods  will  accommo- 
date any  a  value.  Procedures  follow 
closely  those  of  example  1 . 


80  i- 


ERCEPT 


YP  =  66 


800  (X) 
0  (800-X) 


Figure  1 2.— Y  =  a  +  b  (XP  -  X  f. 

0<X<XP,n>1. 


Table  3.— Mathematical  hypothesis  development 


Control 
points 


(XP-X)  (Y-a) 


Overlay 
control 
points 


d(X1)270 


Perform- 
ance of 
Y1 


b(X1)264 


Perform- 
ance of 
Y2 


a  +  Y2 


X 

Y 

X1 

Y1 

X1/XP1 

Y1/YP1 

Y1 

(Y1  -  Y1) 

Y2 

(Y2  -  Y1) 

Y 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

0 

66.0 

800 

66.0 

1.00 

1.00 

66.0 

0.0 

66.0 

0.0 

66.0 

200 

31.4 

600 

31.4 

.75 

.48 

30.4 

-  1.0 

30.9 

-  .5 

30.9 

400 

10.6 

400 

10.6 

.50 

.16 

10.2 

-  .4 

10.6 

.0 

10.6 

600 

1.3 

200 

1.3 

.25 

.02 

1.6 

+  .3 

1.7 

+  .4 

1.7 

800 

.0 

0 

.0 

.00 

.00 

.0 

.0 

.0 

.0 

.0 
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MATHEMATICAL  DESCRIPTIONS  OF  GRAPHED  HYPOTHESES  -  TWO  DIMENSIONS 


Columns  1  and  2:  List  control  points  representative  of  the  curve. 

Shift  the  curve  to  the  upper-right  quadrant  to  match  the  position  of  the  Standards: 

Column  3:  Let  X\  =  (XP -  X)  =  (800 -  X),  and  XP\  =  800. 

Column  4:  Since  a  =  zero  here,  Yl  =  ( Y- a)  =  ( Y—  0)  =  Y,  and  UP\  =  66.0. 

Columns  5-10:  Identical  to  example  1,  table  2,  columns  5-10. 

Column  11:  The  final  hypothesis  is  then, 

Y  =  a  +  b(XP  -  X)n  =  0  +  0. 1430  x  10~ 6(800  -  X)2 ■64,  0  <  X<  800. 
Least  squares  fit.  A 

With  a  fixed  at  any  value;  Y=  c(XT),  where  XT=  b(Xl  )n,c  =  L)XT*  Y\  )/{LXT)2, 
and  Y=a  +  c(b(Xiy). 

With  a  variable;  Y=a+  b(X\)n.  Letting X=  (X\)n  and,  under. conventional  notation, 
b  =  Lxy/Lx2  and  a  =  Y- bX. 


Example  4 

Assume  a  hypothesis  like  that  in  fig- 
ure 13  (flat-  to  concave-upward, 
negative  slope,  with  negative 
F-intercept,  a,  and  crossing  into  the 
left  half  of  two-dimensional  space). 
Mathematical  formulation  proce- 
dures follow  closely,  those  of 
example  1  and  are  summarized  in 
table  4. 


8  0  - 


 56 


vp  =  ee 


=  4® 


-400  0 
800  400 


-10) 


400  (X) 
0  (400-X) 


Figure  13.— Y  =  a  +  (XP-Xf, 
/f<X<XP,n>1. 


Table  4.— Mathematical  hypothesis  development 


Overlay 

d(X1)2-70 

Perform- 

Perform- 

Control 

(XP-X) 

<y-a) 

control 

ance  of 

b(X1)2M 

ance  of 

a  + Y2 

points 

points 

Y1 

Y2 

X 

y 

X1 

Y1 

X1/XP1 

Y1/YP1 

Y1 

(Y1  -  Y1) 

Y2 

(Y2-  Y1) 

Y 

d) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

-400 

56.0 

800 

66.0 

1.00 

1.00 

66.0 

0.0 

66.0 

0.0 

56.0 

-200 

21.4 

600 

31.4 

.75 

.48 

30.4 

-  1.0 

30.9 

-  .5 

20.9 

0 

.6 

400 

10.6 

.50 

.16 

10.2 

-  .4 

10.6 

.0 

.6 

200  - 

8.7 

200 

1.3 

.25 

.02 

1.6 

+  .3 

1.7 

+  .4 

-  8.3 

400  - 

10.0 

0 

.0 

.00 

.00 

.0 

.0 

.0 

.0 

-10.0 
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MATHEMATICAL  DESCRIPTIONS  OF  GRAPHED  HYPOTHESES  -  TWO  DIMENSIONS 


•  Columns  1  and  2:  List  control  points  representative  of  the  curve. 

•  Shift  the  curve  to  the  upper-right  quadrant  to  match  the  position  of  the  Standards: 
Column  3:  Let  XI  =  (XP -X)=  (400-  X),  and  XPl  =  800. 

Column 4:  Let  YI  =(Y-a)  =  (Y-(-  10))  =  (r+  10),  and  >7>1  =  66.0. 

•  Columns  5-10:  Identical  to  example  1,  table  2,  columns  5-10. 

•  Column  11:  The  final  hypothesis  is  then, 

Y=a  +  b(XP-X)n  =  -  10+  1.430X  10'6(400-^)2 ■64,  k<X<400,  and  k  =  the  lower 
applicable  limit  ( -  400  here)  of  the  Arrange. 

•  Least  squares  fit.  ^ 

With  a  fixed  at  any  value;  Y=  c(XT), 

where  XT=  b(Xl)",  c  =  L(XT*Y\)/0:XT)2,  and  Y=a+c(b(XlY). 

With  a  variable,  Y=  a  +  b(Xl)n.  Letting  X  =  (X\)n  and,  under  conventional  notation, 

b  =  Zxy/Zx2  and  a  =  Y-  bX. 


Example  5 

Assume  a  hypothesis  like  that  in  fig- 
ure 14  (flat-  to  convex-upward,  posi- 
tive slope,  and  occupying  any  vertical 
position  in  the  right  half  of  two- 
dimensional  space).  Mathematical 
formulation  for  the  middle  alternative 
is  detailed  in  table  5.  Procedures 
follow  closely  those  of  example  1 . 


-eol  

0  800(X) 
800  0  ( 8 0 0 -  X ) 

Figure  14.—  Y=  YP~b(XP-Xf, 
0<X<XP,n>1. 


Table  5.— Mathematical  hypothesis  development 


Overlay  Perform-  Perform- 

Control  (XP-X)(YP-Y)  control  b(X1)270  ance  of  b(X1)264  a  nee  of  a  +  Y2 
points  =  =  points  =  Y1  =  Y2  = 


X 

Y 

X1 

Y1 

X1/XP1 

Y1/YP1 

Y1 

(Y1  -  Y1) 

Y2 

(Y2-  Y1) 

Y 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

0 

0.0 

800 

66.0 

1.00 

1.00 

66.0 

0.0 

66.0 

0.0 

0.0 

200 

34.6 

600 

31.4 

.75 

.48 

30.4 

-  1.0 

30.9 

-  .5 

35.1 

400 

55.4 

400 

10.6 

.50 

.16 

10.2 

-  .4 

10.6 

.0 

55.4 

600 

64.7 

200 

1.3 

.25 

.02 

1.6 

+  .3 

1.7 

+  .4 

64.3 

800 

66.0 

0 

.0 

.00 

.00 

.0 

.0 

.0 

.0 

66.0 
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MATHEMATICAL  DESCRIPTIONS  OF  GRAPHED  HYPOTHESES  -  TWO  DIMENSIONS 


•  Columns  1  and  2:  List  control  points  representative  of  the  curve. 

•  Shift  the  curve  to  the  upper-right  quadrant  to  match  the  position  of  the  Standards: 
Column  3:  Let  XI  =  (XP -X)  =  (800 - X),  and  XP\  =  800. 

Column  4:  Let  Y\  =  (YP-Y)  =  (66.0 -  Y),  and  YPl  =  66.0. 

•  Columns  5-10:  Identical  to  example  1,  table  2,  columns  5-10. 

•  Column  1 1 :  The  final  hypothesis  is  then, 

Y=  YP-b(XP- X)n  =  66.0-  1.430X  10"6(800- X)2  M,  0<^<800. 

•  Least  squares  fit.  A 

With  YP  fixed  at  any  value;  Y=c(XT),  where  XT=b(X\)",  c  =  (XT*Y\)/(XT)2,  and 

Y=  YP-c(b(X\)"). 

k 

With  YP  variable,  Y=ct+b(Xl)n.  Letting  X=  (XI)"  and,  under  conventional  notation, 
b  =  Lxy/Lx2  and  a=  Y-  bX. 

Example  6 


Assume  a  hypothesis  like  that  in  fig- 
ure 15  (flat-  to  convex-upward,  pos- 
itive slope  with  negative  y-intercept, 
and  crossing  into  the  left  half  of  two- 
dimensional  space).  Mathematical 
formulation  for  this  curve  is  detailed 
in  table  6.  Procedures  closely  follow 
those  of  example  1. 


800     400         O  (400-X) 


Figure  1 5.- Y  =  YP  -  (XP)", 
/c<X<XP,n>1. 


Table  6.— Mathematical  hypothesis  development 


Overlay  Perform-  Perform- 

Control  (XP-X)(YP-Y)  control  b(X1)270  ance  of  b(X1)264  ance  of  a  +  Y2 
points  =  =  points  =  Y1  =  Y2  = 


X  Y  X1  Y1       X1/XP1    Y1/YP1       Y1       (Y1  -  Y1)       Y2       (Y2-Y1)  Y 

(1)  (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10)  (11) 


-400 

-40.0 

800 

66.0 

1.00 

1.00 

66.0 

0.0 

66.0 

0.0 

-40.0 

-200 

-  5.4 

600 

31.4 

.75 

.48 

30.4 

-  1.0 

30.9 

-  .5 

-  4.9 

0 

15.4 

400 

10.6 

.50 

.16 

10.2 

-  .4 

10.6 

.0 

15.4 

200 

24.7 

200 

1.3 

.25 

.02 

1.6 

+  .3 

1.7 

+  .4 

24.3 

400 

26.0 

0 

.0 

.00 

.00 

.0 

.0 

.0 

.0 

26.0 
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MATHEMATICAL  DESCRIPTIONS  OF  GRAPHED  HYPOTHESES  -  TWO  DIMENSIONS 


Columns  1  and  2:  List  control  points  representative  of  the  curve. 

Shift  the  curve  to  the  upper-right  quadrant  to  match  the  position  of  the  Standards: 

Column  3:  Let  X\  =  (XP-X)  =  (400-  X)  and  XP\  =  800. 

Column  4:  Let  Yl  =  (YP-  Y)  =  (66.0-  Y)  and  YP\  =  66.0. 

Columns  5-10:  Identical  to  example  1,  table  2,  columns  5-10. 

Column  1 1 :  The  final  hypothesis  is  then, 

Y=  YP-b(XP-X)n  =  66.0-  1.430X  10"6(400-X)264,  -400<A'<400. 
Least  squares  fit.  A 

With  YP  fixed  at  any  value;  Y=c(XT),  where  XT=b(Xl)\  c  =  L(XT*Yl)/L{XT)2, 
and  Y=  YP-c(b(Xl)n). 

With  YP  variable,  Y=a+  b(Xl)n.  Letting  X=  (X\)n  and,  under  conventional  notation, 
b  =  Lxy/Lx2,  and  a  =  Y-  bX. 


Example  7 

Assume  a  hypothesis  like  that  in 
figure  16  (flat-  to  convex-upward, 
negative  slope  and  occupying  any  ver- 
tical position  in  the  right  half  of  two- 
dimensional  space).  Mathematical 
formulation  for  the  middle  alternative 
is  detailed  in  table  7.  Procedures 
closely  follow  those  of  example  1. 


YP-INTERCEPT 


►  YP-«6 


800 


Figure  16.-Y  =  YP-b(X)", 
0<X<XP,n>1. 


Table  7.— Mathematical  hypothesis  development 


Control 
points 


X  (YP-Y) 


Overlay 
control 
points 


b(X1)270 


Perform- 
ance of 
Y1 


b(X1) 


,2.64 


Perform- 
ance of 
Y2 


a  +  Y2 


X 

y 

X1 

Y1 

X1/XP1 

Y1/YP1 

Y1 

(Y1  -  Y1) 

Y2 

(Y2  -  Y1) 

Y 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

0 

66.0 

0 

0.0 

0.00 

0.00 

0.0 

0.0 

0.0 

0.0 

66.0 

200 

64.7 

200 

1.3 

.25 

.02 

1.6 

+  .3 

1.7 

+  .4 

64.3 

400 

55.4 

400 

10.6 

.50 

.16 

10.2 

-  .4 

10.6 

.0 

55.4 

600 

34.6 

600 

31.4 

.75 

.48 

30.4 

-  1.0 

30.9 

-  .5 

35.1 

800 

.0 

800 

66.0 

1.00 

1.00 

66.0 

.0 

66.0 

.0 

.0 
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Columns  1  and  2:  List  control  points  representative  of  the  curve. 
Shift  the  curve  to  the  upper-right  quadrant  to  match  the  position  of  the  Standards: 
Column  3:  Let  XI  =  X,  and  XP\  =  800. 
Column  4:  Let  Y\  =  (YP-  Y)  =  (66.0  -  Y),  and  YPl  =  66.0. 
Columns  5-10:  Identical  to  example  1,  table  2,  columns  5-10. 
Column  11:  The  final  hypothesis  is  then, 
Y=  YP-b(X)n  =  66.0-  1.430X  \0~6(X)2  (A,  0<^<800. 
•  Least  squares  fit.  A 

With  YP  fixed  at  any  value;  Y=  c(XT),  where  XT=  b{X\)n,  c  =  L(XT*  Y\)/L(XT)2, 

k 

and  r=  YP-c(b(X\)n). 

With  YP  variable,  Y=a+b(XX)n.  Letting  X=  (X\)n  and,  under  conventional  notation, 
b  =  Lxy/Lx2,  and  a  =  Y-  bX. 


Example  8 

Assume  a  hypothesis  like  that  in  fig- 
ure 17  (flat-  to  convex-upward,  nega- 
tive slope  with  negative  F-intercept 
and  crossing  into  the  left  half  of  two- 
dimensional  space).  Mathematical 
formulation  for  this  curve  is  detailed 
in  table  8.  Procedures  closely  follow, 
those  of  example  1 . 


YP-INTERCEPT 


>  YP=66 


-40 


-400         0  400(X) 

O         400  800(X+400) 


Figure  1 7.— Y=  YP-(X~kf, 
lt<X<XP,n>1. 


Table  8.— Mathematical  hypothesis  development 


Overlay 

Perform- 

Perform- 

Control 

(X-K) 

(YP  -  Y) 

control 

b(Xl)2-70 

ance  of 

b(X1)2M 

ance  of 

a  +  Y2 

points 

points 

Y1 

Y2 

X 

y 

X1 

Y1 

X1/XP1 

Y1/YP1 

Y1 

(Y1  -  Y1) 

Y2 

(Y2-  Y1) 

Y 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

0) 

(10) 

(11) 

-400 

26.0 

0 

0.0 

0.00 

0.00 

0.0 

0.0 

0.0 

0.0 

26.0 

-200 

24.7 

200 

1.3 

.25 

.02 

1.6 

+  .3 

1.7 

+  .4 

24.3 

0 

15.4 

400 

10.6 

.50 

.16 

10.2 

-  .4 

10.6 

.0 

15.4 

200  - 

5.4 

600 

31.4 

.75 

.48 

30.4 

-  1.0 

30.9 

-  .5 

-  4.9 

400  - 

40.0 

800 

66.0 

1.00 

1.00 

66.0 

.0 

66.0 

.0 

.0 

•  Columns  1  and  2:  List  control  points  representative  of  the  curve. 

•  Shift  the  curve  to  the  upper-right  quadrant  to  match  the  position  of  the  Standards: 
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Class  £  ~Ky 
(Bell-Shaped  or 
Sigmoidal  Curves) 


Column  3:  Let  X\  -(X-k),  where  k  =  the  smallest  value  of  X  to  be  represented  by 
the  equation,  -  400  here.  Then,  X\  =  (X+  400),  and  XP\  =  800. 
Column  4:  Let  XI  =  ( YP  -  Y)  =  (26.0  -  Y),  and  YPl  =  66.0. 

•  Columns  5-10:  Identical  to  example  1,  table  2,  columns  5-10. 

•  Column  11:  The  final  hypothesis  is  then, 

Y=  YP  -  b(X)n  =  26.0-  1. 430  x  10~6LY)2M,  -400<X<400. 

•  Least  squares  fit.  ^ 

With  YP  fixed  at  any  value;  Y=c(XT),  where  XT=b(Xl)n,  c  =  L(XT*Yl)/L(XT)2, 

ft 

and  Y=  YP  -  db(Xl))n . 

With  YP  variable,  Y=a+b(Xl)n.  Letting  X  =  (AT)"  and,  under  conventional  notation, 
b  =  Lxy/Lx2,  and  a  =  Y-  bX. 

By  reason  of  e  K  symmetry,  shape  parameters  for  this  bell-shaped  function  are  identical 
for  either  sigmoidal  half  (see  e~K  ITS  DERIVATION,  CAPABILITIES,  AND  LIMITA- 
TIONS, appendix  A).  So,  the  e~  A-Standards  are  limited  to  sigmoids  and  these  have  nor- 
malized (proportional)  X-  and  V-scales,  are  in  the  upper  right  quadrant  of  two- 
dimensional  space,  and  have  zero  intercept  and  slope  upward  to  the  right. 

Representative  control  points  (point-coordinates)  are  selected  from  a  graphed  sigmoidal 
hypothesis  and  these  points  are  shifted  to  the  position  of  the  Standards  through  manipula- 
tion of  the  X-  and  F-scales  as  needed.  These  revised  scales  are  normalized  (converted  to 
proportions)  and  the  normalized  control  points  are  plotted  on  graph  paper  at  the  exact 
scale  of  the  Standards.  This  point-representation  of  the  hypothesis  (with  or  without  the 
corresponding  curve  sketched  in)  can  be  laid  over  and  compared  with  alternative  sigmoids 
on  any  page  of  the  Standards.  When  Standards  are  found  that  bracket  the  overlay  curve 
and  are  suitably  similar  to  it,  the  analyst  can  interpolate  between  parameters  of  these 
Standards  to  arrive  at  the  best  sigmoid  provided  by  the  system. 

The  result  is  a  mathematical  hypothesis  that  can  be  fitted  to  new  data  and  evaluated 
statistically.  For  example,  assume  the  graphed  hypothesis  shown  in  figure  18. 
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Recognizing  that  the  shape  of  the  curve  on  each  side  of  the  peak  is  about  the  same,  it 
seems  likely  that  a  symmetrical  curve  fitting  the  left,  most  completely  specified  side  of  the 
relation,  will  also  fit  the  top  segment  of  the  implied  sigmoid  to  the  right  of  the  peak.  So, 
we  only  need  identify  the  mathematical  form  of  the  sigmoid  to  the  left  of  the  peak 
(0<X<  150)  to  describe  the  entire  graphed  hypothesis.  Mathematical  formulation  of  this 
curve  is  detailed  in  table  9. 


Table  9.— Mathematical  hypothesis  development 


Overlay 

Control  X       (Y-a)  control 

points  =  -  points 


n  =  2,/=.67  n  =  2,/=.6701 

XP1=150,  Perform-    XP1=150,  Perform- 

YP1 = 1000  ance  of  YP1 = 1000  ance  of 
a=0,Y1  =          Y1         a=0,Y2=  Y2 


X  Y          X1          Y1      X1/XP1    Y1/YP1   a+YP1*e  K    (Y1  -  Y)    a  +  YP1*e  K  (Y2-Y) 

(1)         (2)         (3)         (4)         (5)         (6)  (7)  (8)  (9)  (10) 


0 

0 

0 

0 

0.00 

0.00 

0.0 

0.0 

0.0 

0.0 

75 

100 

75 

100 

.50 

.10 

100.6 

+  .6 

100.5 

+  .5 

105 

430 

105 

430 

.70 

.43 

437.5 

+  7.5 

437.3 

+  7.3 

120 

700 

120 

700 

.80 

.70 

692.6 

-7.4 

692.4 

-7.6 

150 

1000 

150 

1000 

1.00 

1.00 

1000.0 

.0 

1000.0 

.0 

Procedures  generally  follow  those  of  example  1  for  curves  of  the  A^-class. 

•  Columns  1  and  2:  Five  control  points  representative  of  the  curve  (including  XP,  the 
point  in  X  at  which  Speaks,  and  YP,  the  peak  value  of  Y)  have  been  selected  from  the 
sigmoid  and  are  listed  here.  Use  fewer  control  points  where  the  desired  sensitivity  is  less, 
more  for  greater  sensitivity. 

•  Columns  3  and  4:  The  Standards  for  e~  A are  oriented  at  zero  in  the  upper-right 
quadrant  (as  for  the  ^-Standards)  and,  since  this  curve  already  conforms,  there  are  no 
transformations  of  the  X-  and  F-axes.  So,  X\  =  A' and  71  =  Y,  respectively.  The  X\-  and 
Y\ -variables  will  be  shown  and  used  here  to  be  consistent  with  procedures  for  more  com- 
plex examples  that  follow.  Note  that  XPl  is  the  maximum  value  of  XI  and  that  YP\  is 
the  maximum  value  of  Yl . 

•  Columns  5  and  6:  Compute  ratios  of  X\  and  Yl  to  XPl  and  YPl,  respectively. 
These  arrays  of  ratios  are  the  normalized  (proportional)  scales  for  the  XI-  and  Yl  -control 
points.  The  resulting  overlay  points  can  now  be  plotted  on  a  sheet  of  conventional  graph 
paper  (10  x  10  to  Vi  inch  is  best),  at  the  exact  scale  specified  in  figure  19  (5  inches  for  the 
YP  range  and  Vi  inch  per  1/10  unit  of  Yl,  IVi  inches  for  the  XP  range  and  3/4  inch  per 
1/10  unit  of  XI).  Although  not  necessary,  it  may  help  to  sketch  a  smooth  curve  through 
these  plotted  points  as  shown  in  figure  19. 

•  Use  this  curve  as  an  overlay  for  any  page  of  the  Standards,  being  careful  to  match 
exactly,  the  X-  and  7-axes  of  the  overlay  with  those  of  each  page  examined.  Find  adjoin- 
ing sigmoids  that  are  most  nearly  shaped  like  the  overlay  curve  and  bracket  it  (fig.  20). 

Note  in  the  Standards  that  /  is  the  proportional  distance  of  the  inflection  point  from 
zero  to  1.0  and  n  is  the  power  of  the  negative  exponents  for  e  in  e~K  (appendix).  Curva- 
ture of  the  sigmoids  becomes  more  pronounced  with  increased  n,  and  n  changes  between 
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h  75"  


—  —   -3T-^=T 

0.0  0.1  0  2  0  3  0.4  0  5  0.6  0  7  0  8  0  9  10 

X/XP 

Figure  19.— The  plotted  overlay  points. 


X/XP 


Figure  20.— The  Standards  of  set  (n  =  2)  are  best  matched  to 
the  original  curve.  Several  adjoining  curves  of  this  set 
bracket  the  overlay  curve. 

pages  of  the  Standards.  There  are  nine  sigmoids  on  each  of  the  10  graphs.  In  this  case, 
the  sigmoids  /=  0.6  and  0.7  (fig.  20)  bracket  the  overlay  curve  nicely. 

Having  found  an  acceptable  pair  of  adjoining  sigmoids  that  bracket  the  overlay  curve, 
the  analyst  can  interpolate  between  the  proportional  inflection  points  they  represent  and 
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thus  select  a  sigmoid  more  nearly  matching  the  hypothesis.  Proportional  horizontal  depar- 
ture of  the  overlay  curve  from  the  left  bracketing  sigmoid  at  Y/  YP\  =  0.5  can  be  used  as 
the  basis  for  interpolation.  In  figure  20,  the  overlay  curve  lies  at  about  the  70  percent 
point  between  inflection  points  0.6  and  0.7;  so  the  interpolated  /=  0.67.  (This  same  inter- 
polating principle  can  also  be  used  between  Standards  of  different  sets  for  refined  'V 
estimates  as  the  need  arises.)  Then  substituting  in  the  parent  e  A-equation,  we  have: 


X 
XP 


-  1 


1-/ 


-  e 


\-I 


\-e 


l-I 


X_ 
150 


.33 


.33 


1-e 


.33 


This  is  easily  programed  on  a  small  desk-top  computer  (with  the  number  of  steps  re- 
quired generally  being  less  than  100)  and  solved  for  Fl-values  at  control-point 
X\ -values. 

•  Column  7:  e~K  (ranging  only  from  0  to  1)  is  scaled  to  the  graphed  hypothesis  by 
simply  multiplying  e~  K  by  YPl  ( =  YP  here),  the  height  of  the  sigmoid  at  its  peak.  Add- 
ing the  intercept  a,  produces  Yl. 

•  Column  8:  The  selected  sigmoid  is  reasonably  close  to  that  of  the  hypothesis  with  a 
slight  surplus  of  positive  values  ( +  0.6),  as  shown.  In  the  light  of  the  y-values  involved 
(100,  430,  700),  the  departures  +  0.6,  +  7.5,  and  -  7.4,  respectively,  are  small  and  this 
sigmoid  would  most  likely  be  adopted.  But,  as  might  be  more  appropriate  in  other 
cases,  computer  iteration  can  be  used  to  select  an  even  closer  e~  A-approximation  of  the 
hypothesized  sigmoid. 

•  Columns  9  and  10:  Yl  and  its  performance  are  slightly  better  at  /=  0.6701 .  Assum- 
ing this  to  be  the  final  form  adopted,  we  have: 


-1 


F2  =  0+  1000 


X 
150 
0.3299 


0.001 


1-0.001 


\ 


And,  since  the  right  numerator  is  very  small  and  the  denominator  approaches  1.0, 


72  =  0+  1000*e> 


-  1 


X_ 
150 
0.3299 


,  the  final  hypothesis. 


Note  that  at  the  scale  shown  in  figure  21,  differences  between  the  graphed  hypothesis  and 

A 

Y2  are  not  discernible. 

This  example  produced  a  closer  mathematical  match  for  the  objective  curve  than  one 
might  expect  on  the  average,  but  serves  to  demonstrate  the  procedures  very  well.  As  a 
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1200 


1000 


800 


Y  600- 


400  - 


200- 


ORIGINAL  GRAPHED  FORM 


ORIGINAL 
DATA  POINTS 


POINTS  FROM 
FITTED  MODEL 


20         40         60         80        100        120       140       160  180 


Figure  21.— Points  from  the  fitted  model  have  been 
plotted  over  the  original  graphed  form. 


matter  of  fitting  methods  interest,  the  hypothesis  was  then  fitted  by  least  squares  (forced 
through  the  origin)  back  to  the  data  from  which  it  was  partially  derived. 


And:  Y=cXT,  XT  =  1000* e 
c  =  (L(XT)Y)/UXT)2  = 
almost  no  change. 


—  -i 

150 

0.3299 
0.9960  or  , 


-1 


Y=  0.9960 


\000*e 


X 
150 
0.3299 


which  again,  is  not  distinguishable  from  the  hypothesis  at  the  scale  shown  in  figure  21. 

While  this  result  shows  the  development  process  has  not  deteriorated  the  fit  of  the 
hypothesis  to  the  source  data,  statistical  evaluation  should  be  reserved  for  a  new  data 
set.  The  same  fitting  procedures  are  applicable. 

Alternative  Arrangements       Unique  spacial  arrangement  of  the  bell-shaped  or  sigmoidal  form  may  necessitate  ad- 
in  Space  justment  of  the  F-values  for  an  intercept  and/or  manipulation  of  the  A-scale  to  shift 

the  objective  figure  to  the  upper  right  quadrant  (like  the  Standards).  Commonly  encoun- 
tered spacial  alternatives  and  associated  shifting  techniques  are  shown  in  figures  22-26 
which  follow. 
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•  Here,  the  sigmoid  never  drops  below  the  intercept  a,  regardless  of  whether  a  is 
positive  (fig.  22A)  or  negative  (fig.  22B)  and  is  suitably  adjusted  by  ( Y-a),  as  shown  at 
the  top  of  column  4,  table  9. 
Given  the  same  hypothesis, 


—  -1 
150  

0.3299 

Y  =  0  +  1000*e  ,  but  with  a  being  a  constraint  other  than  zero,  the  fitting 

process  is  identical  to  the  one  just  specified.  If,  for  example,  the  fixed  a=  -20,  the  final 
form  would  be: 


Y=  -20  +  0.9960(1000*? 


X_ 
150 


0.3299 


),  0<^<(1YP=300). 


Where  a  is  variable  (not  fixed),  use  linear  curve-fitting  techniques  to  estimate  both  a  and  b 
in 


Y  =  a  +  b(\000*e 


150 


0.3299 


Then,  letting  X=  1000*e~A,  as  above,  b  =  Lxy/Lx2  and  a=  Y-bX,  as  in  conventional 
notation. 


XP 


2XP 


Figure  23.— Right  half  of  implied  bell- 
shaped  curve  specified. 
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•  Plot  the  mirror  image  of  the  right  side  over  the  range  0<  X<  XP  and  identify  an 
appropriate  sigmoid  using  X-  and  7-values  from  that  range,  as  in  table  9.  The  sigmoid 
adopted  will  apply  equally  on  either  side  of  XP. 

YP  -- 

Y 

0  — 
k 

O 

Figure  24.— 
XP  ±  XP. 

•  Substitute  X\  =  (X-  k )  for  X  in  the  range  XP<X<  (XP+  k ).  Then  follow  table  9 
procedures  to  identify  an  eA -sigmoid  for  the  left  side  of  the  bell-shaped  curve  ( =  right 
side  by  virtue  of  symmetry  here.) 

YP  7'~, 


XP 


Figure  25.— Incomplete  sig- 
moid specified. 


•  Complete  the  implied  sigmoid  in  accord  with  expectation  and  identify  e  A  as  in 
table  9.  The  hypothesis  would  be  presumed  to  apply  in  the  range  0<  X<  XP. 


Figure  26.— Inverted  bell- 
shaped  curve,  or  portion 
thereof. 
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•  Plot  the  negative  departures  of  the  curve  from  the  intercept  as  positive  departures 
on  the  intercept,  that  is,  rotate  the  curve  upward  about  a.  Identify  the  appropriate  sig- 
moid as  shown  previously  and  finally,  Y=a  -  YP*e  K 

Using  Applicable  parts  of        On  occasion,  the  added  effort  required  to  compile  a  multi-component  descriptor  can 
e~  ^-curves  be  avoided  by  utilizing  a  portion  of  a  single-component  curve. 

Assume  the  curve  in  figure  27A  is  to  be  described  mathematically  and  that  no  satis- 
factory matching  curve  has  been  found  using  Xn.  As  an  alternative,  the  analyst  could 
undertake  a  multiple-component  description,  but  another  single-component  possibility 
exists. 

(A)  (B) 


Y  = 

(R)  (R+  500MR) 


(X+10) 

Figure  27.— Using  the  upper  half  of  a  sigmoid:  (A)  the  problem, 
and  (B)  a  solution. 


•  Extend  the  desired  curve  in  figure  27A  such  that  the  resulting  form  approximates 
that  of  a  sigmoid.  Add  a  constant  to  the  abscissa  such  that  XP  ±  XP  covers  the  perti- 
nent range  of  X.  (X+  10)  turns  the  trick  here.  Also,  add  a  suitable  constant,  500  in  this 
case,  to  R  to  shift  the  entire  sigmoid  vertically  to  the  upper-right  quadrant.  For  repre- 
sentative points  from  the  right  half  of  this  curve  (for  example,  figure  27B),  construct  an 
overlay  curve  and  find  a  matching  Standard  for  the  upper  half  of  the  overlay.  The 
adopted  single-component  Standard  fitted  to  the  X-  and  F-values  then  represents  the 
original  curve  in  the  region  0<  X<  10.  The  resulting  descriptor  is 

Y=  1000*e~A  -  500,  0<X<  10, 

or 

(X+  10) 
20 

e  1-/ 
7=1000'  —  — 

1  I. 

•  Procedures  for  any  fraction  of  any  form  are  analogous. 
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Notes  on  e~K  There  may  be  cases  where  the  e~  A-transform  determined  from  the  Standards  is  not 

sufficiently  accurate  for  the  purposes  of  the  analyst.  In  such  cases,  curve-form  descrip- 
tion procedures  beyond  the  scope  of  this  paper  must  be  used. 

The  portion  of  the  curve  to  the  right  of  XP  is  a  mirror  image  of  the  left  in  e  K.  Once 
n  and  the  proportional  inflection  point,  /,  have  been  selected  and  held  constant,  changes 
in  Y  values  depend  only  on  departure  of  X /XP  from  1 .0.  Equal  departures  on  opposite 
sides  of  1 .0  will  result  in  equal  y-values.  For  the  occasions  when  the  analyst's  curve  is 
bell  shaped  and  symmetrical,  we  really  need  only  identify  the  e~  A-transform  for  one 
side  (left),  as  for  figure  18  and  in  table  9  procedures. 

Every  e  A-transform  is  forced  to  zero  at  X  -  0  and  at  X  =  IX P  and  should  only  be 
used  within  these  limits  (appendix,  e  K  DERIVATION). 

When  a  suitable  match  for  a  graphed  two-dimensional  hypothesis  of  the  A""-class  can- 
not be  found  in  the  ^"-Standards,  more  satisfactory  descriptors  can  often  be  developed 
by  matching  unique  portions  of  the  curve  sequentially  and  adding  descriptors  for  the 
parts  to  arrive  at  the  whole.  Also,  development  of  descriptors  for  contiguous  segments 
of  the  X-axis  constitutes  another  flexible  alternative. 

One  example,  a  commonly  encountered  kind  of  curve  form  that  departs  from  the 
single-component  Standards  of  the  A^-class,  is  represented  by  the  graphed  hypothesis  in 
figure  28. 

100  -1 
90 

80  ■ 
70  ■ 
60 
Y     50  ■ 
40 
30 
20 
10 

0  +  1  1  1  ,  1  ,  1  1  1  1 

0  100  200  300  400  500  600  700  800  900  1000 

X 

Figure  28.— A  graphed  hypothesis. 


Assume  this  is  the  curve  to  be  described.  The  discerning  analyst  will  immediately  re- 
ject it  as  a  single-component  X"  form  because  of  the  second  pronounced  bend  in  the 
range  700< X<  1000.  But,  let  us  run  it  through  the  ^''-single  component  procedures  to 
provide  a  contrast  for  the  two-component  form  to  be  developed  subsequently. 


Multiple-Component 
Descriptors 
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Control  points  for  the  curve  are  normalized  in  table  10  and  are  replotted  as  a  5-  by 
7.5-inch  overlay  curve  the  same  size  as  the  Standards  (fig.  29). 

Table  10.— Normalized  control  points  for  figure  28 


Overlay 

Control  X         (Y-a)  control 

points  =  =  points 


X  Y  X1  Y1         X1/XP1  Y1/YP1 

(1)  (2)  (3)  (4)  (5)  (6) 


0 

0 

0 

200 

5.2 

200 

400 

17.0 

400 

600 

33.5 

600 

700 

43.0 

700 

800 

54.8 

800 

900 

71.1 

900 

1000 

100.0 

1000 

0.0 

0.00 

0.00 

5.2 

.20 

.05 

17.0 

.40 

.17 

33.5 

.60 

.34 

43.0 

.70 

.43 

54.8 

.80 

.55 

71.1 

.90 

.71 

100.0 

1.00 

1.00 

X/XP 

Figure  29.— Overlay  curve  for  entire  X-range  compared  to  the 

X"-Standards. 

There  is  no  suitable  match  for  the  overlay  curve  in  the  ^''-Standards;  so  we  resort  to 
description  of  the  graphed  curve  with  two  components.  The  first  component  covers  the 
range  0  <  X<  700  since  that  portion  of  the  curve  in  figure  28  has  only  one  point  of  bend 
and  seems  reasonably  similar  in  conformation  to  the  ^"-Standards.  Next,  X and  Fare 
scaled  to  1. 0  at  .Y=  700  (table  11,  columns  1-4). 
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Table  11.— Mathematical  hypothesis  development 


Overlay  Overlay 
X      (Y-a)      control  points  t       Perform-       control  points         tt  Perform- 

=         =  1st  portion        bX1 65  ance  of  Y1       2nd  portion        cX12  5     ance  of  Y2 


X1         Y1        X1/700      Y1/43.0       Y1       (Y1  -  Y1)    X1/1000  |Col.6|/22.5      Y2      (Y2-  |Col.6|) 

(1)         (2)  (3)  (4)  (5)  (6)  (7)  (8)  (9)  (10) 


0 

0.0 

0.00 

0.00 

0.0 

0.0 

0.0 

0.0 

200 

5.2 

.29 

.12 

5.4 

+  .2 

.0 

-  .2 

400 

17.0 

.57 

.40 

17.1 

+  .1 

.0 

-  .1 

600 

33.5 

.86 

.78 

33.3 

-  .2 

.0 

+  .2 

700 

43.0 

1.00 

1.00 

43.0 

.0 

.3 

+  .3 

800 

54.8 

53.6 

-  1.2 

0.80 

0.05 

1.4 

+  .2 

900 

71.1 

65.1 

-  6.0 

.90 

.27 

6.0 

.0 

1000 

100.0 

77.5 

-22.5 

1.00 

1.00 

22.5 

.0 

t  b5YP1/(XP)1  65  =  43.0/(700)1  65  =  8.69  x  1 0  "  4 

tt  c  =  YP1/(XP)12  5  =  22.5/(10O0)12  5  =  7.1 15  x  10  "  37 

After  plotting  the  first-portion  overlay  curve  and  comparing  it  to  the  ^"-Standards 
(fig.  30),  we  find  that  A1-65  gives  a  good  match  for  the  curve  in  the  0<  X<100  range. 
Scaled  to  43  at  X=  700,  A"1-65  gives  a  fairly  close  approximation  of  the  Y\ -values  (table 
11,  column  6).  As  might  be  expected,  this  relatively  flat  curve  is  too  flat  when  extended 


X/XP 

r-igure  30.— First-portion  overlay  curve  compared  to  the  Xn- 
Standards. 
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to  represent  the  desired  curve  in  the  700<  X<  1000  range  (column  5).  The  differences, 
1 .2,  6.0,  and  22.5,  represent  the  curve  of  values  that  must  still  be  added  to  match  the 
right  side  of  the  desired  curve.  Needed  in  the  descriptor  is  a  second  component  that  is 
essentially  zero  in  the  0  <  X<  700  range  and  about  equal  to  the  curve  of  values  still 
needed  in  the  range  700 <X<  1000. 

Plotting  the  overlay  curve  for  the  second  curve  portion  and  comparing  it  to  the  X" 
-Standards  as  in  figure  31,  we  find  that  Xn  might  match  the  required  curve  fairly  well 


0  .1  .2  .3  .4  .5  .6  .7  .8  .9  1.0 

X/XP 

Figure  31.— Second-portion  overlay  curve  compared  to  the  X"- 
Standards. 

but,  on  trial,  X12-5  performs  a  little  better.  Scaling  Xn  5  to  22.5  at  X=  1000  gives  values 
close  to  those  desired  in  the  700<       1000  range  (table  11,  column  9).  Then  the  com- 
plete descriptor  is: 

y=  (8.691  x  10-4)Arl65  +  (7.115x  lO"37)^125, 

A 

with  final  values  listed  in  table  12.  In  this  case,  Y is  regarded  as  being  satisfactorily 
close  to  Y. 

Table  12.— Model  performance 


.  X 

Y 

y 

0 

0.0 

0.0 

200 

5.2 

5.4 

400 

17.0 

17.1 

600 

33.5 

33.4 

700 

43.0 

43.3 

800 

54.8 

55.0 

900 

71.1 

71.1 

1000 

100.0 

100.0 
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The  computational  ideas  just  shown  are  applicable  to  other  curve  forms,  but  a  larger 
array  of  ideas  may  be  necessary  to  independent  descriptor  development  efforts  by  the 
reader.  The  examples  that  follow  are  presented  with  decreasing  explanatory  detail  as 
seems  appropriate  to  the  stage  of  discussion.  We  will  start  with  an  abbreviated  version 
of  descriptor  procedures  just  presented  in  detail  for  the  hypothesized  curve  in  figure  28. 

Two-Component  A",  No  •  Read  a  set  of  representative  XP  points  from  the  graphed  hypothesis  (fig.  32).  Select 
Inflection  one  A-value  in  this  set,  point  d,  such  that  the  curve  over  the  range  0  to  d  appears  to 

have  about  the  same  general  conformation  (one  point  of  bend)  as  the  X"  Standards. 


Figure  32.— The  graphed  hypothesis. 


•  Compare  the  corresponding  overlay  curve  with  the  A^-Standards  and  select  an  ap- 
propriate n  in  X". 

•  Scale  X"  to /at  d,  b  =  f/d".  Then,  Yl=2  +  bX".  Extend  Y\  over  the  range  d  to  e 
(dashed  line). 

•  Compare  the  overlay  curve  for  the  exponential  curve  of  residuals  over  the  range  0 
to  e  with  the  A"-Standards  and  select  an  appropriate  m  in  X'". 

•  Scale  X7"  to  g  at  e,  so  that  c  =  g/e"1.  Then,  Yl  =  cX'n. 

•  Add  components  Y\  and  Y2.  Then,  Y=  a  +  bX"  +  cX"\  and  a  =  0. 
Alternative  orientations  of  the  above  curve  in  space  can  be  handled  simply,  as 

specified  earlier  in  the  section  entitled  "TWO  DIMENSIONS". 

Two-Component  X",  with  Figure  33  is  a  variant  of  the  preceding  sample.  All  steps  of  the  descriptive  process  still 
Inflection  apply,  except  that  the  second  component  is  subtracted  from  the  first.  Then, 

Y=a+bXn-cXm. 


Again,  alternative  orientations  of  this  curve  in  space  can  be  handled  as  specified  in 
the  AMext,  along  with  appropriate  changes  in  sign  for  the  last  component. 
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Three-Component,  X",  The  descriptor  for  a  curve  that  has  an  extended  flat  midsection  (fig.  34)  can  be 

with  Flat  Central  Segment    treated  as  follows: 

•  Determine  the  intercept,  a,  by  graphic  extension  of  the  flat  segment.  Calculate  the 
slope  of  the  line.  Then, 
Y\=a+bX. 


•  Determine  the  residuals  of  the  curve  from  the  straight  line  for  representative  points 
over  X.  Include  three  or  four  such  points  at  each  end  of  X,  where  the  curve  bends  away 
from  the  straight  line. 

•  For  the  right-end  departures,  plot  an  overlay  curve  and  compare  to  the  X" 
-Standards.  Find  a  suitable  n  in  X"  and  scale  X'  to  h  at  g,  wherein  c  =  h/g".  Then, 

Y2  =  cXn. 

•  For  the  left-end  departures,  plot  an  overlay  curve  and  compare  to  the  X" 
-Standards.  Reverse  the  .Y-axis  by  substituting  K  =  if-X).  Plot  the  left-end  overlay  curve 
and  compare  with  the  A^-Standards.  Find  a  suitable  m  in  Km  and  scale  Km  to  /  at  j; 
such  that,  d=i/jm. 

•  Then,  the  left-end  descriptor  is: 

h  =  dKm. 
Then  the  entire  descriptor  is: 

Y=a  +  bX+cXn  +  dK">,  0<X<g. 
Variants  from  this  curve  include  negative  departures  from  either  or  both  ends  of  the 
straight  line.  In  such  cases,  the  signs  of  the  corresponding  components  simply  become 
negative. 

Alternative  orientations  of  this  figure  in  space  can  be  handled  as  specified  earlier  in 
"TWO  DIMENSIONS". 

Segmentation  with  X"      Where  X"  does  not  provide  a  suitably  accurate  descriptor  for  7  over  the  entire  range 

of  X,  describing  Y=J{X)  for  each  of  several  contiguous  segments  of  the  A'-range,  will 
often  remedy  the  problem.  In  figure  35: 

Y=  1.0  +  6.456X  10  11  (X)5  9,0<X<60,  results  in  departures  from  the  objective 
curve  of  from  +0.03  to  -0.07. 
Y=  1.0,0<^<30,  and 

Y=  1 .0+ 5.708 x  10"4  (X-30)2\  30<X<60,  results  in  departures  of  only  ±0.01. 
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10      20     30      40     50  60 


Figure  35.— The  graphed  hypothesis. 


Straight-Line  Slopes        F°r  any  curve  that  has  an  inflection  and  a  distinctive  peak  (positive  or  negative),  e  A 
(Xn,  n  =  1)  Plus  or       -sigmoids  may  be  appropriate.  (Note  that  sigmoids  may  be  created  with  two  exponen- 
Minus  Sigmoids  t ials  of  opposite  sign,  as  shown  under  multiple  exponentials  with  inflection.)  Sigmoids 

are  shown  here  in  combination  with  linear  effects  for  added  descriptive  power. 

In  this  case  (fig.  36),  we  can  use  a  straight  line  with  slope  b=f/d.  The  intercept  plus 

A 

the  straight  line  is  Y\  =  a  +  bX.  Then  add  to  Y\  a  sigmoid  matching  the  curve  differ- 
ences, Yd,  from  the  linear  function,  scaled  to  g  at  the  point  in  X  (e  here)  where  Yd 

A 

peaks.  (Follow  the  methods  in  table  9.)  The  sigmoid  component  is  Yl-ge~K. 


0  d  X 

X 

Figure  36.— The  graphed  hypothesis. 


The  final  descriptor  is  then: 


Y=a+bX  +  ge  ~K. 


Alternative  arrangements  of  the  sigmoid  about  the  straight  line  would  require  appro- 
priate treatment  (fig.  37a,  b,  and  c). 


x  x 
Figure  37.— The  graphed  hypotheses. 
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For  figure  37a,  XP  =  d  in  the  A'-scale  and  the  sign  of  the  sigmoid  component  would 
become  negative: 

Y=a  +  bX-ge~K. 

For  figure  37b,  P=  (d- X)  would  be  substituted  for  X in  e'K,  XP  =  d in  the  P-scale, 
and  the  sigmoids  would  be  positive: 

Y=a  +  bX+ge~K,0<X<d. 

For  figure  37c,  again,  P=(d-X)  is  substituted  for  X  in  e~K,  XP  =  d  in  the  P-scale,  and 
the  sigmoid  is  subtracted: 

Y=a+bX-ge~K,  0<X<d. 

The  same  sigmoid  alternatives  exist  when  the  slope  of  the  straight  line  is  negative,  the 
only  change  in  the  descriptors  would  be  that  the  sign  of  the  first  component  would  be 
negative. 

X",  n>l,  Plus  Of  course,  exponentials  and  sigmoids  can  be  used  in  combination  (fig.  38).  Where  the 

Sigmoids  left  component  is  an  X"  in  P,  subtracted  from  a,  and  the  right  component  is  a  sigmoid 

in  X  added  to  a,  then, 

Y=a-(f/d")Pn  +  ge-K,  0<X<d. 


o  d  x 

d  0  P  =  (d-X) 

Figure  38.— The  graphed  hypothesis. 


Multiple  Sigmoids 

Unsegmented  Sometimes  the  combination  of  several  sigmoidal  (upright  or  inverted)  components  can 

be  used  effectively  to  describe  asymmetrical  curves.  Each  sigmoid  is  matched  and  scaled 
independently  and  both  are  summed  in  the  final  descriptor. 

For  figure  39a,  P  =  (d—X)  is  substituted  for  X  in  e^K,  for  the  sigmoid  to  the  left  of  c 
and  e^K,  the  sigmoid  to  the  right  of  c,  is  described  in  the  A'-scale.  XP  =  d  in  X,  for  both 
X-  and  P-scales.  Then, 

Y=a+fe[K  +  gexK,  0<X<d. 
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For  figure  39b,  the  descriptor  is  much  the  same  as  for  figure  39a,  but  the  a-value  is 
higher  and  signs  of  the  sigmoidal  effects  are  reversed.  Then, 

Y=  a  -feL  K  -  geR  K,0<X<d. 

Segmented  A  simple  way  of  handling  asymmetry  is  to  break  the  X  range  into  two  parts, 

0<X<XP,  and  XP<X<L  (fig.  40).  Identify  a  suitable  sigmoid  for  each  part  of  the 
skewed,  bell-shaped  curve  and  scale  to  YP  using  the  X-scale  for  the  left  sigmoid  and  the 
(L-X)  scale  for  the  one  on  the  right.  The  result  is  a  two-component  descriptor 

Y=  YPe[K,  0<X<XP 
and 

Y=  YPe^K,  XP<X<L. 


Contiguous  Bell-  Bell-shaped  curves,  in  whole  or  in  part,  can  be  summed  in  a  descriptor.  And,  for 

Shaped  Curves  this  particular  example  (fig.  41),  the  peaks  of  the  two  leftmost  bell-shaped  curves  occur 

at  points  other  than  the  upper  or  lower  extreme  of  the  Arrange,  which  is  not  true  in 
previous  examples  involving  sigmoids.  A  peculiarity  of  e ~K  is  that  its  controlled  use  is 
within  XP  +  XP,  where  XP  is  the  point  in  the  Arrange  at  which  the  K-peak  (either 
positive  or  negative)  occurs,  and  XP±XP  are  the  points  at  which  the  tails  of  the  bell- 
shaped  curve  drop  to  zero.  Beyond  these  points,  the  sigmoid  values  become  negative. 
Then,  where  an  XP  is  located  in  such  a  way  that  XP  ±  XP  does  not  cover  the  pertinent 
range  of  X,  it  becomes  necessary  to  alter  the  A-scale  so  that  it  does.  The  discussion 
below  for  the  foregoing  graph  (fig.  41)  should  clarify  both  the  problem  and  its  solution. 
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A  symmetrical,  bell-shaped  curve  with  values  ranging  from  1 .0  at  XP  to  zero  at 
XP  ±  XP  is  specified  by  e~ K.  With  XP  =  40  for  the  first  bell-shaped  curve,  the  range  of 
application  for  e~K  is  40  ±  40,  or  from  zero  to  80.  Beyond  X- 80,  e~K  would  actually 
dip  below  zero  by  some  unspecified  amount.  The  problem  is  simply  remedied  by  adding 
a  large  enough  constant  to  X  such  that  XP  +  XP  in  the  new  scale  would  cover  the 
original  maximum  of  1 15.  P=  (X+  40)  accomplishes  this  nicely;  so  the  descriptor  would 
be: 

A 

For  the  first  curve,  Y\  =  fe[K,  where  P  replaces  X  in  e~K  and  XP-  80  in  P; 
For  the  second  curve,  725  =  ge2~    where  A' is  used  without  transformation  in  e~K, 
since  XP  =  90  gives  an  applicable  range  for  X  of  zero  to  180,  covering  the  A-extreme 
of  115  as  required; 

For  the  third  curve,  Y3  =  he^K,  where  X  is  also  used  without  transformation  in  e~ K 
since  XP  ±XP  (115  ±  115)  covers  the  pertinent  range  of  A",  0  -  115. 
Then,  the  entire  descriptor  is: 

Y=a+fe^K  +  ge2K+hejK,  0<A<  115 

where:  ef  K  involves  P 
e2"  K  involves  X 
e7K  involves  X. 
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MATHEMATICAL  DESCRIPTION  OF  GRAPHED 
HYPOTHESES  -  THREE  DIMENSIONS 

To  reiterate  what  has  been  said  prevously,  the  analyst's  capability  for  developing  sen- 
sitive multidimensional  hypothesis  models  depends  large  on  his  or  her  2-D  talents.  The 
examples  shown  here  should  provide  enough  procedural  background  to  permit  intuitive 
extension  of  the  system  to  include  more  flexible  curve-form  alternatives.  The  live  data 
examples  included,  serve  to  reinforce  these  ideas  through  application  under  real  data 
circumstances. 


Outline  for  "THREE  DIMENSIONS" 


PAGE 


Sigmoidal  AT-Effect,  with  Constant  Point  of  Peaking  in  XI  and 
Interacting  with  a  Concave-Upward  A'l-Effect 


Symmetrical  Bell-Shaped  Ridge  Over  XI,  with  Moving  Peak  in  X2  and 
Interacting  with  a  Concave-Upward  Al-Effect. 
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Asymmetrical  Bell-Shaped  Ridge  Over  XI,  with  Moving  Peak  in  XI.  This 
is  Described  in  Two  Segments  and  Interacts  with  a  Sigmoidal  A^-Effect. 
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Asymmetrical  Bell-Shaped  Ridge  Over  XI,  with  Moving  Peak  in  XI.  This 
is  Described  in  Three  Segments  and  Interacts  with  a  Sigmoidal  .Y2-Effect. 


Symmetrical,  Inverse  Bell-Shaped  Ridge  Over  XI,  with  Constant  Point  of 
Minimum  in  XI  and  Interacting  with  a  Concave-Upward  A7-Effect. 
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PAGE 


Concave-Upward  Relation  in  Two  Dimensions,  with  Constant  Point  of 
Peaking  in  XI  Over  XI. 


Sigmoidal  Xl-Effect, 
with  Constant  Point 
of  Peaking  in  XI  and 
Interacting  with  a 
Concave-Upward 
^2-Effect 

Contrived  Example  Given  graphed  2-D  hypotheses  over  XI  for  representative  levels  of  X2  (fig.  42),  and 

given  the  basic  2-D  descriptor  for  the  A'l-effect  ( Y=  intercept  +  scalarCA'l  sigmoid))  for 
each  curve  of  the  set,  the  analyst  is  in  a  position  to  formulate  a  descriptor  for  the  im- 
plied surface.  He  simply  expresses  the  changing  intercepts,  scalars,  and  parameter(s)  of 
the  A'l-sigmoid  in  terms  of  X2  and  then  substitutes  appropriately  in  the  basic  2-D 
descriptor.  Pertinent  information  for  the  four  sigmoid  curves  above  has  been  assembled 
in  table  13,  columns  1-6,  by  using  information  read  directly  from  the  graphs  along  with 
the  X"-  and  e~  A-Standards: 


0         2      4      6      8  10 


X  1 

Figure  42.— The  graphed  hypothesis  and  descrip- 
tor performance. 
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Table  13.— Mathematical  hypothesis  development 


Sigmoid  parameters 


X2 

intercept 

XP 

n 

/ 

scalar 

i 

scalar 

(D 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

0 

20 

10 

3.0 

0.90 

25 

0.900 

25.0 

2 

20 

10 

3.0 

.87 

61 

.870 

61.1 

4 

20 

10 

3.0 

.66 

169 

.663 

169.4 

6 

20 

10 

3.0 

.10 

350  ' 

.100 

350.0 

Here,  the  intercept  and  sigmoid  parameters,  XP  and  n,  are  constant  over  X2,  but  the 
scalar  increases  and  /  decreases  exponentially  with  increasing  A7-values.  Using  the  Xn- 
Standards,  /=  0.9  -  0.003704(^2)3  and  scalar  =  25  +  9.028(A'2)2.  These  estimators  per- 
form rather  well  as  may  be  seen  in  columns  7  and  8.  Then  substituting  in  the  basic 

A 

equation  for  Y, 


Y=  intercept  +  scalar(AT  sigmoid), 
Y=  20  +  (25  +  9.028(X2)2)C*T  sigmoid) 


where: 


The  X\  sigmoid : 


CY1/10)-  1 


1-7 


—  e 


\-I 


\-e 


\-I 


and  /,  as  specified  previously,  =0.9-0. 003 704(^2)3. 

Application  limits  for  this  hypothesis  are  0<A'l  <  10,  and  0<  X2<6. 

By  plotting  computer  solutions  for  Y  at  pertinent  combinations  of  X\  and  X2,  we  can 
see  that  the  descriptor  lies  quite  close  to  the  graphed  curves  (fig.  42).  And,  the  entire  pre- 
dicted surface  appears  as  follows  (fig.  43):  .  , 

(370) 


400  f 


8     1 0 


XI 

Figure  43.— The  mathematical  hypothesis. 
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Sigmoidal  .Yl-Effect, 
with  Constant  Point 
of  Peaking  in  XI  and 
Interacting  with  a 
Concave-Upward 
A7-Effect 

Live  Data  Example 


Note  that  the  original  curves  may  be  given  over  X2  instead  of  XI,  in  which  case,  the 
roles  of  XI  and  XI  in  descriptor  development  would  be  interchanged.  Also,  where 
multiple-component  2-D  descriptors  are  involved,  they  are  simply  incorporated,  like  the 
single  ones  here,  in  the  final  descriptor. 

This  descriptor  system  is  applicable  to  virtually  any  surface  and  is  especially  useful  in 
emulating  strong  interactive  relationships.  An  application  to  a  live  data  set  involving 
such  interaction  is  documented  in  the  next  example  as  a  means  of  demonstrating  actual 
use  of  some  of  the  2-D  alternatives  available  to  the  analyst. 


Wildlife  managers  have  expressed  interest  in  the  thermal  cover  provided  large  game 
animals  by  the  live  crowns  of  lodgepole  pine  trees,  as  reported  in  Cole,  D.  M.,  and 
Jensen  1981.  Related  to  such  cover  is  the  height  to  base  of  live  crown  {HBC)  in 
lodgepole  shown  in  figure  44  where  HBC  is  expressed  as  a  function  of  the  average 
height  of  dominant  trees  (//)  in  the  stand  and  stand  density  as  reflected  in  Crown  Com- 
petition Factor  (CCF).  This  hypothesis  reflects  the  expected  effect  for  CCF,  which  is 
positive  and  sigmoidal,  and  that  for  H  is  positive  with  no  strong  expectation  as  to  form. 

These  expected  effects  were  found  to  exist  in  a  data  set  dedicated  to  the  formulation 
of  a  mathematical  hypothesis  for  the  relation.  The  134  observations  involved,  were 
sorted  into  10-foot  height  groups  wherein  HBC  was  plotted  over  CCF(fig.  45). 

The  expected  sigmoidal  CCF-effect  was  clearly  visible  in  the  data  and  appeared  to 
break  downward  from  the  asymptotes  at  about  the  same  point  in  CCF (180-200)  for  all 
height  groups.  The  strongest  expression  of  the  sigmoidal  form  appeared  to  exist  in  the 


100 


0 


40 


80     120    160    200    240    280    320   360  400 


CCF 


Figure  44.— The  mathematical  hypothesis. 
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0  100         200         300         400         500        0  100         200         300         400  500 


CCF  CCF 

Figure  45.— Graphed  hypotheses  and  model  performance  for 
partitioned  data. 


30-  and  50-foot  //-groups.  Following  procedures  specified  in  the  TWO  DIMENSIONS 
text,  control  points  were  read  from  the  approximate  least  deviations-fit,  smoothed 
curves,  and  a  specific  sigmoidal  form  for  each  curve  was  determined  (table  14)  using  the 
e  ^-Standards. 

A  single,  unchanging  sigmoidal  form  was  suitably  accurate  in  describing  CCF  effects 
at  all  heights.  While  the  sigmoid,  described  by  n  =  10,  /=  0.183,  etc.,  provided  the  best 
match  for  the  curve  at  H=  36.5  and  that  sigmoid  with  n  =  10,  /=  0. 190,  etc.,  was  best 
for  the  curve  at  H=  54.4.  An  acceptable  compromise  was  the  sigmoid  with  /  held  cons- 
tant at  0.186  (see  the  last  column  in  table  14).  Then,  on  simplifying, 


HBC 


YP*e 


CCF 
500 


-  1 


1-0.186 


10 


where  YP  is  the  scalar  for  the  sigmoids  over  H.  Plotting  and  smoothing  the  heights  of 
the  graphed  sigmoids  (YP^  at  CCF=  500,  a  slightly  concave-upward  curve  appears  to  ex- 
ist (fig.  46). 

Reading  control  points  from  the  smoothed  curve,  a  matching  ^"-function  was  deter- 
mined (table  15)  from  the  ^"-Standards. 

where,  6  =  48.8/(80.0)' 37,  and  yP  =  0.12056(//)> 37. 
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Table  14.— Mathematical  hypothesis  development 


n  =  10,/  =  0.183,  n  =  10,/  =  0.186, 

Overlay             XP  =  500,  XP  =  500, 

control             control              VP  =18.0,  Performance       YP  =  18.0,  Performance 

points               points            a  =  0,HBC1  =  of  HBC1        a  =  0,HBC2  =        of  HBC2 


CCF   HBC  CCF/500    HBC/18      a  +  YP*eK     (HBCA-HBC)     a+YP'e  K  {HBC2-HBC) 


80 

5.0 

0.16 

0.28 

4.8 

-0.2 

4.6 

-0.4 

120 

10.1 

.24 

.56 

11.1 

+  1.0 

10.9 

+  0.8 

160 

16.1 

.32 

.89 

15.3 

-0.8 

15.3 

-0.8 

500 

18.0 

1.00 

1.00 

18.0 

0.0 

18.0 

0.0 

n  =  10,/  =  0.190,  n  +  10,/  =  0.186, 

Overlay             XP  =  500,  XP  =  500, 

control             control              YP  =  29.5,  Performance       VP  =  29.5,  Performance 

points              points            a  =  0,HBC  =  of  H6C1       a  =  0.HBC2  =        of  HBC2 


CCF    HBC  CCF/500  HBC/29.5     a+YP*e~K     (HBC1  -  HBC)     a  +  YP*eK  (HBC2-HBC) 


54.4 


80  7.7 

120  16.3 

160  25.4 

500  29.5 


0.16 
.24 
.32 

1.00 


0.26 
.55 
.86 

1.00 


7.0 
17.4 
24.8 
29.5 


-0.7 
+  1.1 
-0.6 
0.0 


7.5 
17.8 
25.0 
29.5 


-0.2 
+  1.5 
-0.4 
0.0 


0        10     20     30     40    50     60     70  80 

H 

Figure  46.— Graphed  hypothesis  and  model 
performance  for  YP. 
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Table  15.— Development  of  VP 


Overlay 

Control  control  Performance 

points  points  b(H)1-37        of  YP 


YP        H/80.0       VP/48.8         YP  (YP-YP) 


25.1 

9.3 

0.31 

0.19 

10.0 

'  +0.7 

36.5 

16.4 

.46 

.34 

16.7 

+  .3 

44.6 

21.0 

.56 

.43 

21.9 

+  .9 

54.4 

29.4 

.68 

.60 

28.8 

-  .6 

70.2 

41.5 

.88 

.85 

40.8 

-  .7 

74.8 

44.5 

.94 

.91 

44.5 

.0 

80.0 

48.8 

1.00 

1.00 

48.8 

.0 

The  complete  hypothesis  is  then: 


HBC  =(0.12056(Hy  37)*e 


CCF 
500 


0.814 


with  limits:  0<CCF<500;  0<H<  100. 

The  performance  of  this  model  on  the  data  from  which  it  was  partially  derived 
(/?2  =  0.87)  is  shown  as  a  dashed  line  within  each  //-group  in  figure  45.  The  final  hypo- 
thesis presented  in  figure  44  is  produced  from  model  values  in  table  16. 


Table  16.— HBC  values  from  the  model 


CCF 


AVH 


40        60        80       100       120       140       160       180       200       300  500 


20  1235667788 
30  2  3  6  8        10        11         12        13        13  13 


40 

1 

2 

5 

9 

12 

15 

17 

18 

19 

20 

20 

50 

1 

3 

7 

12 

16 

20 

23 

25 

26 

27 

27 

60 

1 

4 

9 

15 

21 

26 

29 

32 

33 

35 

35 

70 

1 

5 

11 

18 

26 

32 

36 

39 

41 

43 

43 

80 

2 

6 

13 

22 

31 

38 

43 

47 

49 

51 

51 

90 

2 

7 

15 

26 

36 

45 

51 

55 

57 

60 

60 

100 

2 

8 

18 

30 

42 

52 

59 

64 

66 

69 

70 

- 1 


t  Limits  within  which  data  were  included. 
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The  mathematical  hypothesis,  then,  can  be  adjusted  to  data  sets  from  different  popu- 
lations of  trees  through  simple  rescaling  with  b  =  LXY/I,X2,  as  is  indicated  in  the  TWO 
DIMENSIONS  text. 

Symmetrical 
Bell-Shaped  Ridge 
Over  X\  with  Moving 
Peak  in  XI  and 
Interacting  with  a 
Concave-Upward 
^2-Effect 

Live  Data  Example  George  and  Blakely  (1973)  reported  on  the  calibration  of  fire-retardant  dispersal  from 

delivery  aircraft.  On  each  of  16  passes  over  a  level  target  area,  600  gal  of  fire  retardant 
were  dropped  from  an  aircraft.  Square-foot  coverage  and  volume  of  material  reaching 
the  ground  per  100  ft2  were  measured  for  each  aerial  drop  along  with  drop  height, 
windspeed,  wind  direction,  aircraft  speed,  temperature,  and  humidity.  Interest  finally 
centered  on  the  change  in  coverage  of  >2  gal/100  ft2  over  a  controlled  range  of  drop 
heights  within  a  partially  controlled  range  of  windspeeds.  (The  effects  of  all  other  vari- 
ables were  either  negligible  or  unidentifiable  in  this  data  set). 

If  the  plane  were  rolling  along  the  ground  (drop  height  =  3  ft)  at  normal  flying  speed 
and  in  the  absence  of  wind,  it  would  be  expected  that  a  600-gal  drop  would  be  distrib- 
uted over  a  relatively  narrow  strip  of  ground  and  that  the  area  covered  by  >2  gal/100 
ft2  would  be  held  to  some  nominal  value.  Increased  aircraft  (and  drop)  height  should 
result  in  greater  dispersal;  maximum  coverage  should  be  reached  at  some  optimal 
height.  As  drop  height  is  increased  beyond  the  optimum,  dispersal  of  the  retardant  from 
air  friction  and  evaporation  should  become  more  complete;  coverage  should  finally 
reach  zero  at  some  relatively  great  drop  height. 

So  expectation  at  zero  wind  is  for  a  bell-shaped  curve  truncated  to  the  left  of  the 
peak.  With  increasing  wind,  optimal  drop  height  and  peak  coverage  should  diminish, 
since  dispersal  of  the  retardant  would  be  accentuated  by  the  horizontal  shearing  force  of 
the  wind.  Thus,  the  expected  surface  would  be  as  shown  in  figure  47. 


6 

DROP  HEIGHT 
(FEET) 


Figure  47.— The  general  expected  relation. 
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Using  these  expectations  as  guides,  coverage,  C,  was  curved  over  height,  H,  fitting  the 
curves  to  actual  data  points  by  approximate  least  deviations  (Karst  1958)  for  each  of 
three  data  groups  in  wind,  W  (fig.  48). 
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3.2  1.7-4.8 
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H 

Figure  48.— Expected  trends  fitted  through  the  data. 


These  curves,  plotted  at  their  respective  W-means,  were  connected  over  W  to  form  the 
implied  surface,  figure  49.  And,  it  can  be  seen  that  major  features  of  the  expectation 
are  reflected  in  the  data — even  within  the  limited  range  of  H  and  W. 


100  200  300 

-H- 

Figure  49.— The  graphed  hypothesis. 


MATHEMATICAL  DESCRIPTION  OF  GRAPHED  HYPOTHESES  -  THREE  DIMENSIONS 


Figure  49,  then,  is  the  basic  hypothesis  for  which  a  descriptor  was  developed.  Assum- 
ing that  the  curves  over  H  (fig.  48)  could  be  suitably  represented  by  segments  of  sym- 
metrical bell-shaped  forms,  e~K  was  deemed  applicable  and  the  identification  of  match- 
ing ^-transforms  was  begun. 

The  e-transform  is  limited  in  that  it  will  yield  values  of  from  zero  to  one  in  the  full 
bell-shaped  form  only  within  the  range  HP±HP,  where  HP  is  the  point  in  H  at  which 
the  C-values  peak.  And,  it  was  a  unique  feature  of  this  example  that  the  maximum 
//-values  for  all  three  curves  of  the  set  exceed  2HP  (fig.  50). 


toooo 


POINTS  AT  WHICH 
PEAK  COVERAGES 
OCCUR  XP, 


Figure  50.— Graphed  trends  from  figures  49  and  50,  extended. 


The  broadest  curve  (for  W=  0.5)  was  estimated  to  peak  at  //=  273.  In  this  case,  a 
matching  e-transform  would  have  dropped  below  zero  where  //exceeded  1HP,  or 
2(273)  =  546,  long  before  the  approximate  graphed  limit  of  H=  1 150  was  reached.  One 
solution  here  is  to  reverse  the  //-scale  using  as  a  base,  the  largest  //-value  of  the 
broadest  sigmoid  likely  to  be  encountered  in  the  set  of  interest.  Since  an  even  broader 
sigmoid  than  W=  0.5  is  expected  to  exist  at  W=Q,  the  largest  //-value  is  estimated  at 
1150+  100=  1250.  The  //-axis  is  then  reversed  using  XI  =  (1250-//)  where,  as  may  be 
seen  in  figure  50,  XP\  =  977,  1072,  and  1 170  for  W=  0.5,  W=  3.2,  and  W=  9.6,  respec- 
tively. Now,  as  required  for  each  sigmoid,  XP\  ±  XP\  includes  both  the  largest  and 
smallest  value  of  X\  that  will  be  encountered,  and  X\  =  (1250-  X)  is  used  as  the  basis 
for  the  bell-shaped  //-effect  descriptor  (e~K).  The  coverage  estimator  is  then  specified 
as, 


where, 


C=CP1 


(X\/XP\)-  1 


1-/ 


—  e 


-(1/(1-/))" 


-(1/(1-/))" 


c- 

XI-- 
XPl- 
CP\- 

I- 


n  = 


coverage  in  square  feet  (>2  gal/ 100  ft) 
(1250-//) 

the  point  in  XI  at  which  C  peaks 

coverage  peak  at  XP\,  scalar  for  the  sigmoid 

the  proportional  point  in  X\  at  which  the  inflection  point  of  the  bell-shaped 
curve  occurs  (a)  in  the  range  X\  =0  to  XP\  for  the  left  half  of  the  curve  or  (b) 
in  the  range  XPl  to  2XP\  for  the  right  half  (fig.  50). 

the  power  of  the  transform  that  dictates  the  degree  of  curvature  above  and  below 
any  inflection  point. 
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Since  the  API-  and  CP1 -values  read  from  the  three  curves  in  figure  50  varied  with  wind 
in  accord  with  expectation,  they  were  each  expressed  as  a  suitable  function  of  wind. 
Multiple-component  A^-descriptors  were  utilized. 

XPl-  and  CP1-  values  from  the  foregoing  equations  along  with  representative  paired  Al- 
and C-values  read  from  the  left  side  of  each  bell-shaped  curve  were  the  basis  for  overlay 
curves.  Standards  for  e~K  that  were  suitably  similar  to  the  overlay  curves  were  identified., 
and  the  corresponding  n  (held  constant)  and  /-values  were  recorded.  Also,  varying  over 
wind  in  accord  with  expectation,  /  was  described  as  a  function  of  wind  to  complete  the  in- 
puts to  the  coverage  estimator.  The  details  of  these  computations  follow. 

For  X PI  =f{W) 

Paired  W-  and  A'Pl-values  (0.5,  977;  3.2,  1072;  and  9.6,  1 170)  were  read  from  figure 
50.  These  points  were  plotted  and  a  smooth  curve  was  drawn  through  them  (fig.  51) 
that  extended  over  the  relevant  range  of  W(0  —  15  m\/h  ). 
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Figure  51.— Describing  the  XP-hypothesis. 
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The  initial  descriptor  approach  here  was  to  subtract  a  single-component  exponential 
curve  on  the  reversed  W'-scale  from  the  intercept,  1 190  (table  17,  columns  1-7).  W  was 
reversed  in  XI  =  (15  -  W)  to  aline  large  values  of  the  independent  variable  with  those  of 
API.  Next,  the  three  data  point  A'l-values  and  associated  coverage  differences  (ab- 
solute) from  1 190  were  scaled  to  1 .0  at  XI  =  14.5.  The  overlay  curve  was  not  well 
matched  by  any  curve  in  the  A""-Standards;  so  a  two-component  approach  was  adopted 
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Table  17.— Descriptor  development  for  XP1 


Overlay 

t 

w 

XP1 

(15 -IV) 

(1190 -XP1) 

control  points 

b(X1)2-25 

= 

isi  segment 

v 
A 

y 

A  1 

V-i 

Vi  /  V  D1                V-1  IVD1 

V1 
T  1 

(1) 

(2) 

(3) 

(4) 

(5)  (6) 

(7) 

0.0 

15.0 

.5 

977 

14.5 

213 
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1072 

11.8 
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1.00  1.00 
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9.6 

1170 

5.4 

20 

.46  .17 
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(15 -IV)       (Y1-Y)  control  points  c(X1)20  XP1 

=  =  2nd  segment  =  = 


X1  X2  X1/XP2         Y2IYP2  Y2 


(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

15.0 

936 

14.5 

25 

1.00 

1.00 

25 

977 

11.8 

0 

.81 

.00 

0 

1072 

5.4 

0 

.37 

.00 

0 

1170 

.0 

0 

.00 

.00 

0 

1190 

t  b=  118/(11. 8)2  25  =  0.4572 
tt  c  =  25/(1 4.5)20=  1.4811  x  10  "22 
tttXPI  =  1 190 -£>X12 25 -cX120— which  gives  a  perfect  match  for  the  original  XPI-values  at  control  points  0-14.5  in 
wind. 


(table  17,  columns  8-13).  This  time,  the  differences  were  fitted  at  XI  =  0,  5.4,  and  11.8 
to  start  with  (see  the  dotted  line  below  the  intercept,  1 190,  in  figure  51).  Then,  the  dif- 
ference (213  -  188  =  25)  between  the  first -component  curve  and  the  desired  curve  was 
added  in  a  second  component.  Computations  are  summarized  in  table  17,  columns  8-13, 
and  overlay  curves  are  shown  in  figure  52. 
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Y/YP  5 


r     OVERLAY  CURVE 

FOR  SECOND 

COMPONENT 

OVERLAY 

-""  POINT  {.81.  0) 

Figure  52.— First-  and  second-component  overlay  curves  com- 
pared to  the  /"-Standards. 

The  completed  function  was  then 

XPl  =  1190-0.4572(15-  W)2  25 -  1.4811  x  10"22(15-  W)20,  0<  W<  15. 


For  CP!  =f(W) 

Paired  W-  and  CPl-values  (0.5,  8900;  3.2,  8100;  and  9.6,  6970)  were  read  from  figure 
50.  These  points  were  plotted  and  a  smooth  curve  was  drawn  through  them  (fig.  53) 
that  extended  qver  the  relevant  range  of  wind,  0-15  mi/h. 


C  P  i) 


9000 


si 


8000  - 


7000- 

6390 
6000 


©  8900 


0.5 
1 

14.5 


3.2 
i 

11.8 


9.6 
1 

5.4 


Figure  53.— A  flat  curve,  X1 1,  matches  the  lower  end 
(9.6  <  kV<15.0)  of  the  graphed  hypothesis. 


15.0  -  (W) 
0.0   -  (15-W)  =  X1 
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Again,  a  single-component  exponential  was  first  tried  for  the  descriptor  here  (table 
18,  columns  1-7).  Reversing  the  W^axis  to  XI  =  (15  -  W)  to  aline  large  values  of  the  in- 
dependent variable  with  those  of  CPl,  the  XI-  and  CPl -values  were  scaled  to  1.0  at 
XI  =  14.5.  An  overlay  curve  was  constructed  and  compared  to  the  ^"-Standards.  But, 
none  of  the  Standards  were  acceptably  close  to  the  overlay,  so  a  two-component  model 
was  tried  next  (table  18). 

To  start  with,  the  lower  end  of  the  curve  (0.0<  X\  <  15.4)  was  fitted  with  an  intui- 
tively selected  flat  form  (XV'1);  it  was  obvious  without  scaling  and  overlays  that  the 
curve  was  extremely  flat  in  that  range  of  XI  (fig.  53).  The  component  X\ 11  was  then 
scaled  to  the  difference  (580)  between  the  intercept  (6390)  and  the  desired  curve  height 
(6970)  at  XI  =  5.4.  So  this  first  component  (table  18,  column  7),  along  with  the  inter- 
cept, 6390,  constitutes  the  partial  descriptor: 

CPl  =  6390+ b(Xl)lA. 


Table  18.— Descriptor  development  for  CP1 


Overlay  f 
(15 -IV)    (CP1  -6390)        control  points  bX1 1 

=  =  1st  segment  = 


W  CP1  X1  Y1  X1/5.4          Y1/580  Y1 

(D  (2)  (3)  (4)  (5)  (6)  (7) 


0.0 

.5 

8900 

14.5 

2510 

1719 

3.2 

8100 

11.8 

1710 

1370 

9.6 

6970 

5.4 

580 

1.00 

1.00 

580 

15.0 

6390 

.0 

0 

.00 

.00 

0 

Overlay  tt  tt/f 

(15  -IV)       (Y-Y1)  control  points  cX4  CP1 

-  =  2nd  segment  = 


X1  Y1  X1/14.5         Y27791  Y2 


(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

15.0 

9080 

14.5 

791 

1.00 

1.00 

791 

8900 

11.8 

340 

.81 

.43 

347 

8107 

5.4 

0 

.37 

.00 

15 

6985 

.0 

0 

.00 

.00 

0 

6390 

t  b  =  580/(5.4)1  1  =90.739 
tt  c  =  791/(1 4.5)4  =  0.01 7894 
ttt  CP2  =  6390  +  i>X11  1  +cX14 
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This  is  the  dashed  line  above  the  intercept  (6390)  in  figure  53.  The  remaining  differ- 
ences, 340  and  791 ,  and  the  associated  A"l-values  were  listed  and  scaled  to  1 .0  at 
X\  =  14.5  (table  18,  columns  8-11).  The  overlay  curve  was  then  plotted  and  compared  to 
the  ^-Standards  (fig.  54). 


The  curve  n  =  4.0  matched  the  overlay  fairly  well.  So  (AT)4  was  scaled  to  the  difference, 
791,  at  XI  =  14.5  (table  18,  column  12).  This  was  added  to  the  intercept  and  the  first 

A 

component  to  arrive  at  the  two-component  descriptor,  CP1  =  6390  +  90.739(15  -  W)]  A 
+  0.017894(15-  W)4,  <  W<  15.  Note  that  the  related  estimates  in  table  18,  column  13  were 
regarded  as  being  satisfactorily  close  to  the  original  objective  values  in  column  2. 

Matching  the  Bell-Shaped  Curves 

Tables  19  through  21  list  XP\-  and  CPl-values  from  the  equations  developed  in  the 
foregoing  sections  along  with  some  representative  H-  and  C-values  from  the  data-based 
portions  of  the  bell-shaped  curves  in  figure  50.  Also  listed  are  C-value  estimates  ob- 
tained by  using  matching  e~  A-Standards  (see  the  overlay  curves  in  figure  55). 

As  may  be  seen  in  tables  19-21,  the  coveragevalues  of  the  curve  for  each  wind  group 
were  scaled  to  1 .0  by  using  the  corresponding  CP\  as  the  maximum.  The  values  of  X\ 
or  of  (2AP1  -  AT)  for  each  wind  group  were  similarly  scaled  by  using  the  corresponding 

A 

XP\  as  the  maximum.  After  the  paired  proportions  for  //-levels  within  each  wind  group 
were  plotted  and  smoothed,  the  resulting  overlay  curves  were  compared  to  the  e~K 
-Standards  for  n  =  2  (fig.  55).  Although  better  matches  for  some  of  these  overlays  can 
be  found  in  Standards  that  have  n  other  than  2.0,  it  was  decided  to  hold  n  constant 
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Table  19.— Descriptor  development  for  C,  W-0.5 


t 

++ 

1  1 

Overlay 

ttt 

I  I  I 

(1250  -H)  (2XP1  -X1) 

(C-a) 

control  points 

w 

H 

C 

X1 

T 

Y1 

77977 

V1/8900 

Y1 

0.5 

50 

6860 

1200 

754 

6860 

0.77 

0.77 

7137 

100 

7800 

1150 

804 

7800 

.82 

.88 

7793 

150 

8450 

1100 

854 

8450 

.87 

.95 

8322 

200 

8770 

1050 

904 

8770 

.93 

.99 

8692 

273 

I 89Q0 | 

f977l 

(9771 

|  8900 1 

1.00 

1.00 

8900 

t  Temporary  transform  for  plotting  the  overlay  curve  only.  The  temporary  Xlvalues  (7)  must  be  in  ascending  order 
toward  XP1  =977  to  match  the  position  of  the  e"  K-Standards  in  two-space, 
ft  All  sigmoids  of  the  set  are  oriented  at  a  =  zero,  so  (C-a)  =  C  =  V1. 
ttt  Estimated  coverage  using  e~K  with  X1  =(1250-  H),  /  =  0.51,  n  =  2.Q,  XP1  =977,  CP1  =8900,  and  a  zero  intercept, 
tttt  Smoothed  values  from  XP1  and  CP1,  each  a  function  of  wind  (shown  earlier  in  the  text). 


Table  20.— Descriptor  development  for  C,  W  =  3.2 


Overlay  t 
(1250 -  H)  (2XP1  - X1)    (C-a)         control  points  Y1 


w 

H 

C 

X1 

T 

Y1 

T71072  Y1/8107 

3.2 

50 

6730 

1200 

944 

6730 

0.88 

0.83 

6843 

100 

7600 

1150 

944 

7600 

.93 

.94 

7612 

150 

8040 

1100 

1044 

8040 

.97 

.99 

8042 

178 

181071 

1 1072  | 

1 1072  | 

| 8107  | 

1.00 

1.00 

8107 

t  Estimated  coverage  using  e 

"K  with  X1 

=  (1250- 

H),  /  =  0.71,  n 

=  2.0,  XP1  = 

1072,  CP1  = 

8107,  and  a  zero  intercept. 

over  wind.  Then,  the  /-values  corresponding  to  the  best  matching  curve  alternatives  in 
this  set,  were  estimated  to  be  0.51,  0.71,  and  0.88  for  ^=0.5,  3.2,  and  9.6,  respective- 
ly. On  scaling  and  checking  these  and  alternative  curves  (by  means  of  a  small  computer) 
at  the  tabled  //-values  within  H^-group,  some  improvement  was  achieved  with  the  final 
/-array,  0.51,  0.71,  and  0.88.  The  last  column  in  each  of  tables  19-21  shows  the  final 
scaled  coverage  estimates.  There  was  close  proximity  of  these  values  to  actual  coverage, 
C,  so  /= /( W)  was  developed  next. 
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Table  21.— Descriptor  development  for  C,  W  =  9.6 


-  THREE  DIMENSIONS 


t  Overlay  ft 

(1250  -  H)      (C-a)  control  points  Y 


H 

C 

X1 

Y1 

X1/1170 

Y1/6985 

80 

| 6985 | 

|1170| 

| 6985 | 

1.00 

1.00 

6985 

100 

6820 

1150 

6820 

.98 

.98 

6845 

150 

5300 

1100 

5300 

.94 

.76 

5448 

200 

3270 

1050 

3270 

.90 

.47 

3364 

250 

1830 

1000 

1830 

.85 

.26 

1612 

300 

997 

950 

997 

.81 

.14 

600 

t  Since  the  right  half  of  this  curve  (fig.  50)  was  read,  X1 -values  are  already  ascending  toward  XP1  =1170  and  no 
temporary  transform  of  X1  is  necessary  for  plotting  the  overlay  curve. 

tt  Estimated  coverage  using  e~K  with  X1  =(1250-H),  /  =  0.88,  n  =  2.0,  XP1  =  1170,  CP1  =6985,  and  a  zero  intercept. 


0.0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 


X/XP 

Figure  55.— Standards  for  n  =  2.0  and  proportional  departure 
of  inflection  points  in  X,  (X/),  from  XP,  (XIIXP ). 
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4  <> 

I                              I  i 
0.5                       3.2  9.6 
 '  ,  1  r-1  1  .  .  1  1  1  1 — .  1  1  .  ,  1 

0  1  2         3         4         5         6  7  8         9         10        11        12        13        14  15 

•W- 

Figure  56.— Initial  and  final  curves  of  /  over  W. 


For  I =f(W) 

Using  the  paired  wind  and  /-values  just  determined  (0.5,  0.51;  3.2,  0.71;  and  9.6, 
0.88),  /  was  plotted  over  W^and  a  smooth  initial  curve  drawn  through  the  three  data 
points  (fig.  56). 

On  checking  the  performance  of  e~K  over  the  range  of  W,  with  smoothed  CP-,  XP-, 
and  initial  /-values,  it  was  found  that  the  left  edge  of  the  surface  undulated  unaccept- 
ably  (at  //=  50,  0<  W<  5).  /-values  were  adjusted  iteratively  to  minimize  the  problem  and 
the  new  (final)  set  of  points  shown  in  figure  56  were  established.  The  smooth  (final) 
curve  through  these  points  became  the  objective  curve  for  which  a  descriptor  was 
developed. 

This  curve  was  judged  to  require  a  multicomponent  model  because  neither  the  sig- 
moids nor  the  exponentials  have  the  flattened  central  segment  found  there.  The  flat 
segment  was  represented  by  a  straight  line;  the  negative  differences  at  the  left  end  by  an 
exponential;  and  those  at  the  right  by  the  lower  portion  of  a  sigmoid  (fig.  57). 

The  Flat  Segment. — A  straight  line  was  drawn  along  the  final  curve.  The  intercept, 
a  =  0.62,  and  values  necessary  to  the  coefficient  for  W,  (1.054-  0.620)/ 15  =  0.0289, 

A 

were  read  directly  from  the  graph.  Then  71,  the  first  descriptor  component, 
=  0.620  +  0.0289(  W),  (see  table  22,  column  3). 

The  Left  End. — This  residual  portion  of  the  curve  ( Y2)  was  satisfactorily  matched  by 
a  single  exponential  of  the  reversed  H^axis  (Y2),  scaled  at  X\  =  14.5  (table  22,  columns 
4-8). 

So  the  partial  form  at  this  point  is: 

Y (partial)  =  0.62  +  0.0289(  W)  -  4.3703  x  10  21(15  -  W)16  75. 
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.150 

1.1 


1.0 


.9 


.8 


-I-  .7 


.6 


.5 


.4 

.3  I  i  i  i  i  i  i  i 

0         2         4         6         8        10       12  14 

-W- 

Figure  57.— Elements  of  /  =  f  (W),  graphed. 

Table  22.— Descriptor  for  /,  the  left  end  and  central  segment 


I  a  +  b(W)       (15 -IV)      (Y1-Y1)  Overlay  t 

=  =  =  =  control  points  cX16 


W  Y1  Y1  X1  Y2  X1/14.5        Y2/0.124  Y2 

(D  (2)  (3)  (4)  (5)  (6)  (7)  (8) 


0.0 

0.400 

0.620 

15.0 

0.220 

0.219 

.5 

.510 

.634 

14.5 

.124 

1.00 

1.00 

.124 

1.0 

.581 

.649 

14.0 

.068 

.97 

.55 

.069 

2.0 

.687 

.678 

13.0 

.021 

.90 

.17 

.020 

3.2 

.710 

.712 

.11.8 

.002 

.81 

.02 

.004 

t  Although  not  shown  in  the  text  here,  the  curve  implied  by  the  overlay  points  matched  16<n<17  in  the  X"- 
Standards  very  well.  Iterative  checking  showed  n  =  16.75  to  perform  best  within  that  range.  Then,  c  =  0.124/(14.5) 
=  4.3703  x  10"  21 


The  Right  End. — For  this  portion  of  the  curve,  the  differences  (Y3)  from  the  straight 
line  and  associated  scaling  information  are  listed  in  table  23.  It  is  apparent  in  figure  58 
(right  side)  that  the  objective  curve  is  not  well  matched  by  any  simple  power  of  W.  The 
same  held  true  for  a  multiple-component  descriptor  that  was  subsequently  developed.  By 


0.62  +  0.0289  (W) 


"FINAL"  CURVE 
FROM  FIG.  9 


124  .002 
.021 
.068 
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Table  23.— Descriptor  for  /,  the  right  end 


t 

a  +  b(W)     (Y1  - 1)  Overlay 
=  =  control  points 


W  I  Y1  Y3  IV/15  Y3/0.150 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

8.4 

0.863 

0.863 

0.000 

0.56 

0.00 

9.6 

.886 

.897 

.011 

.64 

.07 

10.4 

.897 

.921 

.024 

.69 

.16 

13.2 

.904 

1.001 

.097 

.88 

.65 

15.0 

.904 

1.054 

.150 

1.00 

1.00 

t  The  overlay  curve  is  compared  to  the  X"-Standards  at  the  right  side  of 
figure  58. 


reorienting  the  exponential  curve  in  space,  however,  reversing  Wto  X\  =  (15  -  W)  and 
inverting  Y3  to  YA  =  (0. 1 5  -  Y3),  the  curve  becomes  convex-upward  with  increasing 
values  of  A'l.  Also,  it  is  at  least  roughly  matched  by  the  sigmoid  Standard,  n  =  8  and 
7=0. 108  (see  table  24  and  the  left  overlay  curve  in  fig.  58).  This  sigmoid  was  the  result 
of  an  interpolation  between  the  two  Standards  /=  0. 1  and  0.2  at  n  =  8. 


X/XP 

Figure  58.— Overlay  curves  compared  to  the  Standards. 
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Table  24.— Descriptor  (continued)  for  /,  the  right  end 


„  t  tt 
(15  -W)  (Y1-/)   (0.150  -Y3)        Overlay        de  K 
-  =  control  points  = 


w 

X1 

Y3 

Y4 

X1/15 

Y4/0.150 

Y4 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

0.0 

15.0 

0.150 

1.00 

1.00 

0.150 

8.4 

6.6 

0.000 

.150 

.44 

1.00 

.146 

9.6 

5.4 

.011 

.139 

.36 

.93 

.139 

10.4 

4.6 

.024 

.126 

.31 

.84 

.130 

13.2 

1.8 

.097 

.053 

.12 

35 

.053 

15.0 

.0 

.150 

.000 

.00 

.00 

.000 

f  From  table  23,  column  4. 

tt  e"K  using  X1,  n  =  8,  /  =  0.108,  XP1  =15.0,  VP  =  0.150,  and  the  intercept, 

a  =  0. 


Then,  the  right  end  was  specified  as: 

f4  =  (0.15-r3)  =  0.150(e-i9 
Y3  =  0.150 -0.150(e-*) 


73  =  0.150-0.150 


(15-  W) 
15 


-  1 


1-0.108 


—  e 


1 


1-0.108 


l-e 


1 

'  1-0.108 


and  simplified: 


Y3  =  0.1 6349  e 


(15-  W) 
15 


-  1 


0.892 


The  three  components  were  assembled  for  /. — 


/=  0.62  +  0.0289(  W)  -  4.3703  x  10~  21(15  -        75  -  0. 16349e 


(15-jn 
15 


0.892 


The  overall  performance  of  this  descriptor  car^  be  seen  from  the  tabulation  of  the  objec- 
tive and  estimated  /  in  table  25,  where  /  and  /  were  suitably  close  over  the  range  of  W, 
as  already  verified  by  segment  during  development. 
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Table  25.— The  performance  of  / 


-  THREE  DIMENSIONS 


Objective 

W  I  I 


0.0 

0.400 

0.401 

.5 

.510 

.510 

1.0 

.581 

.580 

2.0 

.657 

.658 

3.2 

.710 

.709 

5.0 

.766 

.764 

7.5 

.838 

.835 

8.4 

.863 

.859 

9.6 

.886 

.886 

10.4 

.897 

.900 

13.2 

.904 

.905 

15.0 

.904 

.904 

Summary  for  C 

At  this  point,  all  necessary  inputs  have  been  derived  and  the  complete  hypothesis  is: 


(X\/XPl)-l 


C  =  CP\ 


l-I 


l-e 


(1/(1  -I))2 


(1/(1  -Df 


where: 


C"     =  coverage  in  square  feet,  estimated  from  the  preleast  squares  fit  motel 

CPl  =  coverage  peak  for  any  specified  windspeed,  W,  and 

CP\  =  6390  +  90.739(15-  W)XA  +0.017894(15-  W)\  0<  W<  15 

X\     =  (1250-//) 

H     =  drop  height,  50<//<300 

XPI  =  point  in  X\  at  which  CPl  occurs,  or 

XPl  =  1190-0.4572(15-  W)2-25-  1.4811  x  10  22(15—  W)10 

I      =  inflection  point  in  X  expressed  as  a  proportion  of  XP,  (XI /XP  in  Matcha- 
curve-1,  Jensen  and  Homeyer  1970), 


/       =  0.62 +  0.0289(^)- 4.3703  x  10"2,(15  —  W)16  75  -  0. 16349<? 


(15  -  W) 
15 

0.892 


The  hypothesis  was  then  adjusted  back  to  the  data  from  which  it  was  partially  derived 
following  the  fitting  procedures  discussed  under  2-D,  X"  and  e~K,  where  C  =  X  and  the 
simple  least  squares  adjustment  coefficient  for  the  entire  hypothesis, 

b  =  LXY/LX2=  1.000569  and  C=bC . 
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The  coefficient,  b,  is  very  near  1 .0000;  so  it  is  evident  that  even  before  least-squares 
adjustment,  the  original  hypothesis  was  well  alined  with  the  data  spatially.  But,  in 

A 

checking  the  form  of  the  hypothesis  by  comparing  predicted  (C)  and  actual  (C) 
coverages  in  figure  59,  at  least  one  anomaly  appears.  Actual  coverage  data  points  for 
the  three  high-wind  points  (W=  8.4,  10.4  and  13.2)  probably  lie  too  close  to  the  surface 
since  substantial  variation  occurs  about  the  surface  at  lower  wind  levels;  i.e.,  reasonable 
variance  about  the  surface  would  be  expected  at  high-wind  levels  also. 

W  -  GROUP 


Figure  59.— Hypothesis  performance  by  wind  group. 

Recall  that  through  the  use  of  controlled  /-values  (fig.  56)  for  curves  passing  close  to 
the  three  high-wind  points,  low  variation  of  these  points  from  the  surface  was  assured. 
Alternative  /-values  could  be  used  that  would  allow  for  variation  from  the  surface,  in 
accord  with  some  assumed  variance  criteria.  This  refinement  was  not  attempted  and  the 
model  was  adopted  as  developed.  Predicted  values  (table  26)  and  the  associated  surface 
(fig.  60)  for  the  model  are  as  shown  below. 

Table  26.— Model  values 


Drop  height 


Wind  50  100  150  200  250  300 


mi/h 

 Feet 

0 

7,198 

7,804 

8,316 

8,710 

8,969 

9,081 

5 

6,830 

7,525 

7,707 

7,337 

6,491 

5,338 

10 

6,617 

6,707 

4,925 

2,619 

1,009 

282 

15 

6,345 

5,663 

3,460 

1,447 

414 

81 
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50  300 
H 

Figure  60.— The  3-D  hypothesis. 


It  is  appropriate  to  review  the  imposing  list  of  constraints  that  can  be  met  in  models 
like  this  one.  Controllable  items  include:  the  intercept  for  the  whole  surface;  the  eleva- 
tion, breadth,  and  positioning  of  the  peaks  over  the  interacting  independent  variables; 
the  trend  in  coverage  at  the  truncated  left  edge  of  the  surface;  the  magnitude  of  curva- 
ture above  and  below  the  inflection  points  in  the  bell-shaped  curves;  and  the  positioning 
of  the  inflection  points  themselves  within  the  range  of  the  independent  variables. 

Asymmetrical 
Bell-Shaped  Ridge 
over  XI,  with  Moving 
Peak  in  XI.  This  is 
Described  in  Two 
Segments  and 
Interacts  with  a 
Sigmoidal  ^-Effect 

Live  Data  Example  In  this  hypothesis  from  Klein,  Parker,  and  Jensen  1978,  tree  mortality  percent  in  a 

western  forest  is  characterized  by  tree  diameter  (d.b.h.)  over  the  course  of  a  beetle 
epidemic  (fig.  61). 

The  graphed  hypothesis  is  applicable  only  at  discrete  points  in  time  and  at  the  mid- 
points of  2-inch  d.b.h.  classes,  but  both  variables  are  treated  as  continuous.  As  was  evi- 
dent in  the  original  data,  the  more-or-less  bell-shaped  trends  over  time  differ  substan- 
tially on  either  side  of  the  central  ridge.  For  example,  observe  strong  asymmetry  at 
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1966  1967  1968  1969  1970  1971 

YEAR  OF  OBSERVATION 

Figure  61.— Tree  mortality  over  the  course  of  a  beetle 
epidemic. 


d.b.h.  =  12  to  18  inches,  made  most  noticeable  perhaps  by  elevational  differences  in  the 
curves  at  their  left  and  right  extremes.  See  also  the  dwindling  breadth  of  curve  crowns 
above  and  below  15  inches  d.b.h.,  an  interactive  change  included  in  the  descriptor  along 
with  asymmetry. 

To  assemble  descriptor  components,  left  and  right  sides  of  the  ridge  were  described 
separately  (fig.  62  and  63),  using  sigmoids  from  the  e~  A-Standards  that  varied  in  shape 
and  scale  according  to  data  trends  over  d.b.h.  Because  the  lower  portion  of  curves  of 
the  left  half  were  asymptotic  at  values  larger  than  zero,  it  was  necessary  to  include  the 
left-edge  intercept,  Int,  or  floor,  upon  which  the  left-side  sigmoids  rested  conceptually. 
Int  is  a  function  of  d.b.h.,  as  is  XP,  the  location  of  the  ridge  in  time. 


66      Yr  ^XP 


1966  1 967  1968  1969  1970  1971 


YEAR  OF  OBSERVATION 

Figure  62.— Tree  mortality  over  the  course  of  a  beetle  epi 
demic,  left  segment. 
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Sigmoids  are  described  within  the  range  XP±  5  years  ( Yr),  since  all  sigmoids  are 
estimated  to  be  fully  expressed  therein.  Use  of  a  constant  maximum  range,  such  as 
XP±  5  years,  keeps  descriptor  components  to  a  minimum  for  the  sigmoidal  time  effects, 
yet  permits  satisfactory  matching  of  the  objective  curves  in  this  case.  Restated  in  general 
terms,  the  whole  descriptor  is  as  follows: 


Mortality  percent  =  Int  +  Lsc(Lsig )  +  Rsc(Rsig ) 
where 

Lsig  =  left  sigmoid  =  fL(Ln,  LI,  XP,  Yr).  See  the  basic  sigmoid  parameters  defined  in 
Jensen  1979. 
Ln  =  Lsig  power  =  /,(d.b.h.) 

LI=  Lsig  inflection  point  as  a  proportion  of  the  range  in  years  from  XP-  5  to  XP, 
=/2(d.b.h.) 

XP=  point  in  time  where  surface  peaks  =/3(d.b.h.) 

Rsig  =  right  sigmoid  =  fR(Rn,  RI,  XP,  Yr) 
Rn  =  Rsig  power  =  /4(d.b.h.) 

RI  =  Rsig  inflection  point  as  a  proportion  of  the  range  in  years  from  XP+  5  to  XP, 

=/5(d.b.h.) 
Int  =  intercept,  left  edge  =  /6(d.b.h.) 
Rsc  =  Rsig  scalar  =  ridgetop  =  /7(d.b.h.) 
Lsc  =  Lsig  scalar  =  (ridgetop)  -  (Int ) 
Segmental  constraints: 

Left  side;     66  <  Yr <  XP,  discrete  values  only 

Right  side;    XP<  Yr<l\,  discrete  values  only 

Either  side;  8<d.b.h.<  18,  midpoints  of  2-inch  d.b.h. 
class  only  ...  8,  10,  12,  etc. 
The  model,  refitted  by  least  squares  to  smoothed  mortality  percent  at  36  control 
points  on  the  original  data  cross  sections  over  time  resulted  in  an  R2  of  0.96.  Used  as  a 
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goodness-of-fit  index,  this  high  /?2-value  attests  to  the  fact  that  the  descriptor  duplicates 
the  objective  graph  with  reasonable  accuracy. 

The  explicit  hypothesis  and  moderate  developmental  detail  follow: 

Mortality  percent  =  K{Int  +  Lsc[Lsig]  +  Rsc[Rsig]) 
where 

K=  least-squares  coefficient  =  0.9877 
Int=  l+0.2321(d.b.h.-8)163 
/&c  =  2.6429(d.b.h.)-5.157 
Lsc  =  Rsc-  Int 


5-\Yr-XP\ 


Lsig  = 


\-LI 


Ln 


-e 


Ln 


\-e 


Ln 


XP=  67.65+1. 2257 e 


20-d.b.h. 
12 
0.65 


-  1 


3.9 


7-|d.b.h.-13.5| 


-1 


17=0.851  -  0.1531<? 


0.45 


Ln=l.5 

Rsig  =  as  for  Lsig,  but  using  RI  and  Rn 
RI=  0.752  +  7.8173  x  10  4|d.b.h.  -  14|2-8 
Rn  =  2 

Limits 

66<  yr<71,  discrete  units  only 
66  <  Yr<XP,  for  Int,  Lsc,  Lsig 
XP<  Yr<l\,  for  Rsc,  Rsig 
8<d.b.h.<18. 


Supplementary  Int  was  estimated  from  leftward  extension  (from  the  ridge)  of  the  mortality  trend  in- 

Explanatory  Detail  dicated  by  data  in  each  of  six  2-inch  d.b.h.  groups.  All  six  trends  reached  lower  asymp- 

totes by  1963  and,  at  that  point,  suggested  a  flat,  concave-upward  curve  over  d.b.h. 
(fig.  64).  This  was  satisfactorily  described  by  using  the  ^"-Standards. 

Rsc  is  simply  the  height  of  the  ridge  above  zero  (no  intercept)  and  serves  as  the  scalar 
for  the  right-half  sigmoids.  Rsc  is  a  linear  function  of  d.b.h.  adopted  to  represent  the 
somewhat  irregular  pattern  of  ridge  values  for  the  six  d.b.h.  groups  (fig.  65).  Note  that 
the  ridge  line  in  figure  62  only  appears  to  be  sigmoidal  by  reason  of  the  sigmoidal 
change  of  point  -of-peaking  in  time  with  a  change  in  diameter. 
Lsc  =  (Rsc-  Int )  and  is  the  scalar  for  the  left-half  sigmoids. 

XP,  the  point  in  time  at  which  the  bell-shaped  curves  peaked,  was  estimated  from  the 
XP's  for  the  six  data-group  cross  sections  (fig.  66).  This  curve  was  estimated  to  asymp- 
tote at  68.88  for  lower  diameters  and  67.75  for  higher  diameters.  A  suitable  match  was 
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Figure  64.— Int  smoothed. 
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Figure  65.— flsc  smoothed. 
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Figure  66.— XP  smoothed. 
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found  in  the  e  ^-Standards  when  the  d.b.h.  scale  was  reversed  to  20-(d.b.h.). 

Lsig  specifies  the  sigmoidal  shape  of  the  left  side  (fig.  62).  After  transforming  Year  to 
5  - 1  Yr-XP\  to  create  A- values  ranging  from  zero  at  XP-  5  to  5  at  XP,  the  six  d.b.h. 
group  cross  sections  (left  halves)  were  each  scaled  to  1 .0  in  A1  and  Y  at  critical  points  in 
X.  Overlay  curves  were  plotted  for  these  cross  sections  (fig.  67).  The  pattern  was  one  of 
wider  sigmoidal  crowns  near  d.b.h.  =  14  inches,  narrowing  with  departure  from  that 
d.b.h.  This  is  reflected  in  the  bell-shaped  function  for  the  inflection  points,  LI.  Also, 
the  e~A-Standards  with  n  =  1.5  were  found  to  represent  this  curve  array  with  reasonable 
accuracy,  so  Ln  was  set  at  the  constant,  1.5. 


1.0  r 


X/XP 

Figure  67.— The  overlay  curves. 

LI,  the  function  for  the  inflection  points,  is  the  dashed  line  in  figure  68  and  repre- 
sents the  inflection  points  adopted  and  plotted  for  each  overlay  curve.  LI  reflects  the 
width-of-curve-crown  trend  noted  under  Lsig  and  in  the  text. 

.9  r 


.6  I  i  i  i  1  1  1  1 

8       10       12      14       16      18  20 

D.B.H. 

Figure  68.— LI  smoothed. 
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Asymmetrical 
Bell-Shaped  Ridge 
over  XI,  with 
Moving  Peak  in  XI. 
This  is  described  in 
Three  Segments  and 
Interacts  with  a 
Sigmoidal  ^2-Effect 

Live  Data  Example 


Rsig,  RI,  and  RN  were  obtained  in  a  similar  fashion.  The  Year  transform  again  was 
5  -  |  Yr-XP\  and  ranged  in  value  from  5.0  at  XP  to  zero  at  XP+  5.  The  right-side  cross 
sections  were  fairly  well  matched  by  e~  ^-Standards  with  n  =  2,  so  Rn  was  set  constant  at 
Rn-2.  RI  showed  the  same  general  widening  of  curve  crown  at  d.b.h.  =  14  inches. 


Quality  score  (0  was  originally  described  as  an  aggregate  of  1 1  planar  regions  over 
flow  {F)  and  stability  (5)  of  gap-graded  road  materials,  a  strongly  segmented  descriptor 
(fig.  69,  from  Lee  and  others  1973).  The  objective  here  was  to  smooth  the  figure  mathe- 
matically. Although  the  problem  is  an  unusual  one,  its  solution  serves  admirably  to 
demonstrate  descriptor  segmentation. 

From  figure  69,  it  can  be  seen  that  opposite  sides  of  the  ridge  differ  substantially  in 
slope  so  that  an  symmetrical,  flat-topped,  bell-shaped  curve  would  be  required  to  de- 
scribe the  cross  section  at  any  point  in  S  while  rounding  the  junctures  of  planes.  A  rela- 
tively simple  descriptor  alternative  involves  segmentation  of  the  (S)  and  (F)  regions  as 
shown,  with  left  and  right  orientation  lines.  Note  that  the  lines  are  parallel  to  their 
respective  sides  of  the  ridge  and  lie  one  unit  in  Fcloser  to  the  ridge  center.  This  permits 
use  of  a  single,  but  different,  sigmoid  cross  section  to  represent  and  smooth  the  corners 
of  each  segment,  left  and  right.  The  sigmoids  will  peak  at  their  respective  orientation 
lines  and  will  be  functional  over  a  constant  distance  from  them  (left  or  right  as  appro- 
priate). The  constant  will  vary  by  side  as  needed.  Since  the  ridge  is  flat  topped,  it  will 
have  a  value  of  t=  10  everywhere  in  the  center  segment. 
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Figure  69.— Quality  of  gap-graded  road  materials:  segmented 
planar  form. 
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Then,  to  this  point,  we  have  described  a  ridge  with  three  segments  specified  in  terms 
of  S  and  F;  the  sigmoidal  effect  to  the  left  of  the  left  orientation  line;  the  sigmoidal  ef- 
fect to  the  right  of  the  right  orientation  line;  and  the  flat  ridge  area  between  these  lines 
at  a  value  of  /=  10. 

Finally,  a  sigmoidal  truncation  of  the  front  end  of  this  ridge  is  achieved  through 
multiplication  of  all  components  of  the  descriptor  by  an  appropriate  sigmoid,  changing 
in  value  from  zero  to  one  within  the  range  0<S<2000,  and  being  applicable  for 
0<S<5000. 

In  general  terms,  we  have: 
Lsig  and  Rsig  =  left-  and  right-segment  sigmoids,  respectively 
LO,  RO  =  left-  and  right-orientation  lines,  respectively 
CC=  center  segment,  constant 
Tsig  =  truncator  sigmoid 

Then 

t  =  Tsig(Lsig  +CC+  Rsig ) 
and 

Tsig  =  AS) 

Lsig=J{LO,F) 

Rsig=J{RO,F) 

LO=AS) 

RO=f{S) 
Limits 

For  Lsig,    LO<F<  28 

Forfts/g,  0<F<RO 

ForCC,  LO<F<RO 
and,  0<S<5000. 

The  final  descriptor  form  is  shown  in  figure  70  along  with  the  original  planar  form 
for  comparison.  It  can  be  seen  that  the  descriptor  does  a  creditable  job  of  emulating  the 
planar  form  while  smoothing  the  corners,  all  with  three  segments  in  place  of  the  original 
1 1  segments. 


FLOW.  (F) 

Figure  70.— Quality  of  gap-graded  road  materials:  segmented 
planar  and  final  descriptor  forms. 
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Explicit  hypothesis  and  moderate  developmental  detail  follow: 

If 

F>LO,  t=  Tsig(Lsig) 
LO<F< RO,  t=  WTsig) 
F<RO,t=Tsig(Rsig) 

where 

F=  soil  flow 

LO=  8.888  +  0.0017036<S) 
RO  =  8.404  +  0.00048076(S ) 
S  =  soil  stability 


Lsig=  We 


(LO  +  20-F) 
20 


-1 


0.23 


2.2 


Tsig  =  1  -  e 


(6000-5) 
6000 
0.185 


-  1 


Rsig=  10  e 


\RO-U\+F 


11 


0.30 


2.4 


Limits 
0<F<28, 
0<S<5000. 


Supplementary  Lsig,  the  left  segment  of  figure  69,  has  a  constant  cross  section  over  F  for 

Explanatory  Detail  1240  <5<  5000,  approximated  in  the  descriptor  by  a  sigmoid  oriented  at  F=  LO  and 

extending  leftward  a  distance  of  20  units  in  F.  This  is  about  the  minimum  operational 
span  for  the  sigmoids  here.  At  the  upper  extreme  of  S,  (S  =  5000),  the  left-segment 
sigmoid  is  expected  to  be  completely  specified  in  the  range  27  >  F>  LO,  as  in  figure  70. 
Setting  28  as  the  upper  level  of  F  within  which  all  left-segment  sigmoids  must  be  com- 
pletely specified,  we  turn  to  the  limiting  case  at  S  =  0.  Here,  LO  =  8.888  and  the 
sigmoidal  range  must  then  be  at  least  (28  -  8.888)  =  19. 1 12;  so  20  was  adopted  as  the 
operational  span.  Note  that  a  larger  span  could  have  been  adopted. 

The  F-scale  is  reversed  to  (LO  +  20  -  F),  as  shown  in  figure  7 1 ,  to  associate  the  largest 
value  of  the  sigmoidal  span,  20,  with  the  peak  of  the  sigmoid  at  LO.  By  scaling  control 
points  from  the  left  cross  section  at  S  =  1240  (fig.  71)  to  1 .0  in  X,  (LO  +  20  -  F),  and  Y, 
(score,  t),  and  making  an  overlay  curve,  an  appropriate  sigmoid  was  identified  from  the 
e~  ^-Standards. 


71 


MATHEMATICAL  DESCRIPTION  OF  GRAPHED  HYPOTHESES  -  THREE  DIMENSIONS 
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Figure  71  .—A  t-cross  section  over  F  at  S  =  5000.  Note:  the 
cross-sectional  shape  of  the  left  segment  is  constant  in  the 
range  1 240  <  S  <  5000;  the  right  segment,  although  different 
in  cross-sectional  shape,  is  also  constant  in  shape  over  the 
range  1240  <S<5000. 


Tsig,  from  (fig.  72),  the  planar  ridge  truncation  ranges  from  /  =  0  at  S  =  840  to  t  =  10 
at  S=  1240.  This  plane,  scaled  to  1.0,  is  matched  and  smoothed  by  Tsig,  as  is  shown  in 
figure  72.  Note  from  the  explicit  Tsig  formula  that  the  inverted  sigmoid  to  the  left  of 
the  truncation  plane  (with  base  =  1 .0  and  peaking  at  0,  0)  was  described  on  the  reversed 
S-axis,  (6000-  S),  and  subtracted  from  1.0  to  arrive  at  Tsig.  This  provided  a  more  ac- 
curate duplication  of  the  truncator  plane  than  did  other  sigmoid  alternatives  in  this  case. 

Although  a  maximum  of  S  =  6000  appears  in  Tsig,  the  applicable  range  is  still  limited 
to  5  =  5000  based  on  the  original  figure  (fig.  69). 

PLANAR  RIDGE 


TRUNCATION  SCALED  TO  1.0 


Figure  72.— Ridge  truncation  sigmoid,  Tsig. 

Rsig,  the  right  segment  of  figure  69,  has  a  constant  cross  section  over  F  for 
1240  <  5<  5000,  approximated  in  the  descriptor  by  a  sigmoid  oriented  at  RO  and  ex- 
tending to  the  right  for  a  distance  of  1 1  units  in  F.  Eleven  is  about  the  minimum  opera- 
tional span  for  the  right-segment  sigmoids.  RO  ranges  from  8.404  at  S  =  0  to  10.808  at 
5  =  5000.  A  range  of  1 1  includes  F=  zero  at  both  extremes  and  so  was  adopted.  The 
Rsigs  were  described  as  a  function  of  the  F-transform  \RO  -  1 1 1  +  F,  as  shown  in  figure 
71 .  Thus,  the  maximum  value  of  1 1  always  occurred  at  the  sigmoidal  peak,  F=  RO,  and 
Rsig  functioned  over  the  range  zero  to  1 1  of  the  F-transform. 

Summary  The  final  surface,  then,  is  simply  the  sum  of  the  three  contiguous  ridge  segments,  all 

of  which  are  truncated  at  appropriate  points  in  S  through  multiplication  by  the  propor- 
tional values  of  Tsig. 
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Symmetrical,  Inverse 
Bell-Shaped  Ridge 
over XI,  with 
Constant  Point  of 
Minimum  in  XI  and 
Interacting  with  a 
Concave-Upward 
^2-Effect 


Live  Data  Example  Here,  the  amount  of  natural  reproduction  (number  of  lodgepole  pine  seedlings  per 

acre)  over  the  Sleeping  Child  burn  was  expected  to  be  affected  by  both  azimuth  and 
slope  (data  from  Lyon  1976).  Further,  azimuth  and  slope  were  expected  to  interact 
strongly,  the  magnitude  of  response  to  azimuth  being  dependent  on  slope  percent. 

A  small  data  set  was  used  to  develop  a  graphic  hypothesis  consistent  with  these  expec- 
tations. It  was  very  close  to  the  mathematical  hypothesis  pictured  in  figure  73. 


0         §0      180      226  360 

AZIMUTH  FROM  NORTH, (DEGREES) 

Figure  73.— The  mathematical  hypothesis. 


Descriptor  Development 
for  Figure  73 


The  Azimuth  (A  )-effect  may  be  seen  as  cross  sections  over  AZ  at  slope  percent 
(SL)  =  20,  30,  35,  and  35.  These  were  screened  for  shape  using  the  e~A-Standards.  Note 
that,  since  the  bell-shaped  curves  were  symmetrical,  either  sigmoidal  half  would  be 
representative  of  the  whole  (the  left-half  was  used).  Also,  negative  departures  from 
either  edge  could  be  used  to  specify  the  depth  of  the  sigmoids  (these  negative  values 
were  ignored  in  the  screening,  but  were  accounted  for  in  the  final  mathematical 
hypothesis).  Sigmoids  were  found  to  be  sufficiently  similar  at  different  SL-levels  to 
adopt  the  common  form, 


e~K=  1.0675*6- 


AZ 
180 


0.655 


2.4 


0.0675, 


where  the  trough  is  at  a  fixed  point  in  AZ  (XP=  180,  here). 
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And,  given  the  subtraction  of  suitably  scaled  sigmoids  from  either  edge  (EDGE), 


AZ 
180 


2.4 


Number  of  seedlings  =  EDGE  -  rP(1.0675*e 


0.655 


0.0675), 


where  EDGE,  was  better  matched  by  the  lower  85  percent  of  a  sigmoid  on  the  reversed 
slope  scale  with  a  maximum  value  of  100  than  by  a  simple  A^-form.  Thus, 


EDGE  =  2400*  e 


(100- SL) 
100 
0.255 


-  1 


2.6 


The  trough  of  the  objective  figure  (MIN)  was  expressed  as  the  height  of  the  negative 
bell-shaped  (or  sigmoidal)  peaks  above  the  3-D  base.  As  with  EDGE,  this  was  described 
as  the  lower  85  percent  of  a  sigmoid  on  the  reversed  slope  scale  with  a  maximum  value  of 
100.  And, 


(100-  SE) 
100 

MIN  =  250*  e  °'29 
Then,  YP=  EDGE  -MIN 
LIMITS:  0<v4Z<360,  20<SLOPE<45. 


-  1 


2.3 


Note  that  inverse  asymmetrical  bell-shaped  functions,  often  typical  of  azimuth  effects  on 
plant  growth,  can  be  described  roughly  as  above,  but  in  two  segments.  Each  segment  then 
has  its  own  sigmoidal  shape.  Also,  it  could  easily  be  that  the  trough,  which  was  at  a  fixed 
point  in  AZ  (XP=  180),  in  the  foregoing  example,  changes  in  AZ  depending  on  slope.  In 
such  a  case,  XP=J{SL)  in  the  sigmoidal  effect  of  AZ. 


Concave-Upward 
Relation  in  Two 
Dimensions,  with 
Constant  Point  of 
Peaking  in  XI 
Over  X  2. 

Live  Data  Example  Soil  temperatures  (ST)  at  various  depths  and  percent  of  duff  reduction  were  studied 

on  a  clearcut  area  that  had  been  subjected  to  a  controlled  burn.  A  data  set  (from 
Shearer  1974  and  1975)  was  devoted  to  hypothesis  development  as  shown  below. 
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Data  were  grouped  into  duff-reduction  (DR )  classes  and  trends  over  "depth  in  soil" 
(DS)  were  smoothed  in  accord  with  expectation.  These  were  screened  over  the  reversed 
DS-scale  in  the  A^-Standards  and  were  all  reasonably  represented  by  &,(7  -  DS  )6. 
Then, 

ST=  45  +  ( rP/(7)6)(7  -  DS  )6. 
And,  YP,  the  smoothed  curve  of  scaling  heights  for  the  individual  Z)5-curves  (at  the 
scaling  point  [7  -  DS]  =  7),  was  concave-upward  over  DR  and  could  have  been  described 
as  an  exponential  function  of  DR.  The  lower  half  of  a  sigmoid,  however,  was  found  to 
mimic  the  objective  curve  more  closely  so, 


DR 
200 


2.2 


YP  =  365.  ll\(e 


0.56 


)-  10.171, 


thus  producing  figure  74.  Limits  of  use  for  this  hypothesis  are:  0<DR<  100  and  0<DS<7. 


O  iu 
Z  u 


DUFF 
REDUCTION, % 


DEPTH  IN  SOIL, (CM) 

Figure  74.— The  mathematical  hypothesis. 


Note  that  only  the  lower  (left)  half  of  the  FP-sigmoid  is  utilized  in  this  relation  as  the 
result  of  doubling  the  DR-scale  (XP=20Q)  and  the  S T  scale  as  well  ( YP  +  intercept  =  400,  or 
YP=  355).  In  the  last  case,  it  would  have  been  more  consistent  to  have  used  YP  =  400,  but 
a  suitable  e~K  was  found  with  355;  so  355  was  retained. 
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MATHEMATICAL  DESCRIPTION  OF  GRAPHED 
HYPOTHESES  -  FOUR  DIMENSIONS 
Outline  for  "FOUR  DPMENSIONS" 


Asymmetrical  Bell-Shaped  A'l-Effect  with  Point  of  Peaking  Moving 
in  XI.  The  Al-Effect  is  Described  in  Two  Segments  and  Interacts 
with  the  Bell-Shaped  A2-Effect  and  the  Concave-Upward  A3-Effect 
(Contrived  Example)   


Y 


LO  MED  HI 

XS 


Asymmetrical  Bell-Shaped  A'l-Effect,  with  Constant  Point  of 
Peaking  in  XI.  The  A'l-Effect  is  Described  in  Two  Segments  and 
Interacts  with  Concave-Upward  XI-  and  A3-Effects 
(Live  Data  Example)   83 


+ 

X1 


Convex-Upward  Al-Effect  with  Peak  Moving  in  A2,  Interacting 

with  Concave-Upward  A2-Effects  and  Slight  Convex-Upward 

A3-Effects  (Live  Data  Example)   86 


X3 


0 

X1 
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PAGE 

Concave-Upward  Al-Effects,  with  Constant  Point  of  Peaking  in  X2 
and  Interacting  with  Convex-Upward  A^2-Effects;  Fourth  Dimension 
Discrete  (Live  Data  Example)   90 


Y 


X1 


77 


MATHEMATICAL  DESCRIPTION  OF  GRAPHED  HYPOTHESES  -  FOUR  DIMENSIONS 


Asymmetrical 
Bell-Shaped  JQ-Effect 
with  Point  of  Peaking 
Moving  in  X2.  The 
.Yl-Effect  is 
Described  in  Two 
Segments  and 
Interacts  with 
the  Bell-Shaped 
JT2-Effect  and  the 
Concave-Upward 
^3-Effect 


Contrived  Example  Assume  that  finite  curve  forms  and  scales  for  a  4-D  relation  have  been  developed 

(within  gross  constraints  of  the  expected  relation)  from  trend  information  contained  in 
half  of  a  current  data  set,  partitioned  over  the  ranges  of  the  three  independent  variables 
The  most  prominent  feature  of  the  resulting  curve  array  is  the  asymmetric  nature  of  the 
AT-effect.  This  is  displayed  in  figure  75  at  three  levels  of  X2,  within  each  of  three  levels 
ofX3. 


X3 


0       10    20     30      0      10     20     30      0      10     20     30  (X1) 
50      60    70      80     50     60     70     80     50     60     70     80  (X1  +  50) 

Figure  75.— The  partial  graphed  hypothesis. 


Also  prominent,  is  the  ^-effect,  bell  shaped,  but  truncated  at  the  high  end  of  the  XI 
scale.  This  is  perhaps  more  fully  appreciated  in  figure  76. 


Y 


0     10     20    30  0     10    20    30  0     10     20  30 

X1  X1  X1 

Figure  76.— The  complete  graphed  hypothesis. 

The  number  of  X2-  and  A'3-partitions  are  kept  to  a  minimum  to  simplify  the  presen- 
tation of  associated  descriptive  methods.  More  partitions  would  probably  be  desirable  in 
a  live  data  application. 


MATHEMATICAL  DESCRIPTION  OF  GRAPHED  HYPOTHESES  -  FOUR  DIMENSIONS 


Descriptor  methods  involved  are  identical  to  those  shown  earlier  for  2-  and  3-D 
hypotheses  and  are  simply  expanded  to  accommodate  another  independent  variable. 
From  figure  76  we  can  anticipate  a  reasonable  mathematical  hypothesis  to  be  as  shown 
below. 

AT-Effect. — Since  the  bell-shaped  curves  are  asymmetric,  it  would  be  efficient  to 
create  a  sigmoidal  descriptor  for  the  left  and  right  sides  separately.  They  would  be 
joined  at  XPl  and  would  be  scaled  to  a  common  height,  YPl.  Thus,  with  XI  adjusted 
to  X\  +  50  to  accommodate  the  curve-breadth  limit,  XP± XP,  the  hypothesis  for  say, 
the  left  half,  would  be: 


YL  =  YPl{e~K)  left 


or, 


YL  =  YP\ 


(X\  +  50) 


XPl 


l-IL 


nL 


nL 


l-IL 


nL 


l-IL 


where  YP\,  XPl,  IL,  and  nL  may  change  over  the  range  of  XI,  X3,  or  both. 

YP\. — Plotting  the  heights  of  the  nine  curves  from  figure  76  in  figure  77  and  smooth- 
ing over  these  points,  YPl  varies  as  a  symmetrical  bell-shaped  curve  (er  A  )  over  XI, 
changes  in  height  ( YPl),  and  perhaps  changes  in  shape  (n  and  /)  over  X3.  Then, 
YP1  =  YPl(e^K). 

X3 


YP1 


0         20        40        60         80       100  (X2) 

10        30         SO        70         00       1  10  (X2  +  10) 

Figure  77.— YP1  smoothed. 

YPl. — From  figure  78,  the  plotted  heights  of  these  curves  appear  to  follow  an  ex- 
ponential function  of  X.  Let  us  assume  that  we  have  made  and  compared  an  overlay 
curve  with  the  ^"-Standards  and  that  we  have  found  a  suitable  descriptor  in 


0  200     400      600  800 

X3 

Figure  78.— YP2  smoothed. 
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XI  in  e2  K.— The  largest  sigmoid  of  the  set  (for  A3  =  600)  reaches  X2  =  90.  Since 
XP2  =  45  here,  all  sigmoids  fall  within  the  required  limits,  XP±  XP.  But,  it  might  be 
reasonable  to  extrapolate  to  X3  =  700  where  the  expected  maximum  might  be  say  100. 
In  which  case,  a  constant  can  be  added  to  X2  in  e~K  to  accommodate  the 
(AT  ±  AT  )-requirement,  using  (X2  +  10)  in  place  of  X2,  AT2  =  55,  and  2XP2  =  1 10. 
Then,  the  broadest  curve  will  fall  within  XP2  ±  XP2. 

XP2  in  e2K. — As  may  be  seen  in  figure  77,  AT2  =  45  in  X2  for  all  levels  of  X3.  And, 
since  (A2  +  10)  is  to  be  used  in  place  of  X2,  XP2  =  55. 

n2  in  e2K. — Assume  that  overlay  curves  have  been  made  for  the  left  half  of  the  three 
symmetrical  bell-shaped  curves  in  figure  77  and  that  comparison  with  the  e~  ^-Standards 
resulted  in  n  =  1.7,  2.0,  and  2.5  for  XI  =  200,  400,  and  600,  respectively.  Without  the 
formality  of  plotting,  n2  will  be  assumed  to  be  an  exponential  function  of  A3  so  that: 
n2  =  a4  +  b4(X3)n*. 

I2  in  e2  K. — Along  with  n2,  above,  assume  that  appropriate  7-,-values  have  been  iden- 
tified as  0.50,  0.45,  and  0.40  for  A3  =  200,  400,  and  600,  respectively/By  inspection,  I2 
is  a  negative  linear  function  of  A3  so  that  I2  =  a5-  b5(X3). 

In  summary, 


JT1  =  YP2{ 


(A3 +10) 
XP2 


-  1 


l-L 


I- 1, 


where, 

YP2  =  b3(X3)"i 
XP2  =  55 
n2=a4+b4(X3)"4 
I2  =  a5  +  b5(X3). 

XI  in  (e~K)  )ef(. — Revision  of  XI  is  necessary,  since  the  smallest  XP  (from  figure  76) 
approaches  zero  and  the  broadest  curve  reaches  XI  =  35  at  X3  =  600.  Assuming  ex- 
trapolation to  X3  =  700,  we  might  expect  the  broadest  curve  to  reach  no  more  than  say 
A3  -  50. 

To  meet  the  constraint  that  all  curves  of  the  set  be  within  XP  ±  XP,  (X\  +  50)  was 
used  in  place  of  A'l .  Then,  where  ATI  =  zero,  (XI  ±  50)  =  50  and,  at  the  high  extreme 
where  A'l  =  35,  (A'l  +  50)  =  85.  And,  using  (A'l  +  50)  in  e~K,  all  curves  will  be  within  the 
minimum  ATI  =  50,  since  ±  2AT1  =  100. 

XP1  in  (e~K)  ,eft. — ATI -values  read  from  figure  76  are  plotted  and  smoothed  in 
figure  79. 

X3 


0         20        40  60 

X2 

Figure  79.— XP1  smoothed. 
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We  will  assume  by  inspection,  that  XPl  is  an  exponential  function  of  XI  and  that 

XPl  =  (YP3/(60)n<>(X2)"(>. 

YP3. — At  the  scaling  point,  X2  =  60,  the  heights  (7P3)  of  the  three  curves  in  figure  77 
are  read  as  11,  14,  and  18,  for  A3  =  200,  400,  and  600,  respectively.  Plotting  these 
points  in  figure  80  and  smoothing,  YP3  can  be  seen  to  be  a  very  flat  exponential  func- 
tion of  X3  so  that 

YP3  =  a  +  b1(X3)"\ 


20 


YP3  10 


O  200  400 

X3 

Figure  80.— YP3  smoothed. 


600 


Also,  we  will  assume  that  the  shapes  of  the  three  curves  in  figure  77  have  been  com- 
pared with  the  A^-Standards  and  have  been  found  to  have  a  constant  n. 

Then,  XPl  =  ( YP3/(60)"<>)(X2)n(> 
where  YP3  =  a1  +  b1(X3f . 

nL  and  IL  in  (e  K )  kft. — We  will  assume  that  the  left  sides  of  the  nine  curves  in 
figure  75  have  been  screened  for  matching  e  ^-Standards  with  the  results  shown  in 
tables  27  and  28. 


Table  27.— nL  in  (e~KL 


X3 


X2 


200 


400 


600 


20 
40 
60 


8J 
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Table  28.— IL  in  (e~K) 


left 


X2 


200 


X3 


400 


600 


20 
40 
60 


0.970 
.965 
.973 


IL 


0.965 
.959 
.955 


0.938 
.943 
.940 


average 


0.969 


0.960 


0.940 


nL. — Referring  to  table  27,  the  set  of  e~  A-Standards  with  n  =  2  was  deemed  ap- 
propriate for  the  set,  all  levels  of  X2  and  X3  included.  So,  NL  =  2,  a  constant. 

IL. — The  /-values  in  table  28  are  quite  large  since  the  left-side  curves  have  very  steep 
slopes  and  are  proportionally  far  to  the  right  in  the  (X\  +  50)  range.  Systematic  varia- 
tion is  associated  primarily  with  changes  in  X3  as  indicated  by  the  average  /L-values. 
These  are  plotted  and  smoothed  in  figure  81. 

1.00  i— 


IL 


®.®@ 


200      400      600  800 
X3 

Figure  81. —IL  smoothed. 

And,  again  assuming  this  relation  can  be  satisfactorily  expressed  by  an  exponential, 
we  have: 

IL  =  0.971  -b^Xl)"1. 
Then,  to  summarize,  the  entire  left  side  of  the  hypothesis  may  be  specified  as: 


or, 


YL  =  YP\ 


YL=YP\{e~K)  left 


(XI  +  50) 


nL 


XPl 


-  1 


nL 


\-IL 


l-e 


\-IL 


1 

\-IL 


nL 
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where  YPl,  XPl,  IL,  and  nL  may  change  over  the  range  of  XI,  A3,  or  both. 


FP1  =  YP2y 


(X2+  10) 
XP2 
1-/, 


1-/, 


l-e 


and 

YPl  =  fe3(A3)n3 

XP1  =  55 

n2  =aA  +  b4(X3)"4 

/,  =  a$  +  b,(X3) 

XPl  =  (YP3/{60)n*)(X2)"<>  and 

yP3  =  fl7  +  67(A3)"7 

nL  =  2 

/l  =  o.97i  -bpcyfi. 

Description  of  the  right  side  of  the  hypothesis  would  be  approached  in  the  same  fashion 
as  above,  to  provide  the  unknowns  in: 


YR  =  YPl 


X 
XPl 


l-IR 


nR 


-e 


\-IR 


nR 


1 


\-IR 


NR 


Note  that  XI  for  the  right  side  curves  could  be  set  =  (50  —  A'l)  and  YPl  and  XPl  are  the 
same  for  the  right  side  as  for  the  left. 

Asymmetrical 
Bell-Shaped 
Al-Effect,  with 
Constant  Point  of 
Peaking  in  X2.  The 
Al-Effect  is  Described 
in  Two  Segments  and 
Interacts  with 
Concave-Upward 
XI-  and  A3-Effects 

Live  Data  Example  The  relation  between  total  terpene  content  of  lodgepole  pine  phloem  by  d.b.h.  (D ) 

phloem  thickness  (P),  and  5-year  radial  growth  (G)  (from  W.  E.  Cole  and  others  1981) 
was  hypothesized  graphically  very  nearly  as  in  figure  82. 

In  explaining  descriptor  development  here,  we  will  focus  on  generalities.  The  details 
of  formulating  the  graphed  hypothesis,  selecting  control  points,  making  overlay  curves 
and  comparing  them  to  the  Standards  are  assumed  to  be  known  from  earlier  text. 
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-  FOUR  DIMENSIONS 


0  4  8  12  16  20  IN 

0  10  2  20.3         30. 5         40.6         50.8  CM 


D.  B.H. 

Figure  82.— The  mathematical  hypothesis. 

In  this  example,  the  phloem  thickness  (P)  effect  was  strong  in  both  expectation  and 
data  trends  (after  testing  more  modest  hypotheses,  a  data  set  was  exploited  for  inter- 
active hypothesis  development).  Its  effect  in  the  data  was  clearly  exponential  within 
d.b.h. /growth  data  groups  and  appeared  to  be  reasonably  constant  at  about  n  =  1.55. 
Then,  where  the  scaling  point  in  P  was  at  G  =  0. 14  (the  approximate  upper  limit  of  the 
data  over  G), 


percent  terpene  content  = 


YPP 
(0.14) 


1.55 


,1.55 


The  heights  of  the  peaks  (YPP)  of  the  asymmetrical  bell-shaped  curves  over  d.b.h.  (D),  at 
the  constant  point  of  peaking  (XP),  10.5  inches,  are  specified  as: 


for  D<10.5 


YPP=  YPD(\.O0\65e 


D  +  8.5 
19 


-  1 


0.395 


0.00165), 


for  D>10.5 


YPP=  YPD(\.01092e 


38.5-D 
28 


0.490 


1.4 


0.07092), 


where 


YPD  =  0.38+  1.0292(G)14. 


Note  that  sigmoids  over  D,  while  having  a  different  shape  on  opposite  sides  of  the 
peak,  were  identified  as  having  consistent  shapes  for  the  three  growth  levels. 
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This  completes  the  formulation  applicable  in  the  range  0<P< 0.20,  0<D<25,  and 
0<G<0.35,  but  explanation  of  the  D-transforms  for  the  sigmoids  may  be  helpful  here 
The  broadest  bell-shaped  curve  expected  (at  P=  0.20)  ranges  from  about  D=  -  8.5  on 
the  left  side  to  D=  38.5  on  the  right.  With  XP=  10.5,  the  D  scale  for  the  left  side  is 
changed  to  (D+  8.5)  and  that  for  the  right  is  (38.5  -D),  as  shown  in  figures  83  and  84. 


-12  -8      -4     0      4       8*12  (O) 


0  8.S  19.0  (D  +  8.5)=X1 

Figure  83.— D-transform  for  the  left  side. 


•    |12     18    20     24    28    32   38    |40  (D) 

28.0  0  (38.5-D)=X1 

Figure  84.— D-transform  for  the  right  side. 


Again,  these  transforms  of  D  insure  that  all  points  of  the  sigmoids  are  within  the 
range  XP\±XP\. 

The  foregoing  hypothesis  was  ultimately  refitted  to  the  data  from  which  it  was  largely 
derived.  This  provided  the  most  realistic  appraisal  of  scale  available  at  the  moment  for 
the  relation  characterized. 

Given  that  total  terpene  estimated  from  the  model  =  X  and  that  terpene  content  ac- 
tually measured  in  the  data  =  Y,  the  least  squares  coefficient,  b,  was 

b  =  LXY/LX2  =  1.10301. 

This  corrects  for  underscaling  in  the  hypothesis  insofar  as  the  data  set  utilized  is  con- 
cerned. Slight  improvement  in  R2  was  noted  for  the  interactive  hypothesis  over  that  for 
the  simple  additive  model  with  only  linear  effects  of  D,  P,  and  G:  R2  =  0.39  versus  0.35, 
respectively.  So,  no  potentially  great  advantage  of  the  interactive  over  the  additive 
model  appears  to  exist  in  this  data  set.  But,  the  hypothesis  still  provides  new  insights  for 
the  area  of  interest  that  should  be  of  value  in  planning  new  studies. 
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Convex-Upward 
JTl-Effect  with  Peak 
Moving  in  X2, 
Interacting  with 
Concave-Upward 
^2-Effects  and  Slight 
Convex-Upward 
^3-Effects 

Live  Data  Example  Laboratory-controlled  calibration  data  for  Peltier  thermocouple  psychrometers  were 

assembled  by  Brown  and  Bartos,  1981,  and  are  shown  smoothed  over  water  potential  in 
figure  85. 

COOLING  TIME, SECONDS 


I  I  I  I  i  i  i  i  i  i  i  i 

-80     -60     -40     -20        0       -80     -60     -40     -20        0       -80     -60     -40     -20  0 

WATER  POTENTIAL, BARS 

Figure  85.— The  graphed  hypothesis  and  model  performance. 


The  combination  of  strong  expectations  for  the  relation,  a  well-conceived  laboratory 
study,  and  strong  data-group  means  left  little  doubt  as  to  the  form  and  scale  of  the  rela- 
tion between  microvolt  output  and  the  independent  variables,  water  potential  in  negative 
bars,  temperature  of  the  measured  medium  and  the  cooling  time.  The  smoothed  curves 
represent  an  exceedingly  strong  hypothesis  and,  with  one  additional  variable,  formed  the 
basis  for  a  published  mathematical  model  intended  for  adjustment  to  other  psychrom- 
eters of  the  same  design.  The  four  dimensional  relation  pictured,  in  fact,  was  described 
mathematically  before  including  the  effects  of  the  fifth  dimension  and  provides  a  unique 
opportunity  to  study  four  dimensional  descriptor  methods. 

The  strategy  here  was  to  describe  microvolt  output  (MW)  as  a  function  of  water 
potential  ( WP),  temperature  (T),  and  cooling  time  (CT),  solving  the  resulting  equation 
for  (WP). 

It  was  noted  at  the  outset  that  interest  in  WP  was  limited  to  XP<  WP<0,  where  XP 
is  the  point  in  WP  where  microvolts  peak.  Also,  it  can  be  seen  in  figure  85  that  XP 
changes  over  both  rand  CT.  The  heights  of  the  curves,  called  upper  intercepts  {(JI) 
here,  change  with  both  Tand  CT,  and  all  curves  over  WP  are  oriented  at  zero. 
The  general  form  of  the  relation  is  estimable  from  figure  85  as: 

MV=  UI\-(UI\/\XP\n\)\XP-  WP\nK 
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This  is  an  exponential  in  the  reversed  ^P-scale,  \XP—  WP\,  scaled  to  the  difference  be- 
tween UI\  and  zero  (=  UI\)  at  W/P  =  zero,  and  subtracted  from  UI\.  Smoothed  values 
of  AT  and  UI\  were  established  and  described  first,  since  they  were  subsequently  used 
in  place  of  actual  XP  and  UI\ -values  to  smooth  AP-effects. 

XP. — Within  each  of  the  three  cooling  times,  the  negative  AP-values  were  read  from 
figure  85  for  each  of  the  five  temperature  curves.  The  absolute  values  of  these  were 
plotted  and  smoothed  over  P  within  CP-groups,  with  the  results  shown  in  figure  86  as 


YP2  CT  UI2 


0         10        20        30  40 


TEMPERATURE, T 

Figure  86.— |XP  |  smoothed. 


an  exponential  effect  in  the  reversed  P-scale,  (40-  P),  scaled  to  the  difference  between 
an  upper  intercept  (672),  and  scaled  to  the  absolute  difference  between  UI2  and  the 
curve  at  P=  zero.  Then,  multiplying  by  ( -  l)  to  return  |AP>|  to  its  original  negative 
state, 


XP=(-l)(UT2-(YP2/{40)"*)(40-  P)"2). 


1)12. — The  672-values  were  plotted  over  CPand  smoothed  as  shown  in  figure  87. 


COOLING  TIME.CT 

Figure  87.— C//2  smoothed. 


This  was  described  as  an  exponential  in  the  reversed  CP  scale,  (60-  CP),  subtracted 
from  the  upper  intercept  of  88.0,  and  scaled  to  the  drop  (7.5)  at  CP=  15.  The  resulting 
form  was 


UI2  =  88  -  (7.5/(45)2  7)(60  -  CP)2-7 
=  88  -  0.0002579  (60-  CP)2  7, 


YP2. — Similarly,  the  FP2-values  were  plotted  over  CPand  smoothed  (fig.  88). 
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lOr 


YP2    9  _ 


COOLING  TIME,CT 

Figure  88.—  YP2  smoothed. 


The  resulting  curve  was  described  as  an  exponential  in  the  reversed  CT scale,  (60  -  CT), 
scaled  to  the  difference,  l.O,  from  the  intercept  at  CT=  15,  and  added  to  the  intercept. 
In  this  case, 


YP2= 8.4  +  (1 .0/(45)3-97)(60  -  Cr)3  97 
=  8.4  +  2.734  x  10"  7(60  -  CT)3-91. 


The  shapes  of  the  curves  over  T  in  figure  86  were  determined  through  overlay  curves 
and  the  ^-Standards  to  be  virtually  constant  at  /?2  =  2.35. 
Then 


XP  =  (-  l)(t//2-(KP2/(40)235)(40-  T)2  35). 

UIl. — Following  the  analytical  pattern  for  XT,  peak  values  (UIl)  of  the  15  curves  in 
figure  85  were  read,  plotted,  and  smoothed  over  T  within  each  of  three  Cr-groups. 
Comparison  of  overlay  curves  with  the  A^-Standards  led  to  adoption  of  a  constant  flat 
form  (n  =  1.1)  for  all  three  curves.  Intercepts  and  scaling  heights  changed  over  CT  so 
that 

UIl  =  INT3  +  ( 7P3/(40) 1  • 1 )( T ) 1  1 

where 


INT3  =  12.1  -  0.003475(60  -  Cr)1-63 
YP3  =39.2  -  0.0004346(60 -Cr)2-45. 


ny — After  graphically  adjusting  the  figure  85  curves  over  WP  for  smoothed  values  of 
XP  and  UIl,  absolute  differences  between  the  UIl  and  the  associated  curves  were  read 
for  seven  points  in  WP  (XP<  WP<zero).  The  resulting  control  points  were  normalized, 
plotted,  and  compared  to  the  ^"-Standards.  From  this,     varied  sigmoidally  over  both 
T  and  CT  such  that 


(40  -T) 


40 


-  1 


a?,  =  1.18  +  0.185 


l-I 


l-I 


l-e 


l-I 
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where 


/=  0.45  +  0.000333(r)  +  (1.9846  x  10~  19)(Cr)10 


m  =  2.5  +  e 


(60  -CT) 
60 


-  1 


0.405 


So,  with  the  foregoing  inputs  and  as  originally  formulated, 


MK=  UI\  -  ( UIl / \XP\" i ) \XP -  WP\"\ 


can  now  be  estimated  with  the  foregoing  inputs.  Performance  of  the  model  is  shown  on 
figure  85  where  it  may  be  seen  that  the  model  is  an  excellent  reproduction  of  the  graphed 
hypothesis.  Modest  smoothing  effects  over  CT  exist  near  the  peaks  at  30  and  60  seconds. 
The  resulting  four  dimensional  model  is  presented  in  figure  89  at  two  points  in  the  fourth 
dimension  (CT=  15,60). 


-80       -60     -40       -20  O 

WATER  POTENTIAL, BARS 

Figure  89.— The  mathematical  hypothesis  for  micro- 
volts. 


Solution  for  WP. — As  stated  previously,  interest  here  focused  on  WP=f(CT,  T,  and 
MV).  Solving  the  MV  estimator  for  WP,  we  have 


WP  = 


\XPp(UI\-MV) 


UIl 


+  XP 


k(-  1),  XP<  WP<0. 
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Concave-Upward 
.Yl-Effects,  with 
Constant  Point  of 
Peaking  in  X2  and 
Interacting  with 
Convex-Upward 
^-Effects;  Fourth 
Dimension  Discrete 


Live  Data  Example  Here,  phloem  thickness  (P )  in  lodgepole  pine  was  expressed  as  a  function  of  electrical 

resistance  of  cambium  (R  ),  d.b.h.  (£>)  and  of  three  discrete  points  in  time,  May,  July, 
and  October.  The  hypothesis  (fig.  90)  was  based  on  a  fairly  large  data  set  from  a  period 
of  active  tree  growth  (July)  and  was  developed  within  the  constraints  of  prior  knowledge 
(data  from  D.  M.  Cole  and  Jensen  1980). 


0  10  20  30  40  50 


RESISTANCE  (kfl) 

Figure  90.— Phloem  thickness  of  lodgepole  pine  in  relation  to 
d.b.h.  and  electrical  resistance  of  cambium,  the  mathematical 
hypothesis. 


Data  for  May  and  October  were  subsequently  obtained  and  appeared  to  be  fairly  well 
represented  by  the  July  surface.  When  the  hypothesis  was  fitted  to  each  of  the  new  data 
sets  independently,  the  resulting  surfaces  bracketed  that  for  July  as  expected.  The 
hypothesis  was  subsequently  tested  with  new,  small  data  sets. 

Explanations  of  descriptor  development  that  follow,  assume  familiarity  with  2-D 
descriptor  methods  discussed  earlier. 

The  July  data  were  sorted  into  d.b.h. -groups.  Within  each  d.b.h. -group,  smoothed 
trends  over  R  were  estimated  in  accord  with  expectation  (shown  schematically  in  figure 
91). 

The  negative  relation  of  P  to  R  was  strongly  displayed  in  the  data,  so  this  was  the 
starting  point  for  the  descriptor  using  the  reversed  /?-scale  on  an  estimated  intercept  of 
0.01.  Then 


P  =  0.01  +  (YP/(45)")(50 -R)n. 
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+ 


Figure  91.— The  schematic  trends. 

Screening  the  resulting  curves  for  n  in  the  ^"-Standards,  n  varied  sigmoidally  from 
1 .9  at  smaller  diameters  to  about  3.2  at  larger  ones.  The  trend  was  estimated  to  be  as 
shown  in  figure  92  where, 

3.5  r- 
3.2 


Figure  92.— n  smoothed. 


n=\.9  +  (l.3)e 


D_ 
22 


0.7 


The  R -effect  was  scaled  at  R  =  5  or  (50 -R)  =  45,  since  this  was  the  lower  limit  of  the 
data  over  R.  The  height  of  the  July  surface  at  this  point  minus  the  intercept  of  0.01  =  YP. 
YP  varied  positively  (convex  upward)  with  D  as  shown  on  the  July  surface  (fig.  90)  at 
R  =  5.  This  effect  (fig.  93)  was  described  as  a  double  component  A^-function  of  the  re- 
0.16,- 

(0.148) 
1«t  COMPONENT 


0.12 

YP 

@ 

0.08 

0.04 

0 

0.65 


Figure  93.—  YP  smoothed. 

versed  D-scale,  subtracted  from  the  uppermost  YP,  (yp  =  0.158-0.010  =  0.l48).  The  first 
component  was  scaled  at  Z)  =  4.99,  the  second  at  0.65.  And 

YP=  0. 148  -  (2.2687  x  10"  5)(22  -  D  )2  72 
-(1.2694X  10-l4)(22-D)9, 

thus  completing  the  inputs  for  the  hypothesis  pictured  in  figure  90.  Limits  of  application 
are:  0<£>22;  0</?<50.  For  D>22,  n  =  3.2  and  for  R>50,  P  =  0.0l. 

Based  on  performance  of  the  hypothesis  for  the  May  and  October  data  sets,  R2  =  0.61 
and  0.78,  svx~0.22  and  0.018,  respectively,  the  July  hypothesis  appears  to  represent  other 
months  within  reason. 
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MATHEMATICAL  DESCRIPTION  OF  GRAPHED 
HYPOTHESES  -  FIVE  DIMENSIONS 
Outline  for  "FTVE  DIMENSIONS" 

PAGE 

Extension  of  the  4-D  Psychrometer  Example  (Pages  86-89)  to  Include 

the  Linear  Interactive  Effects  of  X3  and  X4   94 
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PAGE 


Symmetrical  Bell-Shaped  Al-Effect  with  Moving  Peak  in  A3;  X2,  a 
Sigmoid,  Interacts  with  XI  to  Form  the  Ridge  Shown,  Which  Changes 
in  Height  Sigmoidally  Over  the  Ranges  of  A3  and  X4  (Live  Data 
Example)    97 


Symmetrical  Bell-Shaped  Al-Effect  with  Constant  Point  of  Peaking 

in  XI.  The  Interacting  Effect  of  XI  is  Concave-Upward  and  Those 

of  A3  and  X4  are  Sigmoidal  (Live  Data  Example)   101 
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Extension  of  the 
4-D  Psychrometer 
Example 
(Pages  86-89)  to 
Include  the  Linear 
Interactive  Effects  of 
X3  and  X4 

Live  Data  Example  The  Brown  and  Bartos  (1981)  four-dimensional  thermocouple  example  (fig.  89  and 

94),  ultimately  had  a  fifth  dimension.  It  consisted  of  the  variable,  zero-offset  microvolts 


-80   -60  -40  -20  0 

WATER  POTENTIAL, BARS 

Figure  94.— The  4-D  microvolt  mathematical 
hypothesis. 


(ZOM),  and  is  an  index  of  the  temperature  differential  between  the  thermocouple  sen- 
sor tips.  Such  differential  has  an  effect,  zero-offset  error  (ZOE),  on  the  microvolt  out- 
put that  changes  depending  on  the  average  temperature  level  of  the  measured  medium 
and  on  the  direction  of  the  temperature  difference  between  the  tips,  as  indicated  in 
figure  95. 


ZERO- OFFSET,  MICROVOLTS 

Figure  95.— Zero-offset  error  source  information. 
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ZOE  is  quantified  in  figure  95  only  at  CT=  15  seconds  and  WP-  -22.5.  A  propor- 
tional ratio  of  change,  C,  for  microvolts  estimated  in  the  four-dimensional  model  (MV4) 
is  derived  from  AfK-readings  at  WP  =  -22.5,  as  specified  below: 

MVS  +  (MVs/MV]5)ZOE 

where 

MVs  =  MVat  specified  average  temperature,  ZOM,  and  \S<CT<60 

MVX$- as  above,  but  with  CT=  15. 
Then, 

MVS  =  MV4(C) 
where 

MV4=J{CT,T,  WP),  as  previously  specified,  and 
ZOE^b^ZOM)  +  b2(ZOM)(T) 
where 

b{  and  b2  are  estimable  directly  from  figure  96,  an  annotated  version  of  the  right  side 
of  the  three-dimensional  form  implied  in  figure  95. 


0  60 
ZOM 

Figure  96.— Zero-off  set  error  in  3-D. 
From  figure  96, 

6,  =  (0.9/60)(ZOM)  =  0.015(ZOM) 

b2  =  (3.4/60  x  40))(ZOM)(T)  =  0.001471(ZOM)(D. 


Note  that  the  formula  for  the  right  side  of  figure  95  applies  to  the  left  side  also,  since 

the  sign  of  ZOM  and  Cboth  become  negative  there. 

This  completes  the  five-dimensional  relation,  where  MVS  =  MV4(C)  and 

MV5  =  (un-{un/\xp\"i)\xp-  wpp)c. 

Again,  with  interest  focused  on  WPand  solving  MVS  for  WP,  we  have 


WP 


\XP\n\(UI\xC-MV) 
(7/1  xC 


1  n 


+  XP 


*(-  1),  XP<  WP<0. 
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Discrete-Variable 
Extension  of  the 
4-D  Terpene  Model 

Live  Data  Example  Given  the  4-D  relation  already  established  for  total  terpenes  (fig.  82  and  97),  the 

residual  problem  of  estimating  such  a  relation  for  each  of  three  component  or  compo- 
nent groups  of  the  total  (beta-phelandrene,  beta-pinene  +  3-carene  +  myrcene,  and  alpha- 
pinene)  was  approached.  The  associated  data  were  screened  graphically  over  D,  P,  and 
G,  and  were  found  to  contain  more  relative  variability  than  existed  with  analysis  of  total 
terpenes.  There  did  not  appear  to  be,  however,  any  recognizable  departure  of  basic 
forms  from  those  estimated  to  exist  for  total  terpenes. 

5-YR  RADIAL 


D.B.H.  ,IN  . 

Figure  97.— The  4-D  mathematical  hypothesis. 

Then,  the  obvious  alternative  of  fitting  the  relation  developed  for  total  terpenes  to 
that  for  each  component  or  component  group  by  least  squares  was  adopted.  That  is, 
where  estimated  total  terpene  (T)  =  X and  measured  component  terpene  =  Y,  the  least 
squares  coefficient,  b,  was 

b  =  ZXY/1LX2. 

The  expected  sum  of  the  three  coefficients  independently  obtained  was  the  b  of  1 . 10301 
determined  for  total  terpenes  (from  the  4-D  terpene  example).  So,  the  three  b  values 
were  adjusted  by  the  simple  ratio  of 

b  for  total  terpenes 
Lbl Tor  components  . 

Final  component  coefficients  (and  associated  R2,s)  are  shown  in  table  29.  Holding  G  at 
its  average,  1 .59,  the  resulting  surfaces  appear  as  in  figure  98. 

These  surfaces  would  vary  jointly  as  a  rather  flat  concave-upward  curve  over  G  since, 
from  the  4-D  terpene  example, 

YPD  =  0.38+  1.0292(C)14. 


96 


MATHEMATICAL  DESCRIPTION  OF  GRAPHED  HYPOTHESES  -  FIVE  DIMENSIONS 

Table  29.— 4-D  coefficients  and  f?2's  for  terpene  components 


b  R2 


Beta-phelandrene                                          0.64292  0.402 

Beta-pinene  + 
3-carene  + 

myrcene  30222  .352 

Alpha-pinene  15787  .213 


Total  terpenes  (sum)   1.10301  0.387 


PERCENT  OF  ALL  TERPENES 

i  ■  1 


D.B.H..IN. 


Figure  98.— Discrete  partitions  for  the  4-D  terpene 
hypothesis. 

Also,  since  the  form  of  the  relation  is  identical  for  each  component  terpene  and  the 
total  and  since  component  scales  sum  to  that  for  the  total,  components  will  always  sum 
to  the  total  at  any  point  in  the  relation. 

Symmetrical  Bell- 
Shaped  Al-Effect  with 
Moving  Peak  in  X3; 
X2,  a  Sigmoid,  Inter- 
acts with  XI  to  Form 
the  Ridge  Shown, 
Which  Changes  in 
Height  Sigmoidally 
Over  the  Ranges  of 
X3  and  X4 

Live  Data  Example  Our  hypothesis  development  here,  involves  a  small  data  set  (from  Lee  and  others 

1973)  associated  with  the  gap-graded  road  materials  problem  discussed  under  3-D,  asym- 
metrical bell-shaped  ridges.  These  data  were  exhausted  of  information  graphically  under 
constraints  of  expectation  to  arrive  at  the  mathematical  hypothesis  shown  in  figure  99. 


MATHEMATICAL  DESCRIPTION  OF  GRAPHED  HYPOTHESES  -  FIVE  DIMENSIONS 


>  > 

a  a! 


o  o 


a  a 

5  o 


>  > 

a.  o. 

II  II 

o  o 


Q  O 

5  o 


o 

Q. 


z  x  P 

U  u  = 

o  o  2 

c  c  7 

Ul  UJ  (0 

0-0.  £ 

II  II  ID 


The  explanations  of  descriptor  development  that  follow  presume  familiarity  with  the 
methods  presented  in  the  TWO  DIMENSIONS  section. 

The  effect  of  stability  (S)  on  the  dependent  variable,  ratio  (R ),  was  strongly  ex- 
pressed in  the  data  and  so  was  the  first  effect  specified  in  the  descriptor.  S  appeared  to 
have  a  consistent  sigmoidal  shape  within  partitioned  data  groupings  (some  missing 
because  of  the  small  data  set)  of  the  remaining  three  independent  variables,  flow  (F), 
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percent  voids  ( V),  and  percent  voids  filled  ( VF).  This  is  shown  schematically  in  the  top 
two  graphs  in  figure  100.  The  starting  point  in  the  descriptor  is 


35 


5000 


RATE  =  (yf)e 


0.9 


a  sigmoid  upward  to  the  right  over  S  with  zero  intercept  whose  vertical  scale,  YP,  is  de- 
pendent on  the  remaining  three  variables. 
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Figure  100.— Sequential  development  of  forms  and  scalars. 


MATHEMATICAL  DESCRIPTION  OF  GRAPHED  HYPOTHESES  -  FIVE  DIMENSIONS 


The  scalar  effect  of  F  at  5=  5000  is  bell  shaped  and  symmetrical  as  may  be  seen  in  the 
two  graphs  (second  pair  from  the  top)  in  figure  100.  Adding  20  to  Fto  accommodate  the 
broadest  curve  of  the  set  under  the  constraint  that  e~K  functions  only  between  zero  and 
2XP,  these  curves  are  described  by 


(F+  20) 
PXP 


YP  =  FP 


l-FI 


l-e 


1.6 


1 


1.6 


FI 


1.6 


l-FI 


where 

FP,  the  scalar  for  YP,  is  estimated  (conservatively)  to  vary  sigmoidally  with  both  V  and 
VF,  as  shown  in  the  third  tier  of  graphs  from  the  top  of  figure  100.  The  descriptor  is 


FP=FPF+(FPP)e 


10 


0.78 


50 


and 


FPF  =  2.4  +  (\. 8)  e 


(100-  VF) 
40 
0.42 


-  1 


FPP  =  0.8  + (2.9)  e 


(100-  VF) 
40 
0.42 


where  FPP  adds  (3.2- 2.4  =  )0.8  at  VF=  100,  and  (7.9- (2.4+  1.8  +  0.8)  =  )2.9  at  VF=0. 
PXP  is  mobile  and  is  assumed  to  be  sigmoidal  as  shown  in  the  graph  at  the  bottom  of 
figure  100.  PXP  changes  with  VF  only  and  can  be  identified  as 


PXP=  31  +  (2.8)(? 


(100-  VF) 
40 


-  1 


0.42 


And  the  inflection  points  (FI)  for  the  bell-shaped  yP-curves  vary  with  V,  as  shown  in  the 
left  graph  of  the  two  just  above  the  one  for  PXP.  The  graph  to  the  right  of  that  shows 
associated  sigmoidal  trends  (assumed)  over  KFfor  intercepts  (FIF)  and  peaks  (FIP).  All 
told,  the  descriptor  for  FI  is  as  follows 


FI=FIF+(FIP)e 


(10-  V) 
10 


-  1 


0.27 


100 
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where 


f7F=0.41  +  (0.15)e 


(100-  VF) 
40 
0.42 


FIP=(0.08)e 


(100-  VF) 
40 
0.42 


This  chain  of  descriptor  components  is  converted  to  the  Fortran  IV  statements  shown 
below: 


T=(EXP(-(ABS(((100-  VF)/40- l)/.42)**9))) 

PXP=3l+2.S*T 

FPF=  2.4+  1.8*7 

FPP  =  .8  +  2.9*  T 

FIF=A\  +  .\5*T 

FIP=M*T 

FP  =  FPF+  FPP*(EXP(  -  (ABS((  V/\0-  l)/.78)**50))) 
FI=  FIF+  FIP*(EXP(- (ABS(((10-  V)/10-  1)/.27)**10))) 
LN  =  EXP  (-(ABS(((F+  20)/ PXP-  1)/(1  -F/))**1.6)) 
RN  =  EXP  (-((1/(1  —FT))**  1. 6)) 
YP  =  FP*((LN-  RN)/(\-  RN)) 

RATE=  IM15*YP*(EXP(-(ABS((S/5000- l)/.9)**35))) 

Note  that  the  coefficient  of  1.0315  resulted  from  the  refitting  of  the  mathematical 
model  back  through  the  data  from  which  it  was  derived  by  least  squares.  Given  that 
modeled  rates  =  X and  actual  rates  =  Y,  b  =  LXY/LX2.  Then  bX (or  RATE)  is  the 
model  submitted  to  electronic  computer/plotter  processing  that  resulted  in  the  graphic 
presentation  of  the  mathematical  hypothesis  in  figure  99. 

It  is  of  interest  that  the  additive  model  with  linear  effects  of  S,  F,  V,  and  VF,  reached 
an  R2  of  only  0.42  under  least  squares  fit,  while,  for  the  model  just  derived,  R2  =  0.80. 
There  would  appear  to  be  substantial  benefit  from  utilizing  all  information  at  hand  for 
developing  the  hypothesis  in  this  case. 

Symmetrical 
Bell-Shaped 
*1-Effect  with 
Constant  Point  of 
Peaking  in  XI. 
The  Interacting 
Effect  of  XI  is 
Concave-Upward  and 
Those  of  X3  and  X4 
are  Sigmoidal 


Live  Data  Example  The  flexibility  of  e  K  in  representing  a  complex  relation  is  again  evident  in  the  five- 

dimensional  interaction  pictured  in  figure  101. 
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Here,  data  from  Lyon  1976  and  Lyon  and  Jensen  1980  were  the  basis  for  an  index  to 
intensity  of  deer  use  of  forest  openings  created  by  clearcutting  in  western  Montana.  The 
index  is  expressed  as  a  function  of  opening  size,  height  of  new  vegetation,  depth  of 
slash  in  the  opening,  and  depth  of  dead  and  down  timber  adjacent  to  the  opening. 

The  data  at  hand  were  initially  committed  to  statistical  evaluation  of  the  linear  effects 
of  the  above  variables  and  other  independent  variables  on  the  intensity  index.  These 
were  virtually  the  simplest  regression  hypotheses  that  could  have  been  adopted.  In  this 
case,  the  evaluation  provided  only  weak  support  for  expected  results  and  yielded  little 
new  information. 

To  exploit  the  data  further,  a  hypothesis  was  developed  for  the  set  of  four  independ- 
ent variables  above.  Prior  knowledge  of  the  forms  of  the  relations,  including  inter- 
actions, between  these  variables  was  summarized.  Subject  to  these  constraints,  the 
limited  data  were  then  partitioned  and  graphically  exhausted  of  associated  curve-form 
and  scaling  information.  Data  appeared  to  provide  strong  support  for  the  dynamic 
interaction  anticipated.  Procedures  specified  here  were  used  to  develop  a  functional 
form,  XT,  for  the  graphed  interaction.  This  was  refitted  to  the  data  by  least  squares  in 
the  model  PGR  =  b(XT)  and  b  =  0.9721 .  The  adjusted  model  is  specified  below: 


PGR  =  J{  VI,  SI,  SO,  Acres) 


where, 

PGR  =  number  of  deer  pellet  groups  per  acre  inside  the  opening. 

VI  -  height  of  vegetation  in  feet,  inside  the  opening. 

5/  =  depth  of  logging  slash  in  feet,  inside  the  opening. 
SO  -  depth  of  dead  and  down  timber,  in  feet,  outside  the  opening. 
Acres  =  size  of  opening  in  acres. 


USKPO 


PGR  = 


YP 


(PO-0.6)1 55 


{PO-S/}155  (0.9721)(50) 


where 


YP=  YPS 


Acres  +  200 
260 
I -I 


1.75 


1.75 


1.75 


1-/ 


+  0.23 


\-e 


\-I 


VI 


-1 


2(1 


3 -SO 


-  1 


FPS  =  0.57  +  3.23 


0.76 


0.4942 


YL 
10 


10.5 


PO  =  6.289e 


0.999 


-2.289 
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VI 


- 1 


3.5 


/=  0.897  -  0.061  e 


0.404 


USI>PO 
PGR  =  0 
Limits 

0<S/<4,  0<5O<3,  0<  K/<8,  0<  Acres<300. 


Graphic  development  of  the  model  described  by  this  function  apparently  resulted  in 
great  sensitivity  to  the  interaction  information  contained  in  the  data.  But,  since  the  de- 
grees of  freedom  thereby  sacrificed  were  unknown,  conventional  statistical  parameters 
were  not  estimable.  Some  indication  as  to  the  goodness-of-fit  of  the  functional  form  to 
the  data  set  from  which  it  was  largely  derived  is  provided  by  the  proportion  of  the  total 
sum  of  squares  explained  by  the  model,  R2-0JI  here. 

Contrast  this  with  the  R2  of  0.21  achieved  with  a  minimum-effort  additive  regression 
model  where  linear  ef  fects  of  the  same  four  independent  variables  have  been  fitted  to 
the  data  by  least  squares.  It  would  appear  that  sharp  focus  on  interaction  formulation 
was  justified. 

At  this  point  in  the  paper,  it  is  assumed  that  only  a  minimum  of  explanatory  informa- 
tion will  be  necessary  to  clarify  descriptor  development. 

The  expected  negative  SI  effect  was  strongly  apparent  in  the  data  and  so  was  used  as 
the  starting  point  for  the  descriptor.  This  effect  was  exponential  over  the  reversed  SI- 
axis  where  0.6<SI<PO.  The  lower  extreme  of  the  data  in  SI,  0.6,  was  adopted  as  the 
scaling  point  (largest  PGR-va\ues)  in  the  reversed  S/-scale.  PO,  a  variable  over  VI, 
represented  the  point  in  SI  where  the  PGR-curves  reached  zero  (fig.  102). 


PGR 


2  3 
1.0  0 

1.6  o.eo 

2.5   1.5  0.50 

Figure  102.— The  S/-effect  with  changing 
points  of  origin,  PO. 


PO 

SI 

(PO-3l)@PO  =  3.0 
(PO-SI)  @  PO-3.5 
(PO-SI)  @  PO-4.5 


When  the  point  of  exponential  origin  was  allowed  to  vary,  n  in  X"  approached  a  con- 
stant (1 .55)  over  SI,  but  the  scaling  heights  of  these  curves  (@S/=  0.6)  varied  over 
ACRES,  VI and  SO. 
Then, 


PGR  =  (50xb )( YP/(PO  -  0.6)1  55)(PO  -  SI)155,  S/<  PO 
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where 

50  is  the  constant  necessary  to  the  conversion  of  data  units  to  acres; 

b  is  the  least  squares  coefficient  fitting  the  derived  model  back  through  the  source 

data. 

Note  that  segmentation  is  employed  in  the  exponential  variable  here,  since  beyond  PO 
in  SI,  PGR  is  specified  as  zero.  PO  varied  sigmoidally  over  K/and 


PO  =  6.2S9e 


10  

0.999 


10.5 


2.289 


YP=  0.23  +  YPS 


Acres +  200 
260 
1-/ 


l-e 


1.75 


-  1 


1.75 


-I 


1.75 


1-/ 


XP  appeared  constant  at  260  acres  and,  from  the  e~A-Standards,  /?=  1.75  was  reasonably 
representative  for  all  curves.  /,  however,  varied  sigmoidally  over  VI  and 


YL  _  i  3  5 


/=  0.897  -  0.067  e 


0.404 


YPS,  the  heights  of  the  bell-shaped  curves  over  ACRES,  varied  as  an  interaction  of 
sigmoids,  a  new  descriptor  feature.  YPS  varied  sigmoidally  over  VI  as 


VI 


20 


KPS  =  0.57 +  3.23e 


0.76 


but  only  achieved  this  magnitude  of  effect  at  smaller  SO-values.  The  ^/-effect  became  sig- 
moidally smaller  as  SO  increased.  Noting  that  the  unsealed  sigmoid  over  the  reversed  SO- 
axis  (3- SO)  ranges  in  value  from  zero  to  one,  the  SO-effect, 


3-SO 


0.492 


10 


becomes  a  convenient  multiplier  with  which  to  achieve  the  appropriate  sigmoidal  reduction 
of  the  maximum  ^7-effect. 

All  curve-form  estimates  were  derived  from  data-based  graphic  hypotheses,  normalized 
control  points,  overlay  curves,  and  comparison  to  either  the  X"-  or  e_A:-Standards.  While 
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the  foregoing  mathematical  specifications  may  seem  voluminous,  they  are  quite  easily 
specified  in  a  short  Fortran  IV  program  for  computer  processing,  as  shown: 

PO  =  6.2S9*(EXP(-(ABS((VI/\0-  1)/.999)**10.5)))  -  2.289 
IF(SI.GT.PO  )GO  TO  1 

A  =  EXP(-  (ABS  (((3  -  SO  )/3  -  1)/.492)**10)) 
B  =  EXP(-(ABS((VI/S-  l)/.76)**20)) 
TPS  =  3.23M*£+.  57 

1=  .897-  .067*{EXP(-(ABS((VI/&- 1)/.404)**3.5))) 

LN=EXP(-(ABS(((ACR  +  2O0)/2Q0-  1)/(1  -  7))**1.75)) 

RN=EXP( -((1/(1  -I))**  1.75)) 

YP  =  YPS*((LN- RN)/(l  -  RN))  +  .23 

PGR  =  (YP/{PO- . 6)** L55)*((PO-S/)**1.55)*(0.9721)*50 

RETURN 
1^(5)  =  0 
RETURN 
END 

Note  that  figure  102  was  plotted  electronically  using  a  program  that  incorporated  the 
above  statements. 
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APPENDIX-A:  eK  DERIVATION 

e~^l  Its  Derivation,         The  new  function,  simply  identified  as  e~K,  provides  the  analyst  with  a  finite  source 
Potentials,  and  °f  versatile  transformations  for  use  when  it  is  inefficient  to  search  for  alternatives  or 

Limitations  when  alternatives  are  inadequate.  While  e~ K  is  not  a  panacea  for  all  the  problems 

modelers  face  in  developing  functional  relationships,  it  along  with  curves  of  the  class 
X"(n  >  1)  does  serve  a  broad  spectrum  of  transformation  needs  with  no  particular 
limitations  as  to  curve  shapes  for  which  it  is  most  useful.  Methods  for  developing  math- 
ematical descriptors  using  these  functions  have  been  treated  by  Jensen  and  Homeyer 
(1970,  1971)  and  Jensen  (1973,  1976,  and  1979)  and  are  summarized  and  updated  in  this 
paper. 

A  variant  of  the  Normal  function  e~  K  is  defined  as: 


(X/XP)  -1 

n 

-  e 

1 

n 

(1-/) 

Id-/) 

l-e 

1  ] 

n 

(1-/) 

,  0<X<2Xp 


where, 

e  =  natural  logarithm  base. 
XP  =  pivot  point  in  X  for  e~K  and  is  the  point  in  X  at  which  the  bell-shaped  curve 
peaks. 

XI  =  point  of  sigmoidal  inflection  in  X. 
I  =  proportional  departure  of  XI  from  zero  to  XP  or  from  2XP  to  XP. 
n  =  power  of  negative  exponents  for  e. 

The  left  numerator,  like  the  Normal  function,  generates  a  system  of  bell-shaped 
curves  about  a  central  point  in  X,  XP  in  this  case.  These  curves  reach  the  maximum 
value  of  1 .0  at  XP  and  then  decline  sigmoidally  and  symmetrically  on  either  side  of  XP 
with  increasing  departure  of  X  from  XP. 

User  control  of  the  curve  system  is  enhanced  if  all  curves  of  the  set  range  in  value 
from  zero  to  one  within  a  finite  domain  of  X.  To  achieve  this  property,  curves 
generated  by  the  left  numerator  have  been  forced  through  zero  at  XP±XP  as  follows: 

Let  the  left  numerator  =  e~  T  and  let  e~  r°  =  e~  T  at  X=  0  or  2XP.  Consider  the 
residuals  generated  by  1  —  e~  T.  Expanded  by  the  inverse  (1  -  e~  r°)~ 1  and  subtracted 
from  one,  we  have 

e~K=l-(l-e~T^~l(l-e~T)1 

which  simplifies  to  the  final  form 

e-K=(e-T-e-To)(l-e:T°y]. 

Then  the  right  numerator  and  the  denominator  serve  to  force  curves  of  the  left 
numerator  through  zero  at  XP±  XP.  This  domain,  along  with  XP,  can  be  altered  to  ac- 
commodate the  needs  of  individual  conceptual  models  by  adding  constants  ( +  or  - )  to 
the  A'-scale.  The  apparent  complexity  of  e~K  is  much  reduced  by  the  fact  that  in  appli- 
cation, the  right  numerator  and  the  denominator  simplify  to  constants  or  approach  zero 
and  can  be  deleted. 
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Sigmoids  on  either  side  of  XP  are  forced  into  different  areas  of  two-dimensional 
space  by  shifting  /,  the  proportional  point  of  inflection  in  X,  and  the  slope  of  the 
sigmoidal  face  at  /  changes  with  n. 

Arrays  of  curves  from  the  left  half  of  this  function  (a  mirror  image  of  the  right  half) 
are  shown  in  the  e~  A-Standards  that  follow.  Sets  are  given  for  commonly  used  in- 
crements of  n,  \  .2<n<  10.0.  Each  set  has  curves  progressing  from  7=0.1  at  the  left  to 
0.9  at  the  right.  These  curves  or  any  portions  thereof  provide  an  almost  endless  poten- 
tial for  matching  graphed  curves  and  describing  them  mathematically.  Given  that  YP  is 
the  peak  value  of  y(the  dependent  variable)  on  the  objective  curve,  the  selected  e~K 
-function  (e~K  ranges  from  zero  to  one)  may  be  scaled  to  that  curve  through  multi- 
plication by  YP. 

While  the  controllable  parameters  /,  n,  XP,  and  YP  provide  exceptional  descriptive 
power  in  two  dimensions,  such  power  is  magnified  manyfold  for  hypotheses  involving 
three  or  more  dimensions,  as  well  as  strong  interactions.  For  example,  the  location  of 
peaks  Fthat  move  across  the  independent  variables  in  any  fashion  can  be  described 
with  XP=  /,(ArI).  The  heights  of  these  peaks  in  any  pattern  and  changing  sigmoidal 
shapes  leading  up  to  them  can  again  be  pictured  as  functions  of  the  independent 
variables  YP=  /2(^,),  n  =  /3(Ar/),  and  / = fj^X).  Additionally,  asymmetry  in  bell-shaped 
hypotheses  is  easily  accounted  for  by  describing  sigmoids  on  either  side  of  XP  in- 
dependently (segmented  descriptors),  while  scaling  both  to  the  same  height  at  XP. 

A  descriptor  (XT),  adopted  as  a  hypothesis  in  its  entirety,  may  be  fitted  to  pertinent 
data  by  least  squares  in  the  simple  model,  Y=  (3XT+  e,  where  the  e  are  NID  (0,  a2), 
constant  variance,  and  (3  =  LXTY/LXT2. 

Weighted  regression  procedures  are  recommended  where  the  variance  about  regression 
is  not  uniform  over  the  ranges  of  the  independent  variables.  In  such  cases,  reasonable 
success  in  achieving  constant  variance  has  been  obtained  by  solving  for  departures 
( Yj  -  Y,)2  or  d}  as  a  suitable  function  of  the  related  Y  in  d2  =  b  Y".  Then  the  weight,  w,  for  each 

v    l         I'  l  A  'A  A  A 

observation  is  set  equal  to  \/Yn  and  a  weighted  j3  in  Y=  QY is  estimated  as  (3  =  LwYY/LwY2. 

More  complex,  generally  iterative  fitting  procedures,  such  as  the  Newton-Raphson  methods, 
can  be  used  to  arrive  at  statistical  estimates  of  internal  model  parameters  (see  Damaerschalk 
and  Kozak  1977;  Draper  and  Smith  1966). 
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The  problem  of  over-simplified  regression  hypotheses  is  identified  and  alter- 
natives are  discussed.  Methods  for  developing  structured  hypotheses,  first 
graphically  and  then  mathematically,  are  shown  for  two-  through  five- 
dimensional  cases.  Mathematical  formulation  is  accomplished  through  use  of 
two  families  of  curves:  X",  n>1,  and  e~K,  a  highly  manipulate  function  for  bell- 
shaped  curves  or  portions  thereof. 
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The  Intermountain  Station,  headquartered  in  Ogden,  Utah,  is  one 
of  eight  regional  experiment  stations  charged  with  providing  scien- 
tific knowledge  to  help  resource  managers  meet  human  needs  and 
protect  forest  and  range  ecosystems. 

The  Intermountain  Station  includes  the  States  of  Montana, 
Idaho,  Utah,  Nevada,  and  western  Wyoming.  About  231  million 
acres,  or  85  percent,  of  the  land  area  in  the  Station  territory  are 
classified  as  forest  and  rangeland.  These  lands  include  grass- 
lands, deserts,  shrublands,  alpine  areas,  and  well-stocked  forests. 
They  supply  fiber  for  forest  industries;  minerals  for  energy  and  in- 
dustrial development;  and  water  for  domestic  and  industrial  con- 
sumption. They  also  provide  recreation  opportunities  for  millions 
of  visitors  each  year. 

Field  programs  and  research  work  units  of  the  Station  are  main- 
tained in: 

Boise,  Idaho 

Bozeman,  Montana  (in  cooperation  with  Montana  State 
University) 

Logan,  Utah  (in  cooperation  with  Utah  State  University) 

Missoula,  Montana  (in  cooperation  with  the  University 
of  Montana) 

Moscow,  Idaho  (in  cooperation  with  the  University  of 
Idaho) 

Provo,  Utah  (in  cooperation  with  Brigham  Young  Univer- 
sity) 

Reno,  Nevada  (in  cooperation  with  the  University  of 
Nevada) 


